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Foreword 


This volume contains the proceedings of the Workshop on Observations of Recent 
Comets, 1990, held on June 15 and 16, at the Doubletree Hotel in Albuquerque, New Mexico, 
with the support of the Planetary Atmospheres Program of NASA. There were about 35 
participants, including one from France. Two Japanese groups, who did not attend in person, 
submitted posters. 

One-page abstracts were submitted in advance. These were made available to all partic- 
ipants at the time of the Workshop. All presentations were informal and about ten minutes in 
length, followed by about ten minutes of discussions. A written record of most, although not 
all, of the discussions was made. The presenters were encouraged to submit an expanded ab- 
stract of about five pages in length for inclusion in the Proceedings, The recorded discussions 
follow each of these presentations. 

At the time of the Workshop about three years had passed since the highly successful 
and extensive spacecraft and Earth-based observing campaigns of Comet P/Halley. During 
this time, data have been reanalyzed and reinterpreted. In light of the new findings from the 
in siiu measurements and coordinated Earth-based observations and the somewhat earlier 
in situ measurements from the P/Giacobini-Zinner flyby of the ICE spacecraft, the objective 
of the Workshop was to focus on potential interpretations of the observations of recent comets. 
Since much of the data on the recent comets Brorsen-Metcalf (1989o), Okazaki-Levy- Rudenko 
(1989r), Aarseth-Brewington (1989al), and Austin (1989cl) was only partially reduced and 
not fully analyzed or interpreted, these extended abstracts have not been refereed. Consis- 
tent with this objective, it was felt that the Proceedings should get into the hands of the 
participants, as well as other colleagues who were unable to attend, without delay. This is 
another reason for not refereeing these extended abstracts. Thus the reader is warned of the 
preliminary nature of these Proceedings and is encouraged to contact the authors directly. 
We hope that such contacts will lead to profitable collaborations. 

The editors wish to acknowledge the interest of the Planetary Atmospheres Program of 
NASA in sponsoring the Workshop and the publication of the Proceedings. 

Walter F. Huebner 
Peter Wehinger 
Jurgen Rahe 
Ichisliiro Konno 
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RESULTS FROM THE COMET GIACOBINI-ZINNER AND HALLEY CAMPAIGNS 
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The spacecraft flybys of Comets Giacobini-Zinner and Halley were highly successful. 
They confirmed many untested theories and provided additional clues to decide between 
competing models, but without leading to a clear preference* They also raised new questions. 
Based on the improved knowledge, research efforts will now be much better focused than 
was possible before the encounters. Moreover, we should be better equipped to interpret 
Earth-based (ground-based, air-borne, rocket-borne, and Earth satellite) observations. 

The goals of the workshop are defined in terms of the results from spacecraft measure- 
ments and Earth-based observations of Comets Giacobini-Zinner and Halley. A summary is 
presented of the new knowledge about the nucleus, the coma gas and dust, and the plasma. 
We realize that we must be extremely cautious when generalizing results from the Comet 
Giacobini-Zinner and Halley investigations to comets universally. The many surprises of the 
successful missions, the dangers inherent in generalizations to other comets, and the many 
still open and new questions are the reasons for the need of further spacecraft investigations. 
However, an unquestionably successful start has been made in expanding our knowledge about 
the physics and chemistry of comets and all the related implications, and the question about 
the similarities of and differences between comets can now be explored further. 


1. Summary of the P/Giacobini-Zinner and P/Halley Investigations 


Starting with the P/Halley nucleus, we have learned that its size of approximately 500 
km 3 is about three times larger in volume than estimates before the P/Halley missions had 
indicated (Keller et al., 1987). The larger size, combined with the same response to the 
nongravitational forces acting on it (therefore the same mass of about 1.5 x 10 14 kg), leads 
to a much lower density of about 300 kg/m 3 (Rickman, 1989). With the uncertainties in 
the volume and mass, the density could be as low as 200 kg/m 3 and as high as 700 kg/m 3 . 
The axial ratio of the nucleus is about 2:1:1 with the major axis of the ellipsoidal shape 
being about 16 km; this is much more aspherical than had been anticipated. The surface 
morphology suggests subunits of the nucleus that may have been as large as one to several 
kilometers in size. Consistent with the large surface area of about 400 km 3 , the albedo is only 
about 4%. The abundance of CO and probably some of the other coma gases is much lower 
in the nucleus than earlier coma observations had suggested. The nucleus itself appears to 
be inhomogeneous, as indicated by the relatively small active areas that are responsible for 
the jet-like dust features and most, if not all, of the coma gas. The rotational motion of the 
nucleus has been controversial since the spacecraft flyby of P/Halley, but a possible solution 
may now be in sight. 

More relevant to ground-based observations are the coma and the dust. We have learned 
that the frozen gases in the nucleus of P/Halley consist mainly of water. About 85% of the 
coma gas near the nucleus is water. Approximately half to two-thirds of the CO comes 
from distributed sources in the coma, probably indirectly through the dissociation of complex 
molecules that are released from the dust. Closely related to the decay of these complex 


molecules may be the production of H 2 CO, CH, and CHj. Even other species, such as CN, 
NH, NH 2) C 2 , and possibly C3, may come, at least in part, from distributed sources, as the 
jetjrlike structures of these radicals in the coma suggest. The ratio of 13 C/ 13 C = 65 ± 9, as 
defWmined % Wyckoff et al*. (1989). This ratio, obtained from CN observations, is much 
lower than the older values of 89 and higher, which were obtained from Cj observations in 
comets (Vanysek and Rahe, 1978). The value of 89 is the terrestrial value. The carbon isotope 
ratio obtained from the P/Halley dust analysis varies from very low to very high values (Sole 
et al., 1987). It would be interesting to know whether the carbon isotope ratios from CN 
and C 2 are consistently different and if the ratios from gas-phase molecules are different from 
those of comet dust. The ratio of 0.6 X 10“ 4 < D/H < 4.8 x 10“ 4 obtained by Eberhardt et 
al. (1987) is very similar to that of ocean water and of Titan. Very heavy molecules with 
molecular weight up to about 120 amu were detected (Korth et al., 1987). It is likely that 
these molecules are associated with the organic dust known as CHON particles. Gas “jets” 
have been observed in CN, OH, NH, and other molecular emissions. These are most likely 
also related to CHON dust particles. The detection of para- and ortho-water in IR spectra 
(Mumma et al,, 1987) assumes a special importance: It has the potential to open a new era 
of investigations about the formation temperature of comet nuclei. 

The dust size distribution does not appear to have a lower limit; particles as small as 
10 -16 g (the lower limit of particle sizes measurable by the instruments) were detected in 
large numbers (McDonnell et al., 1987). Could this distribution of small particles continue 
downward to the size of molecular clusters? There is evidence that du3‘i particles fragment in 
the coma. Large amounts of organic dust, rich on carbon, hydrogen, oxygen, and nitrogen, 
the so called CHON particles have been detected (Clark et al., 1987). They also contain small 
amounts of sulfur. The CHON particles and their disintegration products are also responsible 
for the IR emissions in the 3.3 pm range. CHON particles can be divided into six categories: 
Those consisting primarily of (1) H and C, (2) H, C, and 0, (3) H, C, and N, (4) H, C, N, and 
0, (5) C and 0, and (6) H and 0. These divisions are probably not sharp, diffusing from one 
group to the next. Nitrogen is underabundant by a factor of 3 or perhaps somewhat more as 
determined from the combined analyses of gas and dust in P/Halley. 

The relative abundances of the rock-forming elements in the dust are within a factor 
of two of the Cl chondrites, while H, 0, and especially C and N are highly enriched, i.e., 
less depleted in comet dust relative to Cl chondrites. Thus the dust in P/Halley is more 
primitive than Cl chondrites, which testifies to the nearly pristine nature of comet dust. The 
distribution of the Fe/(Fe + Mg) ratio in P/Halley silicate dust resembles the distribution 
found in anhydrous interplanetary dust particles (IDPs). Comparing the ternary Mg-Si-Fe 
system, P/Halley dust differs greatly from layer lattice silicate IDPs and matches better 
with anhydrous IDPs (Jessberger et al,, 1988). But there is an important difference between 
comet dust and anhydrous IDPs: Anhydrous IDPs are practically devoid of samples richej; 
in Mg then those lying on the olivine-iron line, while in comet dust many silicon-poor and 
magnesium-rich grains occur. This is what one would expect from primitive Solar System 
material, i.e., silicates that never had been exposed to high temperatures. The dust-to-gas 
mass ratio is about 1 : 1 with an uncertainty of a factor of about 2. If the dust size distribution 
is extrapolated and integrated to the largest particles that can be lifted through entrainment 
by coma gas against the gravitational attraction of the nucleus, then the dust-to-gas mass 
ratio is about 2 in P/Halley. 

Of equal importance, but less accessible for Earth-based observations, are the plasma 
effects that have been determined by in situ spacecraft measurements. The cometopause, 
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which is defined by a particularly steep increase of mass loading of the solar wind, was 
discovered between the bow shock and the contact surface (Gringauz et al., 1986). It was 
found that the bow shock is not sharp and that there exists a turbulent transition zone, from 
the ICE passage through the onset of the ion tail region of Comet Giacobini-Zinner, it was 
found that: (1) The model for field-draping is basically confirmed, but a magnetosheath forms 
at the ion tail boundary. ( 2 ) A thin plasma sheet with a total thickness of about 2000 km and 
a width of about 1.6 x 10* km exists at the center of the ion tail. (3) The plasma flow velocity 
gradually decreases to zero at tl>tt center of the ion tail. (4) The degree of magnetic field 
turbulence and plasma wave activity was seen to correlate with the production of ions from 
cometary neutrals in the solar wind. The study of plasma instabilities, waves, pickup ions, 
ring and shell distributions in velocity space, energetic ions, and charge exchange processes 
with cometary neutrals is afforded by the natural plasma laboratory that a comet provides 
without the interfering wall effects that are always present in plasma chambers. 

2. Relationship to Earth-based Observations 

The above summary is obviously incomplete, but it provides a starting point for dis- 
cussions in this workshop. Among the questions that we may want to address is: How well 
can we define and interpret Earth-based observations in terms of the above results? Some 
observations of P/Halley, such as the gas “jets” (A’Hearn et al., 1986), were ground-based 
and can obviously be tried again. It would be very important to know how common these 
jet-like gas features are in other comets. Other observations will be much more difficult to 
make by remote sensing from the Earth. Here we summarize some of the questions that might 
be answered through such observations. 

Many attempts have been made to determine the spin period of comet nuclei, but how 
well can one determine the combination of spin and precession of the nucleus from ground- 
based observations? What fraction of a nucleus is active? How is this fraction related to the 
dynamical age of a comet? What causes the development of inactive (mantle) areas on the 
surface of a nucleus? What fraction of the coma gas comes directly from the nucleus (i.e., 
from the ice) and what fraction comes from the distributed (i.e., dust) sources? How does 
the relative abundance of gas to dust change and how does the relative abundance of frozen 
gases (e.g., CO : COj : H 2 O) change from comet to comet? Very important for the origin of 
comets is the question about the ratio of CH 4 /CO and NH 3 /N 2 . What are the parent species 
of C 2 , C 3 , CN, CH, NH, and NH 2 ? Are gases trapped in the ice of the nucleus? If so, Ar 
should also be trapped (Bar-Nun et al., 1985). The relative abundance of Ar may become 
an important test case for the trapping theory. Do clathrates exist in comets? How wide 
does the S 2 abundance vary from comet to comet? Are other sulfur compounds, such as H 2 S, 
SO, SO 2 present in comets? What is the , 2 C/ 13 C ratio in the ice and in the dust? Can this 
isotope ratio, as obtained from CN and C 2 , be used to determine the sources of CN and C 2 ? 
For example, does a narrow range in the isotope ratio indicate that the source is in the ice 
and a broad range indicate that dust is the source? What is the D/H ratio in other comets? 
What is the gas production of comets at large heliocentric distances (say, r > 4 AU) and 
which species can be detected? 

What is the ratio of silicate to CHON particles in dusty vs. “dust-free” comets? Do 
the silicate and CHON particle size distributions differ with respect to each other and as 
a function of cometocentric distance? What are the compositions and structures of CHON 
and mixed silicate-and-CHON particles? If “dust-poor” comets are rich on very fine dust 


(particle radius a < 0.1pm), then the production of some “mother molecules” , such as CO, 
HaCO, etc., may be detectable in distributed sources from the ground (see, e.g., Feldman et 
a!., 1988). Are jet-like gas features frequent phenomena in comets? Do jet-like features in OH 
(Clairemidi and Moreels, 1989; Clairemidi et al., 1990) originate from HaO trapped in dust 
particles? Can IR observations detect fragmentation of dust particles? For example, docs the 
relative intensity profile at different IR wavelengths change with cometocentric distance? Can 
we improve detection of distributed sources from Earth-based observations? Can unidentified 
spectral features be related to heavy molecules with masses larger than about 40 amu? What 
is the relationship between amorphous and crystallized silicates and how do they relate to 
other features in comets (e.g., dust-poor or dust-rich comets, the 3.3 pm C-H stretch bonds 
of the organic grains, etc.)? What is the relationship between C 2 , C3, CN, CO (CO + ), dust, 
etc,, in short-period vs. long-period comets? 

With regard to plasmas, can we detect any evidence for an inner shock? How stable is 
the contact surface? How important is the Alfven critical velocity effect (anomalous ioniza- 
tion)? Can one identify the cometopause from ground observations, and if so, how stable is 
it? What are the energetic ions inside of the contact surface and what is their source? What 
are the causes and dynamics of disconnection events (DEs) in plasma tails? We may need 
simultaneous observation of DEs and spectroscopy of tail ions. 

Many other questions can undoubtedly be asked. Some new questions and perhaps 
even some answers may result from discussions in this workshop. But most important for 
this workshop is to outline some of the thoughts and provide the stimulus for development of 
further research and collaboration. 
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Lewis Snyder 

What is the present status of the formaldehyde polymer identification? 

Walter Huebner 

Laboratory experiments with H 2 CO solutions in H 2 O give a mass spectrum in the vapor phase 
that closely resembles the observed POM peaks, H 2 CO in methanol solution produces the 
mass peaks in the vapor exactly and also shows some neighboring peaks to give a width to 
the spectral features, but the width does not fully account for the observed spectrum. Other 
species must be present. 


Susan Wyckoff 

Do the bulk properties of the Halley grains differ chemically from those analyzed by Anders 
and co-workers in primitive meteorites? 

Walter Huebner 

In layer lattice silicate interplanetary dust particles (IDPs) the ratio of Fe/(Fe + Mg) peaks 
sharply around 0.5 and is similar to a carbonaceous chondrite matrix. In anhydrous IDPs 
Fe/(Fe -f Mg) peaks at low values and then is smeared out over the rest of the range. P/Halley 
silicate dust more closely resembles the latter. When comparing the ternary system Mg-Si-Fe, 
Halley silicate dust again is closer to the anhydrous IDPs, particularly in the pyroxene region. 
However, there is a clear difference between Halley silicates and IDPs: P/Halley silicates are 
relatively poor in Si and rich on Mg [there are no IDPs with Mg/(Mg + Si) > 2/3, i.e., richer 
in Mg than olivines]. 


Humberto Campins 

Comment: The mass spectrometry of Halley’s dust confirms the groundbased 10 fim spec- 
trophotometry (Campins and Ryan, Ap. J. 341, 1059, 1989), which also shows that the 
silicates in Halley were anhydrous; more specifically crystalline olivine. 

Walter Huebner 

Yes, Halley silicate dust more closely resembles anhydrous interplanetary dust particles, par- 
ticularly in the pyroxene region, but most Halley silicates are richer on Mg and poorer on Si. 
Does this mean that they have to be crystalline? 


David Lynch 

Are there models or observations of grain fragmentation which could guide those of us who do 
thermal IR spectroscopy as to the distance range on which fragmentation is expected to occur? 

Walter Huebner 

Observationally there are indications that the ratio of small particles relative to large particles 
increases with cometocentric distance, bursts of small particles have been measured indicat- 


* V . * 

ing fragmentation at large cometocentric distances, and the images of “jets” near the nucleus 
show a behavior that is consistent with fragmentation at small cometocentric distances. To- 
morrow you will hear about progress of modeling fiagmentation, but particle lifetimes are 
only guesses at this time. 
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PHOTOMETRIC NARROWBAND CCD IMAGING OF 
^4 COMETS P/BRORSEN-METCALF AND AUSTIN (1989C1) 

] 

David G. Schleicher, David J. Osip, Robert L. Millis (Lowell Observatory), 
Andrea Thompson and Linda M. Sauter (Mass. Inst. Tech.) 


Abstract 

Wc have obtained simultaneous widc-ficld CCD images and conventional aperture 
photometry of Cornels P/Brorscn-Mctcalf (1989o) and Austin (1989cl). These data 
allow us to directly test our absolute calibration of the CCD images and our ability to 
generate full-coma, continuum-subtracted emission band images. Preliminary 
photometric calibration for a portion of the Brorsen-Mctcalf observations yields 
reduced emission band and continuum fluxes which agree to within about 10% 
with those from the photoelectric observations. Conventional photometry of 
Comet Austin shows unusual variations among the production rates of the 
different species as a function of heliocentric distance. 

Introduction 

Wc have recently embarked on a program to obtain full-coma, narrow-band CCD 
images of comets. In addition to obtaining model-independent total abundances, a 
primary reason for imaging the whole coma is to obtain a complete spatial map of 
each species. These maps, in turn, can be used to improve model parameters such 
as lifetimes and scalclcngths which arc required to properly extrapolate to whole- 

coma values from small aperture measurements made with typical photometers 

and spectrophotometers. 

Qh&giyaiiaits 

In order to obtain digital images of the entire coma of a comet and to adequately 
sample the sky background, wc have coupled a 200 mm Takahashi f/4 wide-angle 
telescope to our 800x800 TI CCD. This system gives a 50 aremin field with 4 

arcscc/pixcl resolution. A set of image-quality IHW filters (2x2-inch to eliminate 

vignetting) isolate the CN, C2, C3, CO + , and H20 + bands, plus three continuum 
points. In the last year, a large number of photometrically calibrated narrowband 

images of Comets P/Brorscn-Mctcalf and Austin were recorded with this 

equipment. A contour plot of a typical image is shown in Figure 1; note that the CN 
coma nearly fills even this large field. 

By fitting model profiles to the actual radial profile extracted from the image, sets 
of viable parent and daughter scalclcngths can be obtained. This family of 
solutions is considerably more restricted than would be the case if only the inner 

portion of the coma were measured. An example is given in Figure 2, where wc 

show both the radial profile of CN (from Figure 1) and a contour plot showing yp- 
for each combination of parent and daughter Hascr scalclcngths. The two 
innermost contours arc the domain of viable fits to the radial profile. The 
intersection of this domain with those found using images obtained at other 
heliocentric distances will progressively narrow the solution. 
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Figure 1. Contour plot of Comet P/Brorsen-Metealf as observed in the CN molecular 
band (3870A). The scale of the image is SO areminutes on a side. 




log Parent Scale Length (km) 


Figure 2. Radial profile of CN (left) extracted from the image shown in Figure 1 
after continuum subtraction, and a contoui plot of % ^ for model fits to the radial 
profile as a function of assumed parent and daughter scalclengths (right). The 
two innermost contours are the domain of viable fits to the radial profile. 
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Figure 3. Contour plot of Comet Austin (1989cl) as observed in the OH molecular 
band (3080A). Image scale is approximately 5 aremin on a side. 


The wide-ficld telescope rides piggy-back on the John S. Hall 42" telescope at 
Lowell Observatory’s Anderson Mesa site. Consequently, it has been possible to 
obtain conventional filter photometry at the same time as the narrow-band, wide- 
ficld images were being exposed. In this way, we have been able to independently 
test the photometric calibration of the images. These images will subsequently be 
used to obtain continuum-subtracted, emission band images reduced to column 
density per unit area. Preliminary reduction of a portion of the Brorscn-Mctcalf 
imaging data has given results within about 10% of those from the aperture 
photometry. We suspect that the remaining discrepancy may be due to slight 
differences in the filter characteristics — a possibility that will be explored 
further in the coming months. We emphasize that good absolute flux calibration 
of emission band images is attainable only when accurate sky, standard-star, and 
extinction measurements arc made under photometric conditions, thus permitting 
accurate continuum subtraction to be performed. 

Since glass elements in the optical system of the Takahashi telescope prevent 
working in the UV, we have also obtained CCD images of Comet Austin at the 
Cassegrain focus of the 42" telescope. These images give a measurement of the 
spatial distribution of the OH molecule (3080A) in the inner coma. A sample is 
shown in Figure 3. These data will also be used to constrain OH model parameters. 


*■ > 
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Figure 4t Production rates of Comet Austin as a function of heliocentric distance. 
Pre-perihelion points are shown as open circles, while post-perihelion points are 
shown as filled circles. Note the asymmetries about perihelion, as well as the 
differing behavior among the species. 


The conventional aperture photometry obtained by ourselves and by Peter Birch 

at Perth Observatory give excellent heliocentric distance coverage of Comet 
Austin. These data clearly show that the initial rapid increase in brightness from 
December to January (which led to predictions of the "comet of the decade") slowed 
considerably during February and March (see Figure 4). While the production of 
dust by Comet Austin was essentially constant prior to perihelion (possibly 

tracking the OH production), our observations since late-April show the dust 
production dropping at a rate similar to those of all of the measured gas species 
except OH, which held constant for several weeks before it began to drop. The 

apparent minimum in the dust production near the end of the interval covered by 

our observations may be an aperture effect because minimum geocentric distance 
occurred at this time. 

Comet Jialley . r J91Q 

In an unrelated development, we report the discovery by Schleicher and S. J. Bus 
of a 7.4 day periodicity in the brightness of Comet Halley during the comet's 1910 
apparition. Although this period was previously identified in photometry obtained 
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during the 1986 apparition (Millis and Schleicher 1986), considerable debate has 
cu?iicd as to whether this period directly reflects Halley's rotation period, or 

whether the nucleus is moving in a more complex manner. In fact, the initial 

evidence for a much snorter period of near 2.2 days emerged from the analysis of 

the morphology of dust jets recorded photographically during 1910 (Sekanina and 

Larson 1984; 1986). Reduction and analysis by Schleicher and Bus of brightness 
measurements of a homogeneous set of photographs taken oh a nightly basis at the 
National Observatory of Argentina at Cordoba (Perrine et al. 1934) has yielded a 

period of approximately 7.4 days in each of the two monthly intervals covered by 

the 1910 data. Further, a simple phase shift of the 1910 phased lightcurve with 

respect to 1986 data yield a believable match in the shapes of the lightcurves. 
Although not conclusive due to the uncertainties in the 1910 data, this result 
strongly suggests that not only was a period of 7.4 days present in 1910, but that 
the same source regions were active. A paper discussing these results has been 

submitted to the Astronomical Journal. 
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ABSTRACT 


Near-nucleus imaging observations of comet Austin 1989cl were 
carried out by the Japanese CCD imaging team. Six telescopes were used 
to monitor the time variation of the near -nucleus images in C 2 , CN, 
H 2 O*. Na. continuum in the optical region, and in J, H, K bands in the 
near- I R region. A featureless, round shape of the comet was revealed 
in all images. Although some of the jet features are recognized by 
using image enhancement technique, the azimuthal difference of the 
intensity distribution is about 10%. The images in the H20 + band show 
complex ion structures near the nucleus. 



1. INTRODUCTION 

** . ; 

' * 

• * '* ** 

%. * 

Comet Austin 1990cl was one of the important object for cometary 
researches after the apparition of comet P/Halley, because of its 
brightness along with its good observational condition. In order to 
study the near-nueieus phenomena, we observed this comet with CCD 
imaging systems attached to six telescopes in Japan. This poster paper 
shows some of preliminary results in these observations. 


2„ OBSERVATIONS 

The observations were carried out after the perihelion passage of 
the comet. Telescopes and CCD systems are listed in Table i. Two 
observatories, C.R.L. and S.K.A.O., were engaged for imaging 
observations of this comet from April 20 through May 31! However, we 
could observe the comet only 6 nights in this period because of bad 
weather condition. 

Near-nuceus images in the optical wavelength were taken from 
April 26 through May 26. We used the IHW standard filters of CN, C 2 , 
H 2 0* and correponding continuum bands. The R and I bands used at the 
K.O. and O.A.O. observatories are those in the Kron photometric 
system. The B and V at the K.O. are those in the Johnson system. The 
interference filters for Na D-lines used at the C.R.L. and T.O. are 
those of X c =589nm with A X j/ 2 = 4nm. Near-nucleus images in the 
near-IR wavelength were obtained at the C.R.L. observatory on April 
29, 30. May 9, 26, 27. The photometric bands are J, H, and K. 
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3. PRELIMINARY RESULTS 


The reduction is proceeding with using an image processing 
software library SPIRAL! Ichikawa et ai. 1987). Both the bias 
subtraction and the flat-field correction were applied with some 
relevant data, which were also taken during the observation. 

Figure 1 shows a C 2 image on April 30 taken at the T.O. This 
image shows a featureless, round shape. However, a broad, northward 
jet-like structure is recognized in figure 2, which is the result of 
contrast enhanced image in figure I, using the ring-masking 
method (Cosmov i ci et al. 1988). The azimuthal difference of the 
intensity distribution is about 10% at the distance of 2X 10 4 km from 
the nucleus. On the next night, the position angle of this jet changed 
counterclockwise by 30* -40*. This gives information on the nucleus 
rotation(Sekani na et al. 1981). The images of CN are so poor that we 
can obtain the CN-jet structureXA'Hearn et al. 1986). 

We are grateful to Dr. I. Konno for his arrangement of this paper 
presentation in the workshop. 
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A composition comparison between comets P/Halley and P/Brorsen-Metcalf 
Michael A. DiSanti* and Uwe Fink 
Lunar and Planetary Laboratory, University of Arizona 
Tucson, Arizona 8S72I 


The appearance of comet P/Brorsen-Metcalf in 1989 following only a few yea* a 
on the heels of comet P/Halley in 1986 gave us an unusual opportunity to compare 
two comets of similar orbital elements but potentially different physical and 
chemical characteristics. The period, perihelion distance and eccentricity for 
comet P/Halley are 76.0 years, 0.587 AU and 0.967 while for comet P/Brorsen- 
Metcalf they are 70.6 years, 0.479 AU and 0.972. 

We were able to acquire .photometric spectral data for P/Brorsen-Metcalf on 
two dates: 1989 July 13 and 1989 August 02. The observations were carried out 
with the 154 cm Catalina Telescope of the University of Arizona Observatories. 

The instrumentation consisted of our wide field CCD camera with reducing optics 
and a fast spectrograph covering 5200 - 10400 A, at a dispersion of 7.21 A/pixel 
and a slit length of 200" at a spatial scale of 1 .49"/pixel. Both the camera and 
spectrograph used an 800 x 800 Texas Instrument chip. Because of the reduced 
brightness of comet P/Brorsen-Metcalf we opened the slit to 5.0 " compared to the 
2,5" used throughout the comet P/Halley observations. This increased the ac- 
quired signal considerably while affecting the spectral resolution in only a 
minor way. 

The spectrum for 1989 July 13 is compared to a P/Halley spectrum in Fig: 1. 
The general appearance of the spectra is very similar, both comets showing 
emissions due to C 2 , NH 2 , OI (*D) and CN. Neither displays any emission features 
that are exclusive to only one of the two. The most notable difference between 
the two comets is the much lower continuum level of P/Brorsen-Metcalf as compared 
to P/Halley . Besides the differences in continuum, the spectra also show 
considerably weaker NH 2 emission for comet P/Brorsen-Metcalf as well as weaker CN 
emission and somewhat weaker C 2 emission when compared to the OI *D 6300 A 
emission line. 

Since water is now reasonably well established as the major driving force in 
the activity of a comet, we calculate the H 2 0 production rate for comet 
P/Brorsen-Metcalf in Table I . We followed the same procedure recently applied to 
comet P/Halley (Fink and DiSanti 1990) and therefore do not repeat a discussion 
of our data reduction. The resulting water production rate for comet P/Brorsen- 
Metcalf is 16% of P/Halley at r-1.36 AU and 20% of P/Halley at r=1.01 AU (Table 
I). . Both measurements refer to pre-perihelion dates. The P/Halley water 
production rate in Table 2 is obtained from the relationship Q(H 2 0)=4.36 x 10 29 
r -z.M derived in Fink and DiSanti (1990). For completeness we list the 
continuum flux, in Table 1, for the same aperture used for the water production 
rate. However, to effect a comparison in continuum fluxes, we turn to the values 
in Table 2 which used common apertures for both the P/Halley and P/Brorsen- 
Metcalf reductions. While for P/Halley the ratio of the continuum to OI flux is 
-0.18A' 1 for P/Brorsen-Metcalf it is -0.0085A" 1 . Thus P/Brorsen-Metcalf exhibits 
only about 1/20 of the continuum of P/Halley when both are compared to their 
respective H 2 0 production rates. 


■We next present the relative ratios between the production rates of the 
various species which can be obtained from our data. The dilution factor 
correcting the measured intensity to the total flux emitted into 4 steradians 
is the same for all species of a given comet. Secondly for comets observed at 
the same heliocentric distance the scale lengths are identical and the correction 
factors from the finite observing apertures to an infinite observing radius are 
the same, provided the radii of the observing aperture (i.e, slit lengths for us) 
are chosen to be the same. To meet this criterion, we picked common apertures of 
31 Q00 km for (he observation pair P/Brorsen-Metcalf 1989 July 13 and P/Halley 
1985 Dec. 08/09 and 43 000 km for the pair 1989 Aug. 02 and 1986 Jan. 10 (Table 
2). For the first pair we were fortunate that both objects were at the same 
heliocentric distance r*1.36 AU. For the second P/Brorsen-Metcalf observation we 
interpolated our P/Halley data for 1985 Dec. 08/09 (r»1.36 AU) and 1986 Jan. 10 
(r-0.87 AU) to the heliocentric distance of 1.01 AU. 

Our results are summarized itHTable 2. All species are compared to the H 2 0 
production rate of the two comets. We set P/Halley « 100% for each flux ratio 
which we argued above is equal to their production rate ratios. The largest 
difference between the two comets occurs for their relative NH 2 production. For 
both observation dates P/Brorsen-Metcalf exhibits about one third the production 
rate ratio of comet P/Halley. The parent of NH 2 is thus definitely less abundant 
in P/Brorsen-Metcalf. While comet P/Brorsen-Metcalf is the first comet for which 
we have carried out such a comparison quantitatively, visual inspection of our 
other comet spectra (e.g. Johnson, Fink, and Larson 1984) shows similar 
variability for the NH 2 emissions. C 2 (at 80%) seems to be only slightly lower 
at r-1.36 AU but considerably lower at r»l.01 AU (40%). Cyanogen is down a 
factor of two but participates slightly in the above trend being 65% at r*1.36 AU 
and only 48% at r«1.0l AU. 

The results presented in Table 2 give indications that C 2 originates from an 
extended source of dust particles or small polymers in addition to parent 
molecules. The lower production rates of C 2 and CN in P/Brorsen-Metcalf 
correlate well with both the lower dust production and their smaller increase as 
the comet approaches the sun. While the trend is there, we nevertheless note 
that the comparative dust production for P/Halley is down by a factor of twenty, 
but the production of C 2 and CN is only slightly off. 

Even without a model, our spectroscopic observations of P/Brorsen-Metcalf 
allowed us to draw significant compositional conclusions./ The NH 2 abundance is 
lower by about a factor of three which appears to be independent of heliocentric 
distance. Both C 2 and CN are less abundant in P/Brorsen-Metcalf than in P/Halley 
and do not seem to participate as much in the strong increase as the comet 
approaches the sun. The lower H 2 0 production rate implies that the active 
surface area of comet P/Brorsen-Metcalf is much smaller, or equivalently the 
comet as a whole is smaller. The substantially lower continuum level for 
P/Brorsen-Metcalf cannot readily be explained by smaller or larger active surface 
area, and must be intrinsic to the comet's embedded dust-to-ice ratio, or to an 
overall smaller-sized dust grain population in P/Brorsen-Metcalf. 
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Table I. M 2 () production rale and continuum flux for P/Halley and P/Brorsen- Metcalf 



Halley 

Brorsen- 

Metcalf 

Halley 

Brorsen- 

Metcalf 

Date 

85 Dec 08/09 

89 Jul 13 

85 Dec/86 Jan 

89 Aug 02 

Heliocentric distance (AU) 

1 36 

1.36 

1 01 

1.01 

Aperture size (arc-sec) 

151 x 2 5 

151 x 50 

151 x 2.5 

151 x 50 

(km) 

6300 A 

77800 x 1300 

97200 x 3200 

140000 x 2400 

71000 x 2400 

Ol Flux (photons/s m l ) 
6250 A 

1760 1 

490 

4000 * 

2600 

Cont. Flux (photons/s m* A) 

270 1 

6.4 

820 * 

32 

Luminosity Ol* (s'*) 

12.5 x I0 ,T 

2.0 x 10” 

26.6 x 10” 

5.4 x 10” 

H a O Prod Rate *•« (s l ) 

20.0 x 10” 

3.17 x 10” 

42.5 x 10” 

8.60 x 10” 

Q(H 2 0) Halley- 1 00% 

100% 

16% 

100% 

20% 


Notes: I. Average from Halley observations Dec 08 and 09 (Fink and DiSanti 1990). 

2. Interpolated between our Halley 85 Dec 08/09 observations (r - 1.36) 

and our 86 Jan 10/11 observations (r - 0.86). 


3. 6300 A Ol luminosity and H 2 0 production rate are independent of 
aperture used since they have been determined via a model. 

4. For the production rate of Halley we used the relationship Q(H 2 0) - 4.36 x 10” 

r -J sc (pi n |< an( j DiSanti 1990). 


Table 2. Comparison of emission fluxes between P/Halley and P/Brorsen-Metcalf 



Halley 1 

Brorsen- 

Metcalf 

Halley * 

Brorsen- 

Metcalf 

Date 

85 Dec 08/09 

89 Jul 13 

85 Dec/86 Jan 

89 Aug 02 

Heliocentric distance (AU) 

1 36 

1.36 

1.01 

1.01 

Aperture size (km) 
6300 A 

31000 x 1300 

31000 x 3200 43000 x 2400 

43000 x 2400 

Ol Flux (photons/s m l ) 
6250 A 

1450 

380 

3800 

2430 

Cont. Flux (photons/sm* A) 

230 

3.2 

770 

20 

F(continuum/OI) (A' 1 ) 

0.16 

0.0085 

0.20 

0.0085 


100% 

5.3% 

100% 

4.2% 

F(C,/OI) v— 1 

100% 

80% 

•00% 

42% 

F(NH 2 /OI) 0.10,0 

100% 

34% 

100% 

29% 

0.8,0 

100% 

32% 

100% 

32% 

0,6,0 

100% 

40% 

100% 

36% 

Avg. 

100% 

36% 

100% 

32% 

F(CN/OI) 2-0 

100% 

55% 

100% 

47% 

1-0 

100% 

73% 

100% 

50% 

Avg. 

100% 

65% 

100% 

48% 
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McDonald Observatory, The University of Texas, Austin, TX 78712 


Introduction 


A common feature in the visual spectra of comets is the A 2 A — X 2 II transition of CH at 4314A (At; = 0). 
This feature is very weak and is easily "lost” in the continuum or noise. However, the oscillator strength is 
very small, so even a weak feature indicates a large amount of CH. 

Understanding the parent of the CH in the comae of comets is very important for understanding the 
chemistry of comets and, thus, the early solar nebula. Potential parents for the CH are CH* or more 
complex hydrocarbons, including polyoxymethalenes (POMs). The ratio of CH* to H 2 O is indicative of 
possible sites of formation of the comets. 

Past observations of comets have shown that the CH is present in the coma. However, because of the 
intrinsic weakness of the feature and the ease of photodissociating CH, the distribution of the CH in the 
coma has only been traced in the inner coma, where the fall-off is still where p is the impact parameter. 

This paper presents the first data sets tracing the distribution of CII to large distances from the nucleus 
in two comets. We also will present some preliminary models pertaining to the parent of the CH. 


Observations 

The observations presented here are of periodic comets Halley and Brorsen-Metcalf. Both of these comets 
are Jupiter family comets with periods of around 76 years for Halley and 70 years for Brorsen-Metcalf. The 
inclination of Halley is 162° so the orbit is retrograde. Brorsen-Metcalf has an inclination of 19° and a 
prograde orbit. 

The observations were obtained on the 2.7 m telescope at McDonald Observatory using the Large Cass 
Spectrograph (LCS), a long slit spectrograph with a CCD detector. For the Halley observations, the CCD 
was an RCA 320 X 512 element chip with a relatively high read-noise of 50 electrons. For the Brorsen-Metcalf 
observations, the detector was a TI 800 x 800 element CCD with a read-noise of 12 electrons. An observing 
log is included in Table 1. 

For all of the observations, the spectral resolution (2 pixels) was around 7A. For the Halley observations, 
the data cover the range from 3650A-6800A while the spectral range for the Brorsen-Metcalf spectra was 
either 3200A-5950A or 3000A-5750A. For most of the observations the slit was 2arcsecs wide. 

The data were reduced from raw counts to fluxes using the IRAF data reduction package. The emission 
from the bands was then integrated and column densities computed. The conversion from integrated CH 
band intensity to column density was done using formula 1 of Cochran and Barker (1985). No Swings effect 
was taken into account. 

Table 1 notes that it was photometric for all of the Halley observations but not for any of the Brorsen- 
Metcalf observations. The Brorsen-Metcalf observations point out the advantage of using a long slit spec- 
trograph for observing comets; even without it being photometric, we are able to obtain profiles of the gas 
within the coma which are quite accurate since all of the portions of the coma are observed simultaneously 
in one spectral image. Although the column densities are not absolute, the relative values at different 
positions and between species are very accurate. In fact, for Brorsen-Metcalf, we obtained more than one 
spectral image per night, moving the telescope to observe different parts of the coma and were able to match 
the spectral images because we always obtained significant positional overlaps between adjacent spectral 
images. 

Figure 1 shows the distribution of the CH in the coma for the four nights of observations of comet 
Brorsen-Metcalf. The lack of scatter in the data gives a good indication that the errors are quite low. In 
most cases the CH information extends to 30,000 km. Figure 2 shows two representative nights of the Halley 
data (the data on the left hand and right hand panels is identical and the difference between the panels 


Table It Observing Log 


Comet 

Date 

Ha 

A 

Photometric? 

Halley 

9 Dec 85 

1.36 

0.72 

y 


20 Mar 86 

0.99 

0.81 

y 


21 Mar 86 

1.01 

0.79 

y 


22 Mar 86 

1.02 

0.76 

y 


23 Mar 86 

1.04 

0.74 

y 


24 Mar 86 

1.05 

0.71 

y 


10 Apr 86 

1.32 

0.42 

y 


2 May 86 

1.63 

0.80 

y 


9 May 86 

1.75 

1.05 

y 


10 May 86 

1.77 

1.08 

y 

Brorsen-Metcalf 

16 Jul 89 

1.31 

0.85 

n 


1 Aug 89 

1.04 

0.65 

n 


9 Aug 89 

0.91 

0.62 

n 


11 Aug 89 

0.87 

0.63 

n 


will be explained in the next section). These two nights were chosen because they represent two extremes 
in heliocentric distance. The Halley data show much more scatter than the Brorsen-Metcalf data. There 
are two reasons for this: 1) Halley has a much higher dust/gas ratio than Brorsen-Metcalf so continuum 
subtraction introduces a much higher error; 2) The Brorsen-Metcalf data were obtained with the much lower 
noise TI chip. 

Models 

The lifetime of CH against photodissociation is very short, ^lOOsec at 1 AU (Allen et al y 1989), because 
the threshold for photodissociation of CH is at 3500A. Thus, wherever one observes the Cl! in the coma, 
the direct parent must be very near. We have assumed that the direct parent of CH must be CII2. The 
photodissociative lifetime of CII 2 is 8,33 x 10 3 sec at 1 AU (Allen et a/, 1989), If we further assume that 
CH* is the direct parent of CII2 (the possibly logical daughter of CII4, CII3, is not seen in the lab), the 
photodissociative lifetime of CH4 at I AU is 1.35 x 10 5 sec (Allen cl «/, 1987). We can derive scale lengths 
from these lifetimes if we assume some sort of velocity law. For this work we have assumed the velocity 
law of Delsemme (1982). Utilizing these derived scale lengths, we have fit Haser models to the data. These 
Baser models predict the distribution of the CHa but, as explained above, the CII 2 should be very near 
(~100 km) the CH, so the models also predict the distribution of the CH gas which we can then compare 
to the data. 

Figure 1 shows the comparison of the models and the data. The fits are excellent. Note that the only 
assumptions we have made are a velocity law and the reaction pathway. The lifetimes are from straight 
physics. Thus, for Brorsen-Metcalf, the reaction pathway is plausible (but not necessarily unique). As 
has been often noted, the Haser model lacks for physical insight, but can empirically describe data well. 
This is certainly the case for CH. However, we also ran models of the physically more realistic Vectorial 
model (Festou, 1981). These models are almost indistinguishable from the Ifaser model. The only significant 
deviations between the models occur for data points whore the optocenter is within the spatial bin, probably 
due to how the integration is handled by the two codes. This is a region which also sufTers errors due to 
guiding and seeing efFects. 

We then tried to fit Haser models to all of the Halley data sets (not just the Two shown) using the 
parameters discussed above. These fits were significantly and consistently worse than the Brorsen-Metcalf 
fits, with the greatest problem seen in the inner coma. I11 the left hand panels of figure 2, note that 
the models over-predict the CH in the inner coma. However, the Delsemme velocity law predicts that 
the outflow velocity should have been 0.61km/scc at the time of the spacecraft encounter. The Giotto 
spacecraft measured a velocity of 0.9km/sec (Krankowsky et al y 1986). If one assumes that the scaling of 
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log Column Density (mol cm 2 ) 


Brorsen-Metcalf 



0 2 4 0 2 4 

log Radial Position (km) 


Figure 1: The distribution of Oil gas for the A nights of observations of Brorsen-Metcalf. The Xs represent 
the data The small scatter band indicates the high quality of the data The line is the model described in 
the section on models 




a 


fV-1 01 au 
A-0 56 au 
Doisemma Velocity 


21 Mar 1086 
R*-1 01 au 
A-0 56 au 
Spacecraft Vetoaty 


10 May 1966 
R*-1 77 au 
A-1 06 au 
Delsemme Veloaty 


10 May 1966 
R^-1.77 au 


A-1 06 au 


Spacecraft Veloaty 


log Radial Position (km) 


Figure 2: The distribution of the Cll gas for 2 representative nights (smallest and largest heliocentric 
distance) of comet Halley. The Xs represent the data The data in the left hand and right hand panels are 
identical The line represents the models described in the text: The left hand panels utilizes the Delsemme 
velocity law while the right hand panels utilizes a law tied to the spacecraft velocity. 






the Delsemme law is correct, but that the zero point is incorrect for Halley, one can derive a new scaling 
law of 

v = 0.85fl*°’ s 

for the Halley data. If these new velocities are used along with the lifetimes described above, when can 
derive new scale lengths for use with the Hascr model. The fits with the spacecraft derived velocity scaling 
law are shown in the right hand panels of figure 2. Although these fits are not as good as the Brorsen-Metcalf 
fits, they are a significant improvement over the fits with a Delsemme law. This is true of all of the nights* 
data, not just the two shown in figure 2. It should be noted that there are actually several derived outflow 
velocities for Halley and these are summarized in Lammerzahl et al (1987). 

Thus, there is a Haser model which fits each of the comets using the same lifetimes, but assuming 
different velocity laws. This is not the first time that we have noted that Halley has a different outflow 
velocity than other comets. This appears to be true for fits of NHj also (Cochran et al, 1989). However, 
Halley is the only comet we have studied in detail which does not fit the Delsemme velocity law. It should 
be noted here that the Brorsen-Metcalf data are totally incompatible with the higher Halley spacecraft 
velocities. 

Conclusions 

We have presented some of the first profiles of CH in the coma of a comet which are traced to large radial 
positions. These gas distribution profiles can be fit well using a Haser or Vectorial model and the following 
reaction pathway: 

CH\ -f hu — CH 2 + H 2 
CH 2 + hu -* CH+H 
CH + hu -» C + H 

The lifetimes for these reactions are the same for both comets and are given in the previous section. However, 
the velocities used to convert the lifetimes to scale lengths are different for the two comets: Brorsen-Metcalf 
velocities are derived from the Delsemme velocity law while Halley velocities are derived from a scaled 
spacecraft velocity. 

Although the Haser model has no physical basis, it does a good job of describing the data with the scale 
lengths which we derived. This shows that the reaction pathway which we have employed is a plausible, 
but not necessarily unique, reaction pathway. We cannot rule out more complex hydrocarbons, although 
we do not need to invoke them either. More complex reaction pathways or parents are possible, as long as 
the scale lengths for these more complex reactions are comparable to the scale lengths for the CH« pathway 
given above. 

This work was supported by NASA Grant NAGW 1477 and NSF Grant AST-8912869. 
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Patrick Palmer 

It was interesting that the Delsemme velocity law fits comet Brorsen-Metcalf, but the constant 
multiplier had to be changed by ~ 50% for Comet Halley to agree with spacecraft instruments 
and to make your data fit for P/Halley. Could you give an estimate for tht error ‘in'mulitplier 
of this velocity law, as found by your measurements of Comet Brorsen-Metcalf? 

Anita Cochran 

In the past, we have found the Delsemme velocity law to hold for most comets and have found 
other species which point to the higher spacecraft velocity for Halley. The spacecraft velocities 
are 1.47 times higher than the Delsemme law. Our Brorsen-Metcalf data set is so clean one 
can easily rule out Halley-like velocities for Brorsen-Metcalf. I haven’t tried putting a num- 
ber to it, but I would imagine that one can say that Brorsen-Metcalf follows the Delsemme 
velocity law to about 10%. Of course, that assumes we have the lifetimes correct! 


Susan Wyckoff 

Ferro has also modeled the CH« photolysis. A comparison of a very preliminary version of his 
model with our observed P/Halley CH profile using van Dishoeck’s photodissociation lifetimes 
does not indicate close agreement. Perhaps the difference in our results reflects the difference 
in the photodissociation lifetimes that we used in modeling the CH profiles. 

This work is in progress. 

Anita Cochran 

We did a quick comparison of Ferro’s model to our data during the coffee break and his model 
deviates from our data about the same way as it deviates from your data. An eyeball com- 
parison of the data sets show them to compare favorably (except where yours turn back up 
• something which is non-physical). Clearly, we will have to check your scale lengths against 
our model more formally and you will have to check your model carefully. 


Lewis Snyder 

How do we reconcile the CO, HjCO and CN extent of 1.5 x 10 4 km found for Comet Halley 
with the opacity radius on Comet Brorsen-Metcalf? 

Anita Cochran 

We do not see any conflict. First of all, Brorsen-Metcalf has an extremely high gas-to-dust 
ratio. Secondly, the spacecraft were able to image the nucleus of Comet Halley, even at 1300 
km. The phenomena you talk about with molecules such as CO is’ not optical depth but is 
an extended source. 


AMMONIA AND NITROGEN ABUNDANCES IN COMETS 
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Department of Physics and Astronomy 
Arizona State University 

Introduction 

Comets consist of pristine material preserved from an earlier galactic epoch. Deter- 
mination of the molecular, elemental and isotopic abundances in the dust and volatile 
components of comet nuclei provide vital clues to the chemical evolution of both inter- 
stellar and solar nebula matter. Here the abundances of nitrogen-containing molecules in 
comets are considered. The molecular abundances of NH3 in four comets are summarized 
(Tegler 1990, Wyckoff, Tegler and Engel 1990). From an inventory of nitrogen-containing 
compounds (Wyckoff, Engel and Tegler 1990, Wyckoff, Engel, Womack, Ferro, Tegler and 
Peterson 1990), an estimate of the elemental N abundance is also presented, 

Ammonia Abundance 

The hypothesis that NH3 photolysis can account entirely for the NH2 and NH observed 
in comets Halley and Thiele is tested by calculating the spatial distributions of NH and 
NH2 assuming an NH3 photolysis model, 

NH 3 -* NH 2 — ♦ NH -+ N 

\ / 

NH 

In this photodissociation sequence, NH3 is the only source of both NH and NH2. The 
branching ratio into the NH2 sequence is 97%, while direct photodissociation into NH 
occurs 3% of the time. A multi-generational Monte Carlo model (Ferro 1990) used to 
compute the NH distribution for comet Halley is compared with observations in Figure 1. 
In Figure 2 we present the NH2 profile observed in comet Halley together with the profile 
computed using the vectorial model (Festou 1981), As input parameters to the models 
we have used photodissociation lifetimes calculated using the Solar Maximum Explorer 
satellite EUV flux for the dates of observation. In the vectorial model we have assumed 
an NH 3 outflow velocity of 0.8 km/s and that for NH2 of 1.2 km/s. In the Monte Carlo 
model the velocities were calculated explicitly assuming conservation of momentum and 
energy. We assumed that 10% of the residual photon energy upon photodissociation of 
a molecule is transformed into internal energy, and 90% is converted into kinetic energy 
of the photodissociation products, following Tegler(1990). The fits of the models shown 


in Figures 1 and 2 correspond to equal production rates of NHj, namely Q(NH 3 )= 3.0 x 
10 20 s -1 . We not only find that the shapes of the NH and NH2 spatial profiles are well- 
predicted with the NH 3 photolysis model, but we also find that the absolute scaling of the 
column densities of both NH and NH2 is consistent with a common source, namely, NH3. 
Therefore the NH3 photolysis model is self-consistent and provides strong evidence that 
both NH and NH 3 are the predominant photodissociation products of NH 3 . 

The agreement in the productions rates, Q(NH 3 )~Q(NH2) ~Q(NH) not only indicates 
that the NH 3 photolysis model adequately explains the NH and NH2 observations of comet 
Halley, but also puts the ammonia abundances derived from the optical spectra on a very 
firm basis. The quality of the model fits to the data also indicates that no other species 
(e. g. hydrazine, N2H4) contribute significantly to the observed profiles, as Tegler (1990) 
has previously also demonstrated. 

An analysis of four comets (Thiele, Halley, Hartley-Good and Borrelly) yields an av- 
erage ammonia abundance, NH3/H2O =0.0013 (Wyckoff, Engel and Tegler 1990). In ad- 
dition Tegler (1990) has found that the post-perihelion ammonia abundance agrees within 
a factor ~2 with the pre-perihelion value. From this we conclude that the ammonia abun- 
dances in this small sample of dynamically heterogeneous comets is remarkably uniform, 
and that comet Halley is chemically homogeneous, at least on a scale of ~6 meters which 
was the pre- to post- perihelion erosion depth (Keller et al. 1987, Tegler 1990). 

Inventory of Elemental Nitrogen 

An inventory of nitrogen-containing compounds in comet Halley has been compiled 
(Wyckoff, Tegler and Engel 1990). We find that NH3 and CN are the predominant gas 
carriers of N, and that N2 is a very minor source of the elemental nitrogen in the gas coma 
of this comet. Using the analysis of the comet Halley dust by Jessberger and Kissel (1990), 
we find that ~95% of the nitrogen is partitioned into the dust grains in comet Halley and 
only ~5% of the nitrogen is observed in the gas. Some nitrogen could be hidden from 
detection in large molecules such as cyanide polymers (Matthews and Ludicky 1986). In 
fact a weak feature detected in the spectrum of comet Halley at 4.45 fim has been attributed 
to a -C=N stretch transition (Encrenaz and Knacke 1990). Estimates have been made that 
<10% of coma material is hidden in undetected grains or large molecules (e. g. Ip et al. 
1990, Encrenaz and Knacke 1990). Thus even if all of the hidden(unobserved) material 
were in such nitrogen compounds, estimates of the amount of material unaccounted for 
in our inventory indicate that the large depletion of elemental nitrogen found for comet 
Halley would be virtually unchanged. 

The elemental abundances in comet Halley, including the nitrogen abundance deter- 
mined from our analysis (Wyckoff, Tegler and Engel 1990) are shown in Figure 3. As the 
figure indicates the nitrogen abundance in the gas component of the coma is ~l/40 N©, 
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and the elemental nitrogen abundance in the combined gas -f dust components of the coma 
of comet Halley is depleted relative to the sun, N(gas -f dust) ~ 1/10 N©. Here we have 
assumed a dust/gas mass ratio of two. The elemental abundances in the dust component 
were taken from Jessberger and Kissel (1990), and the abundances of light elements other 
than N were taken from Encrenaz and Knacke (1990) and Wegmann et al. (1987). The 
abundances in Figure 3 have error bars corresponding to factors of two. Analyses of the 
Giotto data have previously indicated deficiencies, although smaller, in the elemental ni- 
trogen abundance in comet Halley (Geiss 1987, Encrenaz, d’Hendecourt and Puget 1988, 
Ip et al. 1990). These previous analyses adopted NH3/H2O ~0.01-0.02, which is a factor 
~5-10 larger than used here, which accounts for the difference between this and previous 
N inventories in comet Halley. An estimate of the N2 abundance in comet Halley confirms 
that it is a minor species with, N2~l/10 NH3 (Wyckoff, Tegler and Engel 1990). 

A likely explanation for the nitrogen depletion in the comet is that a large amount 
of N2 was excluded from the condensation process in the solar nebula. While N2 can 
be incorporated into a clathrate hydrate, the process is very inefficient compared to en- 
clathration of CO and CH 4 . Thus N2 would not have competed effectively with CO and 
CH 4 for the H2O clathrate cage sites in the comet forming environment. Therefore unless 
temperatures in the comet-forming region were <20 K, N2 is expected to have condensed 
into comet nuclei very inefficiently. In this scenario, an N2-rich gaseous residue would have 
accumulated in the solar nebula as condensation and coagulation of cometesimals pro- 
ceeded. Photodissociation, photoionization of the N2 gas and subsequent expulsion from 
the solar system by radiation pressure, solar wind or other processes would have cleared 
the excess gas from the proto-solar accretion disk, leaving the comet nucleus depleted in 
elemental N. If on the other hand, the comet formed in temperatures <20 K where N2 
gas could have condensed as pure N2 ice, then subsequent diffusion due to heating of the 
comet nucleus would have selectively depleted the highly volatile gas, N2* In this case, the 
more thermally evolved a comet, the lower its N2 abundance would be. Of course, some 
combination of condensation fractionation and thermal diffusion could explain the N defi- 
ciency. The fact that CO is less depleted than N2 can be explained by the relative facility 
with which CO can bond with other compounds, especially to form organic molecules. In 
this regard we note that most(~2/3) of the CO observed in comet Halley originated in a 
distributed source, and not the nucleus. Thus the abundance of CO would not be expected 
to correlate very tightly with that of N2. The above scenario assumes that the comet ices 
condensed from the solar nebula, and that the bulk of the volatiles were not preserved as 
ice from their presolar environment. 
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Summary 

We have demonstrated that the observations of NH and NH 2 in comets can be un- 
derstood by assuming production of both species solely by NH 3 photolysis. The average 
ammonia abundance found for four comets was NH 3 /H 2 O ~0.0013 with an estimated un- 
certainty of a factor ~2. An inventory of the nitrogen compounds in comet Halley indicates 
that CN and NH 3 contain virtually all of the N in the gas coma, and that the dust coma 
contains ~95% of the total elemental nitrogen in the comet. We find that the elemental 
nitrogen is depleted by a factor ~40 in the gas coma, and a factor ~10 in the combined 
gas and dust coma. The depletion of elemental nitrogen in the comet may be the result 
of condensation fractionation of N 2 in the solar nebula or preferential diffusion of N 2 from 
the comet nucleus, or both. In a related program, we are surveying star-forming regions in 
the galaxy to determine the N 2 abundances in gas being assimilated into protoplanetary 
accretion disks (Womack 1990, Womack, Ziurys and Wyckoff 1990). 
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Walter Huebner 

How did you determine the N? abundance? 



Comment: Paul Feldman 

In response to W. Huebner’s question: on N 2 : Loss of N 2 in Earth’s atmosphere is principally 
through ionization (not dissociation). Rates for N 2 and CO are similar, so that to first order, 
[N 2 ]/[CO] is probably not too different from [N 2 ]/[CO + ] measured at a given point in the tail. 

Susan Wyckoff 

The estimate of the N 2 /NH 3 we have given for Comet Halley (Wyckoff and Theobald 1989, 
Adv. Space Sci 9) was based on the above suggestion by Paul Feldman to measure the 
N 2 /CO + ratio in our spectra. The ratio of the ion fluxes is relatively model independent, 
and with the ionization rates for N 2 and CO give an estimate of the N 2 /CO ratio. We used 
C0/H 2 0 ~ 0.1 to convert to an abundance, N 2 /H 2 O ~ 10“ 4 . The result is acknowledged to 
be an estimate only with an accuracy of a factor of several. Even allowing for such a large 
error, the elemental nitrogen abundance would still be deficient by an order of magnitude in 
Comet Halley. 


Paul Weissman 

What is the condensation temperature for N 2 in the solar nebula? 

Susan Wyckoff 

The condensation temperature for pure N 2 ice is very low 20 K. (See for example Prinn 
and Fegley, in “Origin and Evolution of Planetary and Satellite Atmospheres,” U. Arizona 
press, p. 78, 1989). 
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/ „• UNIDENTIFIED IONS IN COMETS 

. Lisa Engel 
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Optical spectra were taken of comets Halley, West and Brorsen-Metcalf in the wave- 
length region 3200-6400 A. These spectra were offset ~ 10 -5 — 10“ 6 km tailward from 
the nucleus so that the features detected were all ions with the exception of a very small 
residual Co emission at 5165 A. The full labelled tail spectrum of comet Halley is shown 
in Figure 1. While most of the features detected are attributable to CO + , H20 + , COf, 
CH + and OH + , there are three moderately strong bands in the spectra of comets Hal- 
ley and Brorsen-Metcalf which remain unidentified. These features were not detected in 
comet West. All three spectra in the 4800-5400 A region are shown in Figure 2, in order of 
decreasing gas-to-dust ratio, Brorsen-Metcalf being the gassiest of the three. The central 
wavelengths of the unidentified features are 4930 A, 5300 A and 6000 A. In an effort to 
identify the ions responsible for these features as well as to confirm previous identifications, 
laboratory spectra of ions were compared with the comet spectra. 

HoO + and CO + are known to have extensive emission in the 4800-6200 A region. 
Therefore, the possibility that these unidentified features are due to H20 + and CO'“ u was 
investigated. A search for the following bands was conducted: H20 + (11-0), H20 + (12-0), 
H 2 0+ (13-1), CO + (1-1) and CO + (0-0). CO + (1-1) and CO + (0-0) have been previously 
identified and are present in all the spectra examined here. No positive identifications of 
the H,0 + (12-0) and (13-1) bands could be made. However, the H20 + (11-0) band was 
identified. Comparison of the H 2 0 + (11-0) synthetic spectrum with comet Brorsen-Metcalf 
data is shown in 'Figure-3-. The relative fluxes of the blended H 2 0 + lines in the synthetic 
spectrum match those of the cometary data. The central peak at 5198 A may be due to 
residual night sky [NI] contamination in this region of the data. 

After eliminating and CO + as significant contributors to the stronger uniden- 

tified features, spectra of ions which have not yet been identified in the optical region of 
comets were compared to the data. These ions are NH + , CS + and C f '(Figures 4, 5, 6). 
None of these ions appears to be a significant contributor to the optical spectra of these 
three comets. 

In conclusion, H^O^ (11-0) has been unambiguously identified in cometary optical 
spectra. However, no other molecular ions for which there are laboratory spectra can 
explain the unidentified features found in the ion tails of comets Halley and Brorsen- 
Metcalf. 
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Figure 3. Comet Brorsen-Metcalf ion tail spectrum compared with the synthetic 
spectrum of the HjO -4- (11*0) band. 
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Figure 4. Comparison of NH + synthetic spectrum with comet Halley ion tail spec- 
trum. Note the lack of wavelength coincidence between the NH+ (0-0) feature and the 
actual data. The Swings effect has been taken into account. 
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Walter Huebner ** * ■ * 

Did you compare the unidentified lines with the very detailed spectrum of Greenstein and 
collaborator (1950’s) which contained many unidentified lines? 

Lisa Engel 

That is the Comet Mrkos spectrum which covers a spectral region in the ultraviolet. We have 
concentrated here on the optical region. 

William Cochran 

You rejected several possible identifications of lines based on the presence of ions in the the- 
oretical spectra that were absent in the data. However, the Swings effect can easily reduce 
the intensity of a single line by a large factor. Did your theoretical calculations include the 
Swings effect? 


Lisa Engel 

Yes, we did include the Swings effect. There are, however, problems with several coincidences. 


Humberto Campins 

Comment: There’s a 1965 or 1966 spectrum by C. Slaughter of Comet Ikeya-Seki. 

Lisa Engel 

There are several unidentified features in this spectrum that are of interest to this work. 
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SIMULTANEOUS IMAGING OF 
OPTICAL CN LINES AND RADIO HCN LINES 
IN COMET AUSTIN 


Patrick Palmer (U. of Chicago) 

M. F. A ’Hearn (U. of Maryland) 

Imke de Pater (U. of California - Berkeley) 
James J. Klavetter (U. of Maryland) 
David Mehringer (U. of Chicago) 

F. Peter Schloerb (U. of Massachusetts) 
Lewis E. Snyder (U. of Illinois) 

D. Wilner (l ? . of California - Berkeley) 


Introduction 

The parent molecule of cometary CN has been the subject of speculation for a long 
time. When HCN was detected at 3.4 - mm in Comet Kohoutek (Huebner etal. 1974), the 
problem seemed to be solved. But the much more detailed work on Comet Halley 
(Bockelee- Morvan etal. 1986; Schloerb etal. 1987) raised the quantitative question of 
whether HCN could be the only parent. Therefore comparative CN/HCN studies are 
vital for understanding the origin of cometary CN. 

The striking observation of CN ’’jets* in Comet Halley (A ’Hearn etal. 1986) raised 
another very interesting question about the origin of CN. Traditional theory permits only 
dust features to remain well defined far from the nucleus. Therefore the CN jets were 
interpreted as arising from submicron sized dust particles, perhaps CHON particles. An 
estimated 10 - 50 % of the CN in the comet was in the jets in the Halley observations. 

One can make several hypotheses: 1) some of the CN originates from the dust and has 
nothing to do with HCN ; or, 2) at least some of the HCN is also produced from the dust 
in the coma rather than directly from the nucleus. (In the second hypothesis, whether 
HCN is or is not the parent of the CN associated with the dust would need to be 
established.) The extended scalelengths found for CO by Ebernardt etal. (1987) and for 
HiCO by Snyder etal. (1989) also support ideas like hypothesis 2)-. A third hypothesis 
should be mentioned: contrary to the usual theory, the gas flow' does not become isotropic 
(Combi 1987). 

For all of the reasons mentioned above, it is essential to make a detailed comparison of 
the spatial distributions of CN and HCN. Furthermore, because of the variability of 
cometary emissions, it is necessary to make measurements simultaneously. 

Instrumental Considerat ions 

If one wishes to compare radio data with data like that of A ’Hearn etal. , rather 
stringent conditions are placed on the instrumental parameters. The images shown by 
A’Hearn etal. are ~ 5' (= 10 5 km) square. To obtain even 100 radio pixels in this area 
would require better than 30 " resolution. We also would like to have a resolution that 
corresponds to a fraction of the photodissociation scalelength of HCN (~ 40 " in this case). 

The radio project was designed to produce images with ~ 15 " resolution over field 2'.3. 




(During the period of observation this corresponded to ~ 6500 km resolution over a field 
~ 60000 km in diameter.) In order to do this two instruments were used: the Five College 
Radio Observatory (FCRAO) 14 - m antenna and the Berkeley - Illinois - Maryland Array 
(BIMA) of three 6 - m antennas. The FCRAO antenna has 58 w resolution at this 
wavelength, the BIMA array has a *2!3 field of view; its resolution depends on the locations 
of the antennas. 

The instrumental limitation of the single dish antenna is that information is collected 
for only one beamwidth at a time. In order to map an area, the dish must be pointed at a 
sequence of different positions. In sum, one has only one resolution available at each 
frequency, and one must sacrifice signal - to - noise ratio in proportion to the size of area 
mapped. 

Arrays have a different set of limitations. To understand the limitations, one must 
recall first that one measures the Fourier transform of what one wishes to know. The 
variable conjugate to the angular displacement in the sky is the projected baseline length. 
Therefore, because there is a maximum baseline, there is a minimum possible resolution. 
There is also a minimum projected baseline - the dish diameter - (otherwise one of the 
dishes shadows the other), so there is also a maximum size of feature to which an array is 
sensitive. Within these limitations, the array can make an image of a field of view as large 
as the primary beam of the individual antennas. Second, one does not measure the entire 
2-dimensional Fourier transform of the sky brightness in any finite experiment, but only 
some of the Fourier components. The problem then is to reconstruct an estimate of the sky 
brightness from the finite set of data points. For an array like the VLA with 27 antennas, 
yielding 351 simultaneous independent measurements, in one day enough data is obtained 
to permit accurate reconstruction of the sky brightness over the range of sizes to which the 
array is sensitive. All existing mm - arrays have too few telescopes to permit this. With the 
BIMA array, one has 3 antennas, and therefore only 3 simultaneous independent 
measurements. The successful scientific programs carried out with mm - arrays so far 
typically use 5 configurations of antennas. The data acquisition is, then, spread out over a 
minimum of 5 days (for practical reasons, more typically a few months). For a comet, 
where we cannot trust that what we see will stay the same from day to day, one cannot use 
multiple configurations. Therefore, none of the existing mm- arrays obtain data fast enough 
to permit the construction of reliable images of comets. 

The compromise approach that we chose was to use both a single dish and an array. 
This has several advantages. First, the single dish lias no problem with sensitivity to large 
scale structures. Therefore, no part of the HCN distribution will fail to be observed 
because it is too smooth. Second, by judicious arrangement of the 3 antennas in the BIMA 
array, complementary data can be obtained which allow the resolution to be improved by a 
factor of 2 - 3 beyond that of the FCRAO antenna alone. Numerical simulations were 
carried out to verify the result of combining the single dish and the array data. 
(Unfortunately, the simulations were based on Comet Austin being at least as strong an 
HCN emitter as Cornet Halley.) 

Observations 

The radio observations were scheduled at FCRAO and BIMA for 1990 April 24 
through 1990 April 30. Because of bad weather in New England and other problems, valid 
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Table 1*. Comet Austin Observations 


Date 

r 

A 

scale(20 h ) 

^HCN 

T t 

Apr 27 

.598 AU 

.626 AU 

9080 km 

19000 km 

0.13 day 

Apr 28 

.619 

.606 

8790 

20400 

0.13 

Apr 29 

.640 

.587 

8510 

21800 

0.12 

Apr 30 

.660 

.568 

8240 

23200 

0.12 


data was obtained only April 27 through April 30. Some relevant parameters are 
summarized in Table 1. The first three columns are self - explanatory, the fourth is the 
linear size at the comet that corresponds to 20 ", the fifth is the photodissociation 
scalelength of HCN, and the final column is the minimum time for HCN to move one 
resolution element. (An expansion velocity of 0.8 kms -1 was assumed in the last two 
columns.). The velocity resolution was 0.25 kms -1 at BIMA and 0.34 kms -1 at FCRAO. 
The total velocity range covered was ± 15 kms -1 at BIMA and ± 42 kms -1 at FCRAO. 
(Velocities are referred to the comet.) 

Because Comet Austin was fainter than anticipated, it was decided not to attempt to 
map the entire field of view of the BIMA array with the FCRAO telescope. Therefore the 
chance of detecting the comet was improved, but the final field of view will be reduced to 
about 1 '. 

The optical observations were made at Lowell Observatory using the 42 inch J. S. Hall 
reflector equipped with the U. of Maryland CCD camera. CN images were obtained on 
April 20.5, 22.5, 23.5, 27.5 and May 9.5. Therefore, on only one night was data obtained 
simultaneously with both radio telescopes and the optical telescope. 

Data reduction is still in progress. We have examined images made from each day’s 
data at BIMA, and there are no miracles. Similarly, the raw spectra at FCRAO were 
examined; and, at best, there is a possible detection in this first pass through the data. 
Later this summer we expect to have the final calibrated and edited data from the two 
instruments so that we can combine them. 

The Future 


One may infer (correctly) that because of the faintness of Comet Austin, the most 
valuable feature of the current experiment may turn out to be that it provided a ’’dry run” 
for a future, better comet. It is also instructive to examine how the instrumental situation 
will improve for the radio part of the experiment. 

First, the ’’speed” of the single dish observations will be improved enormously with 
the use of the very recently completed array feed at FCRAO. This feed provides 15 
simultaneous independent beams for the 14 - m antenna. Therefore the entire field of view 
of the BIMA array can be imaged with only 5 pointings of the telescope (rather than 75 to 
do the same thing when we observed Comet Austin). 

Second, more antennas are being added to the BIMA array. Mechanical assembly has 
just been completed for 4 new antennas. Therefore, we will have a 6 element array (usable 
to 1 - mm wavelength) within 2 years. (Three additional antennas are expected to be in 
operation shortly afterward.) With 6 antennas, one has 15 baselines, i.e. as many as are 
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obtained in many multi - configuration studies now. Therefore we will be at the stage of 
obtaining useful images in a single day from the array alone. 

Third, the National Radio Astronomy Observatory has proposed building a mm - 
array of 40 8 - m antennas, to be operated in an innovative way so that the problem of a 
maximum size scale will be greatly relaxed. With 780 simultaneous baselines, one can 
obtain excellent images in a very short time. The total collecting area will be ~ 2000 m 2 , 
i.e. larger than that of any other mm - wavelength instrument in the world (and 12 times 
that of the 6 element BIMA array), so the resulting sensitivity will be correspondingly 
greater. Although many other astronomers might have something to say about the division 
of observing time, it is clear to me that this will be a wonderful comet telescope. 
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Anita Cochran 

Your title talks of simultaneous optical images. What was Mike A ’Hearn able to get from 
Lowell Observatory? 

Patrick Palmer 

Mike obtained CN images on April 20.5, 22.5, 23.5, 27.5, and May 9.5. Therefore one was 
simultaneous with the radio observations. 
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FABRY-PEROT OBSERVATIONS OF COMET AUSTIN 
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We present preliminary results of a program to observe Comet Austin (1990cl) from 
16 April to 4 May and from 11 May to 27 May 1990 using the West Auxiliary of the 
McMath Solar Telescope on Kitt Peak, Arizona. The observations were made with a 
15-cm dual-etalon Fabry-Perot scanning and imaging spectrometer with two modes 
of operation: a high resolution mode with a velocity resolution of 1.2 km/s and a 
medium resolution mode with a velocity resolution of 10 km/s. Scanning data was 
obtained with an RCA C31034A photomultiplier tube and imaging data was 
obtained with a Photometries LN 2 cooled CCD camera with a 516 by 516 Ford chip. 
The results include: (1) information on the coma outflow velocity from high resolu- 
tion spectral profiles of [O IJ6300 and NH 2 emissions, (2) gaseous water production 
rates from medium resolution observations of [O I]6300, (3) spectra of H20 + emis- 
sions in order to study the ionized component of the coma, (4) spatial distribution of 
H 2 0 + emission features from sequences of velocity resolved images (datacubes), and 
(5) spatial distribution of [O I]6300 and NH 2 emission from medium resolution 
images. The field of view on the sky was 10.5 areminutes in diameter. In the imag- 
ing mode the CCD was binned 4 by 4 resulting in 7". 6 per pixel and a subarray 
readout for a field of view of 10'.5. 

i 

High-resolution scans were obtained for heliocentric distances r h of 0.46, 0.58, and 
0.94 AU (see figure 1.) The [O IJ6300 lrnewidth w (FWHM) as a function, of 
heliocentric distance varied approximately as rj, . Also, the [O I]6300 and NH 2 


^Visiting Astronomer, National Solar Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under contract with the National Science Foundation. 
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6298.6 linewidths were significantly less than those of Comet Halley's at 
comparable values of r^. 


[O I] and NH Linewidths 



ln(w) = (0.7310.05) + (-0.5810.09) ln(i^) 

Figure 1: [O IJ6300 and NH 2 6298.6 Linewidths 

Post-perihelion water production rates were determined from medium-resolution 
[O IJ6300 scans at heliocentric distances ranging from 0.44 to 1.11 AU (see figure 2.) 
The precision of these rates will be increased after careful study of [O IJ6300 imag- 
ing data to determine reliable aperture corrections for the amount of emission out- 
side our field of view. The [O IJ6300 is cleanly resolved from any of the nearby NH 2 
lines. 

Medium-resolution scans of the 6159 were obtained throughout the observa- 

tional program. On four nights, both H 2 0 + 6159 and 6147 doublets were scanned. 
The intensity ratio of these two doublets (see table 1) can be used to test models of 
cometary ionospheres. In theory, our spectrometer could resolve the 6147 

doublet (6146. 80 A and 6147. 38 A), but the cometary lines were broad enough 

to blend the two components. 
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Preliminary Comet Austin 
Water Production Rates 



Figure 2: Preliminary H 2 O Production Rates Based 
on [O IJ6300 Intensity Measurements. 


rp, , 1. u 1(6 147) 
Table l.H 2 0 I(6159) 


Date 

1(6147) 

(arbitrary) 

1(6159) 

(arbitrary) 

6147 

6159 

21 May 

191114 

170120 

1.110.2 

23 May 

97113 

104116 

0.910.3 

25 May 

96112 

60130 

1.511.0 

26 May 

115114 

73118 

1.610.5 


A large amount of time in our program was used to obtain velocity-resolved images 
(datacubes) of cometary H20 + emission in the medium resolution mode. The band- 
pass was 0.2A and the images were taken 0.2A apart. Datacubes of the 6159 doublet 
were obtained on 28 April and 3, 19, and 20 May. Datacubes of the 6147 doublet 
were obtained on 21, 23, 25, 26, and 27 May. Each image in a datacube was 
obtained in a five to fifteen minute exposure, with a complete sequence taking one 
to two hours to complete. A potential problem in analyzing these images will be the 
possible time variability of the comet ion tail. H 2 0 + images taken through narrow- 
band interference filters, such as those obtained by the Joint Observatory for 
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Cometary Research, will be extremely useful in separating the time variation from 
the velocity structure of the tail. An interesting feature appearing in almost all the 
datacubes was a two pronged tail at one wavelength and a single pronged tail at an- 
other wavelength, an indication that the outer regions of the ion tail may move with 
a different velocity than the inner regions (see figure 3.) The velocities of each image 
in figure 3 are in an arbitrary geocentric frame of reference. 
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Anita Cochran 

How long did it take to complete each data cube? 

Daniel Schultz 

Each image in a data cube took 5 to 15 minutes, with a whole sequence taking from one to 
two hours to complete. 


Paul Feldman 

Why are your values for Q(HjO) so high when Comet Austin was near the Sun? 

Daniel Schultz 

For these preliminary results, we estimated the aperture corrections. It is quite possible that 
our estimate was high for the first few points. We will analyze (01) imaging data to obtain 
more accurate and precise aperture correction. 
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Abstract 

We report results of spectrophotometric observations of comets P/Brorsen-Metcalf and Austin 
from 3000 to 3600A at a spectral resolution of about 1.8A, The strongest features are the OH 
(A*X) 0-0 and 1-1 bands, and the NH (A-X) 0-0 bands. We have clearly found, for the first time, 
the OH (A-X) 0-1 band. We verify the existence of, and quantitatively ■measure., the CN (B-X) 
2-1 and 3-2 bands, We detect a feature at 3258A that was first seen in uncalibrated spectra, and 
suggest identification as the NH singlet (c-a) 0-0 transition. We have also firmly identified the 
COj features at 3378, 3504, and 3512A, This ion was reported as being present in the tail of 
Comet Bester (1948 I) by Swings and Page (1950). We propose identification of a weak feature at 
3547A as the fundamental transition of H 2 CO, which would make this the first optical cometary 
detection of this molecule which is very abundant in giant molecular clouds. 

Introduction 

The near ultraviolet (3000-3600A) spectral region has not yet been fully exploited for the study of 
comets. Observations in this spectral region are difficult due to the high atmospheric extinction 
and the low throughput and sensitivity of spectrographs and detectors at these wavelengths. 
The first detailed study of the near-ujtraviolet was by Swings, Elvey, and Babcock (1941), who 
first detected the strong OH (A-X) 0-0 band at 3090A in Comet Cunningham (1941 I). One of 
the best ground-based ultraviolet spectra is that of Comet Kohoutek by Lane et al, (1975), The 
development of rocket and satellite observations has emphasized the importance of the OH bands, 
but with the exception of NH, the remainder of this spectral region remains largely unexplored. 
In this paper, we present spectra of comets P/Brorsen-Metcalf and Austin (1989cl) in the 3000- 
3600A spectral region. Detailed results for P/Brorsen-Metcalf have been submitted for publication 
in The Astrophyscial Journal. 

Observations 

All observations were made using the Large Cassegrain Spectrograph mounted on the 2.7-m tele- 
scope at McDonald Observatory. The detector was a thinned, back-side illuminated TI 800x800 
CCD. Prior to each of the observing runs, the CCD was soaked in dry NO gas at one atmosphere 
pressure and room temperature for about 15 minutes, in order to induce a negative charge layer 
on the back side of the chip. This process resulted in a 50-fold improvement in the response of 
the CCD at 3000A. The spectrograph was used in long slit mode, with a scale at the detector of 
1.28 arcsec per readout pixel, and a dispersion of 0.89A per pixel. A log of the observations of 
both comets is given below. 


Comet 

Date 

Exposure 

Rh 

A 

P/Brorsen-Metcalf 

30 July 1989 

15 min 

1.08 AU 

0.66 A U 

P/Brorsen-Metcalf 

10 August 1989 

15 min 

0.89 AU 

0.63 AU 

P/Brorsen-Metcalf 

10 August 1989 

60 min 

0.89 AU 

0.63 A U 

Austin (1989cl) 

30 May 1990 

109 min 

1.24 AU 

0.26 AU 

Austin (1989cl) 

2 June 1990 

65 min 

1.29 AU 

0.28 A U 


Separate exposures were made of spectrophotometric standard stars for absolute flux calibra- 
tion, For cometary spectroscopy, proper subtraction of the sky spectrum is essential. In August 
1989, a separate sky spectrum of 10 minutes was made after that of the comet. In July 1989 a 
set of planetary nebulae were observed before and after the comet. Since these objects did not fill 
the entrance slit length, the outer slit regions were used for determination of the sky spectrum. 
For the May 1990 observations, a separate sky exposure of 60 minutes was obtained. The average 
along the slit of the 10 August 1989 one hour spectrum of P/Brorsen-Metr^Jf is shown in Figure 1, 
and the central 15 pixels of the 30 May 1990 spectrum of Comet Austin in Figure 2. Note that 
the P/Brorsen-Metcalf spectrum has been calibrated to flux units, while the Austin spectrum is 
still in detector counts. 

Identification of Spectral Features 

OH. The well known OH (A-X) 0-0 and 1-1 bands, along with the Nil (A-X) 0-0 band, are 
the strongest transitions in the spectra. Our spectra are at sufficiently high spectral resolution to 
resolve the rotational structure of these bands. The relative intensities of the lines within the 0-0 
band should vary with heliocentric velocity due to the Swings effect. Our observed (0-0) line ratios 
are in good agreement with the theoretical values calculated by Schleicher (1983). The observed 
ratio of the integrated flux in the 1-1 to 0-0 OH bands is a good test of our understanding of 
the physics of cometary comae, Schleicher and A’Hearn’s (1988) prediction of the OH 1-0 to 0-0 
band ratio, together with their listed transition probabilities, leads us to expect a 1-1 to 0-0 band 
ratio of about 0.09. This value is only slightly outside the error bars of our measured value for 
P/Brorsen-Metcalf of 0.12±0.02. An independent measurement of this ratio in P/Brorsen-Metcalf 
obtained with the IUE at about the same time by P. Feldman (personal communication) gives a 
flux ratio of 0.082 ± 0.012. It is not yet obvious if there is any real discrepancy between theory 
and observations. 

In the 3450-3485A region we see an emission band which we identify as the first cometary 
detection of the v' = 0 to v" = 1 band of the (A-X) transition of OH, based on the wavelength 
identification from Dieke and Crosswhite (1962). The relative intensities of the lines are the same 
as those in the 0-0 band, since they arise from the same upper level. 

NH. The identification of the lines of the NH (A-X) 0-0 band at 3360A are in good agreement 
with the results of Kim, A’Hearn, and Cochran (1989). We see no evidence, however, to; the NH 
(A-X) 1-1 band, which should be shifted about 10.4A longward of the 0-0 band. 

We see a weak feature in the P/Brorsen-Metcalf spectrum at 3255A, which was also see in 
Kohoutek by Lane et al. (1975). This feature is confirmed in our comet Austin spectrum, along 
with an additional feature at 3266A. We reject an identification of these lines as the SII (A- 
X) electronic transition due to wavelength discrepancies with the observed feature. Instead, we 
propose the identification of the (c *n - a *A) 0-0 band of NH, According to Pearse (1934), the 
strongest features of this band would be the Q branch at 3255A and the P branch near 3267A. 
This would be the first identification of the NH singlet transitions in a comet. 
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Wavelength (Angstroms) 


Figure 1: The spectrum of comet P/Brorsen- Metcalf on 10 August 1989. This 60-minute exposure 
has been calibrated to flux units. The upper panel is scaled so that the entire OH (0-0) band is 
visible, while the lower panel shows the weak features between the strong OH and NH bands. 








Wavelength (Angstroms) 


Figure 2: The spectrum of comet Austin on 30 May 1990. This 109-minute exposure has not yet 
been calibrated to flux, so the units are CCD counts. Thus, the decrease atmospheric transmission 
and CCD response near 3000A serves to depress the strong OH (0-0) and (1-1) bands relative to 
the features at longer wavelength. The upper panel is scaled so that the entire spectrum is visible, 
while the lower panel shows the weak features between the strong OH and Nil bands. 








COj. We detect a broad emission feature at 3512A and a narrow feature at 3504A. These 
features are present in all of our spectra of both P/Brorsen-Metcaif and Austin. We identify 
these features as COj • The laboratory study of Mrozowski (1941) indicates that in a comet, the 
strongest COj features should be a blend of the first members of the R and P branches of the 0-0 
2 n^ u to 2 n^ at 3512A, and the same lines of the 0-0 2 II^ U to 2 JI| y at 3504A. The next strongest 

feature should be the R and P branches of the 1-0 2 II^ U to 2 II^ a at 3378A. Careful examination 

of the red end of the NH (A-X) 0-0 band at 3360A reveals a line at 3378A which can not be 
identified as NH. This line is most likely CO*. 

CN. A series of lines is detected at 3560-3590A, which we identify as the At = -f-1 band of 
the CN violet (B-X) electronic band. We find good agreement with the expected wavelengths 
of the P and R branches of the 1-0, 2-1, and 3-2 vibrational bands. This band was previously 
reported in Comet Cunningham (194i I) by Swings et al. (1941). 

H 2 CO. The spectrum of P/Brorsen- Metcalf shows a weak peak at 3547A, which is highly 
confined to the nucleus. This feature is also seen in the Comet Austin spectra, along with 
an additional peak at 353lA. We identify these as the ultraviolet electronic bands of H 2 CO. 
Laboratory data would lead us to expect a tight red-degraded band at 3546. 9A, along with an 
additional band at 353lA. Cometary H 2 CO has previously been detected in the IR and radio 
spectral region, but this would be the first optical detection of this molecule. 

Conclusions 

In this paper, we believe that we have established that the far ultraviolet end of the optical window 
is an excellent region for studying cometary physics. The cometary scattered continuum light from 
the Sun is low intensity and the major emission lines are grouped into narrow bands, in contrast 
to the visual and red regions which are dominated by scattered solar continuum and spread-out 
emissions such as C 2 and NH 2 . We have found evidence for several new molecular emissions, and 
have confirmed tentative reports of other species. We believe that improved spectra in the near 
future will result in the firm identification of these species and probably will reveal additional 
emissions. 

This research was supported by NASA Contract NAG5-29451 with Rice University and by 
NASA Grant NAGW-1477 and NSF Grant AST-8912869 to The University of Texas at Austin. 

References 

Dieke, G. H. and Crosswhite, H. M., 1962, J. Quant. Spec. Rad. Trans., 2, 97. 

Kim, S. J., A’Hearn, M. F., and Cochran, W. D., 1989, Icarus, 77, 98. 

Lane, A. L., Stockton, A. N., and Mies, F. H., 1975, in NASA SP-355, Comet Kohoutek, ed. 

G. A. Gary, (Washington, D. C.: NASA), p. 87. 

Mrozowski, S., 1941, Phys. Rev., 60, 730. 

Pearse, R. W. B., 1934, Proc. Roy. Soc. London A, 143, 112. 

Schleicher, D. G., 1983, PhD thesis, University of Maryland, College Park, Maryland. 

Schleicher, D. G. and A’Hearn, M. F., 1988, Ap. J., 331, 1058. 

Swings, P., Elvey, C. T., and Babcock, H. W., 1941, Ap. J., 94 , 320. 

Swings, P. and Page, T., 1950, Ap. J., Ill, 530. 


54 




N9 1 - 21023 
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Long slit ultraviolet spectra of Comet Austin (1989c,) were obtained from a 
sounding rocket experiment launched from White Sands Missile Range, New Mexico, 
on 21 April 1990 at 1015 UT. The instrument, known as the Faint Object Telescope, 
consisted of an 1715.5 Dall-Kirkham telescope, a Rowland Circle spectrograph and a 
microchannel plate intensifier coupled to a Reticon diode array. An onboard television 
camera transmitted images to the ground during flight, which permitted real-time 
maneuvers in order to center the comet in the entrance slit. The comet was held near 
the center of the slit for approximately 270 seconds by sending six pointing 
corrections, all but one of which was less than 30 arcseconds. The comet’s parameters 
at the time of observation, along with the characteristics of the instrument, are 
shown in Table I. 


Table I Comet Austin parameters for NASA 36.073 UG 


Heliocentric Distance 

0.476 AU 


Heliocentric Velocity 

31.44 km / sec 


Geocentric Distance 

0.756 AU 


Geocentric Velocity 

-37.63 km / sec 


Slit length at the comet 

8.1 arcmin 

(2.7 x 10 6 km) 

Slit width at the comet 

7.25 arcsec 

(4.0 x 10 3 km) 

Spectral Resolution 

5.3 A 


Wavelength Range 

1150 - 1850 A 


Time in slit 

T+250 to T+520 sec 


Pitch/Yaw maneuvers 

<, 30" 




The data taken by the instrument consists of the (x,y) coordinates of the 
photons which fall on the Reticon diode array. Figure 1 shows a contour 
representation of the raw data (including detector dark counts) obtained during the 
time that the comet was in the slit. In this representation, wavelength runs along the 
horizontal axis, from approximately 1150 A to 1850 A over the 800 pixels, while the 
vertical axis is the spatial 
direction, with a span of 
approximately 3.7 x 10 6 km 
across the 544 pixels shown; 
the slit is 400 pixels long. The 
coordinates shown are actual 
photon positions on the face of 
the Reticon diode array after 
being corrected for a tilt of the 
entrance slit. Because of the 
centroiding algorithm used, 
the position could be 
determined to an eighth of a 
Reticon diode, although the 
instrument resolution is closer 
to V£ of a diode. 

The three brightest lines in the spectrum, O I (1304 A) and C I (1561 A and 
1657 A) are clearly visible in this picture, while others, including S I (1807, 1820 and 
1826 A) and the CO fourth positive bands (1380 - 1600 A) can also seen in other 
graphical representations of the same data. A CaF 2 filter in front of the entrance slit 
attenuates all wavelengths below approximately 1240 A by several orders of 
magnitude, which prevents H I (1216 A) from being seen. The 15,808 counts in the 
spectrum shown include the detector dark count in addition to the photons from the 
comet. 



Figure 1 Contour plot representation of raw data from 
NASA 36.073 UG 


Figure 2 shows a spectrum obtained from the above data by summing over 4 
diodes in the spatial direction (40 arcseconds, or 21,500 km at the comet) around the 
center of brightness. Although the spatial resolution of the instrument is actually 10 
arcseconds, because of pointing jitter and drift the spatial data is smeared out here 
over approximately 30 arcseconds. It will be possible to remove much of this 
degradation in resolution by taking into account the actual pointing at each instant 
in time during the flight, but this has not been done yet. 

The data shown has been converted from raw counts to Rayleighs by folding 
in the instrument effective area at each wavelength. It has also been binned to V6 of 
a Reticon diode, had the detector dark counts subtracted, and then smoothed over 2 
bins (7 A). Only those wavelengths above 1250 A have been included due to the 
effects of the CaF 2 filter mentioned above. The wavelength scale was determined by 




Figure 2 Spectrum of comet Austin in central 21,500 km of the slit. 


measurement of the positions of the brightest lines seen. The slit is approximately 
l.S bins wide on the scale of this plot. 

The brightest spectral features are labeled in the figure, and their preliminary 
brightnesses are tabulated for several projected areas on the comet in Table II. The 
errors associated with these values are mostly due to the 20% uncertainty in the 
absolute calibration of the instrument, although for the CO bands, the counting 
statistics become a significant source of error. This table also includes a comparison 
to the brightnesses measured in the central coma of Comet Halley (1986 III) with a 
similar rocket-borne instrument on 26 February and 13 March 1986 (Woods et al. 
1986). The wavelength range did not extend above 1750 A in the 1986 flights; thus 
the S I lines were not observed. 

The 0 I (1304 A) multiplet is much weaker relative to Halley due to the lower 
heliocentric velocity which shifts the absorption wavelength of the cometary oxygen 
atoms outside the solar line profile (Dymond et al. 1989) and results in excitation 
through Bowen fluorescence via H I Lyman-fl (Feldman 1982). 

The measured CO brightnesses from the 40" column in Table II, when used in 
a simple radial outflow model, give a production rate of -1.7 x 10 27 s’ 1 , compared to 
2.4 x 10 29 s' 1 for the first Halley flight, and 1 x 10 29 s' 1 for the second. Assuming a 
water production rate of -1 x 10 29 s’ 1 , based on ground-based optical and radio 



Table II Brightnesses measured for Comet Austin compared with Comet Halley 




Comet Austin (1989c,) 

| Comet Halley (1986 HI) f 



21 April 
1990 

26 February 
1986 

13 March 
1986 



0.48 AU 

0.70 AU 

0.90 AU 



7.25" x 20" 

7.25" x 40” 

7.25" x 80“ 

10" x 16" | 

01 

1304 A 

80 ± 16 R 

100 * 20 R 

105 ± 20 R 

1460 ± 320 R 

280 * 65 R 

01 

1356 A 

40 ±8 R 

25 ±5 R 

10 ±2 R 

60 R 

< 10 R 

CO 

1478 A 

24 ± 12 R 

21 ± 11R 

- 

175 ± 50 R 

65 ± 25 R 

CO 

1510 A 

24 ± 12 R 

22 ± 11 R 

- 

140 ± 40 R 

60 ± 25 R 

Cl 

1561 A 

130 ± 25 R 

100 ± 20 R 

70 ± 15 R 

480 ± 110 R 

150 ± 35 R it 

CO 

1597 A 

33 ± 17 R 

15 ±8 R 

- 

90 ± 25 R 

30 ± 15 R | 

Cl 

1657 A 

400 ± 80 R 

360 ± 70 R 

300 ± 60 R 

1600 ± 350 R 

520 ± 115 R | 

SI 

1807 A 

- 

130 ± 25 R 

75 ± 15 R 

Not Observed | 


observations made within a few days of the rocket flight and the relative constancy 
of the visual magnitude during this period, this value indicates an abundance of CO 
relative to water of approximately 2 percent. 

This work was supported by NASA grant NAG 5-619. 
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IUE OBSERVATIONS OF COMETS P/BRORSEN-METCALF 

AND OKAZAKI-LEVY-RUDENKO £ '/¥ 

E. E. Roettger, P. D. Feldman, S. A. Budzien (Johns Hopkins University) ^ S 
M. F. A’Hearn (University of Maryland) 
and M. C. Festou (Observatoire de Besanson) 

Comet P/Brorsen-Metcalf (1989o) w as observed with the International Ultraviolet 
Explorer (IUE) on several dates from 27 July 1989 through 5 August 1989; comet Okazaki- 
Levy-Rudenko (1989r) was observed in low dispersion with IUE on 2 and 22 December 
1989. The water production rates for P/Brorsen-Metcalf derived from OH measurements 
using the vectorial model of Festou (1981) were 4.2, 4.6, 5.3, 4.9, 6.2, and 6.7xl0 28 on 27.2 
and 31.1 July, 1.8, 2.4, 3.9, and 5.2 August 1989, respectively. These are plotted in Fig> 1. 

For perspective, water production rates estimated from postperihelion visual magnitudes 
using the method of Newburn (1981; see also Roettger et a/., 1990) are also plotted. Fig. 

2 shows the water production rates (squares) for comet Okazaki-Levy-Rudenko, which 

were 5.3 and 4.0xl0 28 on 21.9 and 22.2 December 1989. The 21 December measurement 
has a greater uncertainty, as it was derived from the OH (1,1) band, in a noisier region 
of the IUE camera. Results using a newer version code for the vectorial model are also 
plotted (diamonds); they are within a few percent of the previous results. The numbers 
listed are from the older code, to maintain consistency with published results for other 
comets. Also plotted in Fig. 2 are two estimates of the water production rate from the Ol 
emission around 1304 A (octagons).. The higher point is from a Haser model, and should 
be an overestimation as oxygen is primarily a granddaughter product. The lower point 
was obtained by analogy with observations of comet Bradfield 1979X, a comet with similar 
IUE spectra and geometric parameters. These should be considered rough estimates only. 

The ultraviolet continuum was not reliably measurable in comet P/Brorsen-Metcalf 
in most of the spectra, but we derive Afp ~ 30 cm for 2 August. For comet Okazaki-Levy- 
Rudenko, we derive Afp = 50 cm on 22 December. Afp is the product of albedo, filling 
factor of grains in the aperture, and projected aperture radius (A’Hearn et al., 1984); it is 
proportional to the dust production rate if the production is steady-state, and the grains 
move at constant velocity without fragmentation or albedo changes. Both comets had very 
little dust compared to other recently observed comets, as can be seen in Fig. 3. 

r • . 

The ratio of CS to water production rates are shown in Fig. 4. Relative to water, both 
comets had less CS than P/Halley in the same range of heliocentric distances. The CS 
production rales were obtained from a Haser model, assuming a short-lived parent (CS 2 ). 

For model-independent comparison, the column densities of observed species are listed in 
T»bi<e- 1. Because the dynamic range of the IUE spectrographs are limited, not all species 
can be measured in each spectrum. The P/Brorsen-Metcalf spectra from 31 July and 3 
August were obtained from the NSSDC Archives. Due to the observing technique used for 


these spectra, the exposure time and hence the brightnesses and derived quantities have a 
small uncertainty associated with them* 

Mention was made previously (see O'Dell et a/., this volume) of a possible inconsis- 
tency of the OH y-factors from Schleicher and A’Hearn (1988) and the OH (1,1) to (0,0) 
band ratios. The (1,1) band is in a noisy region of the IUE camera, and so those measure- 
ments generally have a greater uncertainty than the OH (0,0) or (1,0) band measurements. 
Although individual band ratios vary from their theoretical values by typically 15%, the 
variation is in both directions and the average of the theoretical to measured ratios differ 
from the ideal value of 1.00 by only 2%, This is easily within the uncertainty of either 
the y-factors or the WE. The nucleus and offset exposures were then averaged separately; 
the deviations from theory are higher but still within the expected uncertainties. These 
results and related calculations for the OH (1,0) band are presented in Table II. 

Acknowledge me nt s 

This work was based on observations obtained at GSFC and Vilspa and supported 
by NASA under Grants NSG-5393 to the Johns Hopkins University and NAG 5-252 to 
the University of Maryland. Additional data courtesy of the WE Science Working Group, 
National Space Science Distribution Center (NSSDC). 

References 

A’Hearn, M. F., D. G. Schleicher, P. D. Feldman, R. L. Millis, and D. T. Thompson 1984* 
Astr. J. 89 , 579. 

Festou, M. C. 1981. Astr. Ap. 95, 69. 

Newburn, R. L, Jr. 1981. In The Comet Halley Dust and Gas Environment Workshop 
(ESA SP-174) (The Netherlands: ESA Publications Division), p. 3. 

O’Dell, C. R., W. D. Cochran, C. O. Miller, A. L. Cochran, C. B. Opal, and E. S. Barker 
1990. This volume. 

Roettger, E. E., P. D. Feldman, M. F. A’Hearn, and M, C. Festou 1990. Icarus, in press. 
Schleicher, D. G., and M. F. A’Hearn 1988. Ap. J. 331, 1058. 


Table L Column densities (9 M xl6 M aperture) 


P/Brorsen-Metcalf Q-L-R 


27 Jul 31 Jul 

1 Aug 

2 Aug 

3 Aug 

5 Aug 

2 Dec 

21 Dec 

22 Dec 

r(AU) 

1,12 

1.06 

1.03 

1.02 

0.99 

0.97 

0.78 

1.04 

1.04 

A (AU) 

0,69 

0.65 

0.64 

0.64 

0.63 

0.63 

0.52 

0.81 

0.81 

f (km 8 *) 

-29.4 

•29.6 

-29.7 

-29.7 

-29.8 

-29.8 

+20.0 

+25.5 

+25.6 

from 

GSFC NSSDCA Vil»p« 

GSFC 

NSSDCA 

GSFC 

GSFC GSFC 

GSFC 

B,///(G,0) (kR) 

11,2 


20,0 

19.3 


31.8 



32.5 

B„//(1,0) (R) 


740 



1250 





B f ,//(i, i) (R) 








1380 


Now (10 13 ) 

3.0 

3.7 

4.4 

4.2 

5.5 

6.3 


4.2 

3.2 

B cs (R) 

150 

210 

(240) 

270 

300 

360 

930 

150 

190 

N C f (10 11 ) 

2.8 

3.4 

3.6 

4.0 

4,2 

4.9 

8.1 

2,3 

2.9 

Hcs/Boh (10- 3 ) 

9.2 

9.3 

8.2 

8.6 

7.6 

7.8 


5.4 

9.3 

Bco^ (R) 

50 

no 

(60) 

70 

100 

120 

280 

100 

120 

N cot ( 10U ) 

2.9 

4.7 

2.9 

3,5 

4.8 

5.4 

2.2 

4,9 

6.2 

* co* Wok ( 10 “ 3 ) 

9.6 

13. 

6.6 

8.3 

8.6 

8.6 


12 

20 

Bo/ (R) 

17 


33 




370 



No flO 13 ) 

1.7 


2.9 




2,7 



No /Non 

0.58 


0.65 







Be/ (R) 

32 


36 




260 



lie (10 ,Z ) 

3.9 


3.7 




4.3 



Nc/Noh 

0.13 


0.084 







Bs/ (R) 

47 

60 

72 




310 



Ns (10 12 ) 

1.4 

1.6 

1.8 




5.0 



Ns/No// (10' 2 ) 

4.6 

4.2 

4.0 







Be, (R) 







250 



Nc- (10 u ) 







4.2 




Table II. 

OH Ratios for Comet P/Brorsen-Metcalf 1989o 



9u/9i hi Bu/B,m B/g' 


9ii/9m B\ 

11/Bin 

B/g' 

9n(i/9l't 

Bi)i>/B 10 

B/g" 


Avg0.0816 0.0803 0.980 


1.75 

1.84 

1.05 

21.33 

21.29 

0.998 


dev 0,001 0,126 0.151 


0.00 

0.30 

0.173 

0.19 

2.31 

0.106 


off by 

2 . 0 % 




5.4% 



0.2% 





NUCLEUS exposures only 





Avg 

1,028 




1.12 



0.997 


dev 

0.173 




0.17 



0.041 


off by 

2.8% 




12% 



0.3% 





OFFSET exposures only 





Avg 

0.923 




0.966 



0.999 


dev 

0.112 




0.127 



0.160 


off by 

7.7% 




4.4% 



0.2% 



‘ The ratios (B|i/B« M i)/(pii/^oo), (£ii/Rm)/(9ii/<7in), and (Ba n / Bw)/{gw/gm) were averaged after divi- 
sion. 


61 



0.5 0.6 0.7 0.6 0.9 1 


HELIOCENTRIC DISTANCE (AU) 

Figure 1 . Preperihelion water production rate* derived from OH emissions (pierced symbols) for comet 
P/Brorsen-Metcalf are plotted with production rates derived from visual magnitude estimates. 
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Figure 2. For comet O-L-R, postperihelion water production rates derived from OH emissions using older 
(squares) and newer (diamonds) vectorial codes are plotted with rates derived from visual magnitudes. Two 
estimates of the water production rate from the Ol feature at 1304 A are also plotted (hexagons); the upper 
is from a Haser model, the lower from comparison with comet Bradfield 1979X. 
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Figure 9. The ga*-to-dust ratios of both comet P/Brorsen Metcalf and comet Okasaki-Levy-Rudenko 
are much higher than those of other recent comets -4/p, a quantity approximately proportional to the dust 
production rate, is derived from reflected solar continuum measurements 


B P/BRORSEN -METCALF 
0 OKAZAKI -LEVY -RUDENKO 

♦ P/TEMPEL 2 

♦ WILSON 
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HELIOCENTRIC DISTANCE (AU) 

Figure 4. The relative CS production rate was derived from CS emissions around 2576 A using a Haser 
model. Comets Brorsen- Metcalf and O-L-R both have lower Qcs IQ ratios than comet Hailev at similar 
heliocentric distances. 
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IUE OBSERVATIONS OP COMET AUSTIN (1989ci) 


Observations of comet Austin (1989ci ) were made using the International Ultraviolet 
Explorer satellite on December 22, 1989 (preperihelion) and on nine dates between May 7 and 
24, 1990 (postperihelion). A summary of observation dates and parameters is given in Table I. 
The ultraviolet spectrum is characteristic of that of non-dusty comets. Here . we present some 
preliminary results for OH and CS prodution rates and dust-reflected solar continuum emission 
based upon the long wavelength, low resolution mode spectra; photometry using the IUE Fine 
Error Sensor; a detection of the CO fourth positive system in the short wavlength, low dispersion 
spectra; and an upper limit on S 2 . 

Figure 1 depicts the H 2 O production rates both pre- and postperihelion, based upon 
OH solar resonance fluorescence and a vectorial model (Festou, 1981) which was circulated by 
Festou in 1988. The preperihelion water production rate was 1.4 x 1O 20 s -1 at a heliocentric 
distance of 2.2 AU, and the postperihelion values are 1.0, 0.96, 0.85, 0.78, and 0.60 x 10 29 for 
May 9.4, 13.7, 16.8, 19.7, and 24.3, 1990 at heliocentric distances from 0.80-1.13 AU. Since 
the longer exposures were saturated in the OH (0-0) band, for some dates the H 2 O production 
rate was derived from the OH (1-0) band, scaling by the observed (0— 0)/(l— 0) band ratio. 
The post -perihelion production rate roughly follows an r -1,6 power law during the period of 
observations. Production rates from visible photometry based on a Haser model (Schleicher, 
1990) are also plotted, scaled by 1.7 to approximate Vectorial model values. 

Figure 2 shows the behavior of the H 2 O and CS production rates, the Afp dust 
parameter at 2950 A, and the FES photometry from May 7-24. The CS prodution rate was 
derived using a Haser model (Haser, 1957), because the parent CS 2 has a dissociative lifetime 
(~ 500 s) chort compared to CS (~ 10 5 s). The CS column densities were based upon either the 
CS (0-0) or CS (0-1) band, depending upon exposure length, using the CS (0-0) g-factor from 
Jackson et al., 1982. The FES is sensitive to wavelengths 380-650 run, which include emissions 
of CN, C 2 , and continuum, and has an effective aperture approximately 18 arcseconds square. 
The FES counts follow an r' -3-8 power law when normalized by heliocentric distance squared and 
geocentric distance to scale as a production rate. The Afp parameter, which is the product of 
grain albedo, aperture filling factor, and projected aperture radius, assumes an r -2 dust density 
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distribution and is proportional to a dust grain production rate. The flattening of Afp at larger 
heliocentric distances may be due to a combination of the decreasing geocentric distance and 
a dust distribution which falls off more sharply than r~ 2 . Using the Qmo/Afp as a measure 
of the gas-to-dust ratio, we find Austin at 1.0 x 10 27 s~ 1 cnT' 1 on May 16.8 is comparable to 
that of Okazaki-Levy-Rudenko and about 15 times higher than Halley at similar heliocentric 
distances (Roettger ei a/., 1990). 

Finally, we used multiple linear regression and the ratio-of- variances methods to fit 
synthetic spectra to observed spectra to determine CO and S 2 abundances. The CO fourth pos- 
itive bands were observed in several short wavelength exposures with various signal-to-noise lev- 
els. From the May 9.4 exposure (Figure 3) we derive a CO production rate of 2.0:t0.9x 10 27 s~ l , 
or 2% the water production rate, based on a Haser model, and the solar maximum g-factors 
of Durrance (1981). From a 440 minute, long wavelength exposure on May 24.5, we obtained 
an S 2 production rate upper limit of upper limit of 2 x 10 24 s* 1 , or 3 x 10“ 5 that of water. 
For comparison, this abundance is an order of magnitude below the S 2 abundance of comet 
IRAS-Araki-Alcock ( 1983 VII) when the species was discovered (A ’Hearn ct al . , 1983). 
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We have been engaged in an observational search for cometary parent molecules 
using infrared spectroscopy in the 1-5 pm region. The investigation involved two 
different observing programs, one at moderate spectral resolution and one at fairly high 
resolution. 

The lower resolution program has been used to study cometary spectra in the 
vicinity of 3.5 pm at k/A^lO 3 . Comets P/Brorsen-Metcalf (1989o), Okazaki-Levy- 
Rudenko (I989r) and Austin (1990cl) have been observed with the Cryogenic Spec- 
trometer (CRSP) at Kitt Peak. The detector incorporated an InSb array with 58 spatial 
elements, each 2.7" on the sky, and 62 spectral channels per spatial element. 

An, as yet, unidentified feature has been detected at approximately 3.52 pm in 
Comet Austin (on 1990 May 4, 5 and 6). The feature is possibly present in 
P/Brorsen-Metcalf (observed on 1989 August 23 and 25), as well. Comet Okazaki- 
Levy-Rudenko exhibited continuum emission only in this spectral region at the time of 
our-observations (1989 November 14 and 16). The data will be presented, and the 
relationship between the 3.52 pm feature and cometary activity (e.g. water production 
rate, visibility of the 3.4 pm emission feature) will be discussed. 

The high-resolution (A./AX~2xl0 4 ) program probed comet Austin in the 4.8 pm 
region. These observations, conducted UT 1990 May 14-17 with the University of 
Texas Infrared Echelle Spectrometer (IRSHELL) on the IRTF atop Mauna Kea, were 
used to search for emission lines comprising the (1-0) vibration-rotation band of the 
ground electronic state of CO. 

The instrument used a Si-As impurity band array detector with 10 spatial ele- 
ments, each 1" on the sky, and 64 spectral channels per spatial element. Two grating 
settings were employed. One was centered between the P2, P3 lines and the other 
between R2, R3. Retrieval of these lines allows a probe of the population distribution 
of levels J'=l through 4 of the excited (v'=l) vibrational state within the ground elec- 
tronic state of CO. Knowledge of this distribution can be used to constrain the rota- 
tional temperature. Preliminary analysis suggests the P3 line was present UT May 16 
at roughly the 5a level. Results concerning the existence of other lines, and physical 
conditions inferred therefrom will be discussed. 
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I. Introduction 

Of the periodic comets, only comets Halley and Brorsen-Metcalf have periods of the 
order of 75 years and are bright enough to study spectroscopically in the infrared. In this 
paper we present thermal IR spectroscopy of comet p/Brorsen-Metcalf near perihelion and 
compare the results to comet p/Halley. 


II. Observations 

Spectra were obtained on comet p/Brorsen-Metcalf in the 8 to 13 Jim region with a 
resolving power SR = X/AX of about 50 using the 3.0 meter Infrared Telescope Facility on 
Mauna Kea with the liquid He-cooled facility circular variable filter (Tokunaga 1986) on 
August 28, 29, 30 and Sept. 4, 5, and 6, 1989 UT (14-6 days before perihelion). North- 
south chopper throws of 20 arc sec were used with a 7 arc sec aperture (5050 km at the 
comet) in August and 60 arc sec N-S chop with a 9.5 arc sec aperture in September (6860 
km at the comet). 


ILL Results 


In all, six 10 pm spectra of comet Brorsen-Metcalf were obtained. Within the 3-8% 
uncertainties, there was no variation in the shape of the spectra. Figure 1 shows one 
spectrum obtained on September 5, 1990 UT (R = 0.50, A = 1.00). The most obvious 
aspect of the spectra is their smoothness. Noticeably absent from the spectrum of Comet 
p/Brorsen-Metcalf in Figure 1 is any evidence for silicate emission. Indeed, no spectral 
feature was present at the 5% level. Also shown is a greybody fit to the data, which gives 
T = 399 K, consistent with the temperatures obtained by IR photometry (Lynch et al. 

1990). 


Figure 2 shows the spectrum of comet p/Halley obtained by Campins and Ryan 
(1989) taken with the facility CVF at the IRTF on Jan 16.08 UT (R = 0.79, A = 1.42) 
with an 8,2 arc sec aperture. Comet Halley displayed a strong, structured emission feature 
that Campins and Ryan identified with crystalline olivine. 


IV. Discussion 

In contrast to comet p/Halley which displayed significant photometric outbursts and 
brightness changes (Tokunaga et al. 1986; Lynch and Russell 1988), comet p/Brorsen- 
Metcalf showed little photometric variability. We note, however, that estimates of 
variability depend on the regularity and duration of observation, Halley was scrutinized 
closely for over a year whereas comet B-M was only observed episodically. 

Like virtually ail comets, both Halley and Brorsen-Metcalf showed IR color 
temperatures that exceeded the equilibrium by ~5% and -20% respectively. The 
enhancement in temperature must occur via a small CL, in the JR and the usual assumption 
involves small grains. Grains small enough to show a significant temperature excess are, 
almost by definition, optically thin and will therefore show any band structure present in 
their emissivity curves (as opposed to optically thick particles which cannot show band 
structure in comets). If we adopt the small grain hypothesis as the origin of the 
temperature excess, then based on the featureless 8-13 pm spectrum and the detection of 
the 3.4-3.6 pm emission features by Tokunaga arid Brooke (private communication) which 
are usually associated with carbon-rich species, then small carbon grains are the most likely 
candidates for the cometary dust population in Comet Brorsen-Metcalf. In contrast, the 
dust in comet p/Halley exhibited strong silicate emission, usually associated with oxygen- 
rich dust, as well as highly variable emission in the IR. 


V. Summary and Conclusions 

Except for their absolute brightness and similar orbits, comets p/Brorsen-Metcalf and 
p/Halley could hardly be more dissimilar. This is especially true of their IR spectra: 

Halley showed strong structured silicate emission, B-M showed no silicate emission. 

Halley was highly variable, B-M displayed no outburst activity or any variation which 
departed from a brightness variation that was monotonic with heliocentric distance. Indeed, 
if we were to ask what comets, periodic or otherwise, most resembled these two comets, 



we would have to say that IRAS-Araki-Alcock (Hanner et ah 1985a; Feierberg et al. 1984), 
Grigg-Skjellerup (Hanner et al, 1984) and Churyumov-Gerasimenko (Hanner et al. 1985b) 
were comparable to B-M, while p/Bradfield 1987s (Lynch et al. 1989; Hanner et al, 1990; 
Russell et al. 1990) was most like p/Halley (Campins and Ryan 1989). The marked 
dissimilarities between the two comets suggests that the periodic comets do qq! comprise 
one class, but rather have very different histories. 
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Walter Huebner 

How well do the laboratory data of amorphous and crystalline olivine fit the Comet Bradfield 
observations? How does one produce crystalline olivine? Are there other materials, particu- 
larly non-silicates that could also explain the observed spectrum? 

David Lynch 

The laboratory data, combined with particle size distributions (which lead to temperature 
distribution) fit the observed data reasonably well. How does one produce crystalline olivine? 
There are lots of ways to do that. I think the real question is why should crystalline olivine 
dominate over amorphous olivine. A match in spectral shapes between laboratory and comet 
data does not nail down an identification; it only suggests one. Given the many processes that 
can alter the crystalline structure in comets (Cosmic rays, UV, chemistry, thermal cycling) I 
am amazed that the observers (including myself!) have been so eager to make mineralogical 
identification. There arc many minerals that give plausible fits, mostly silicates. I would like 
to note that the best fit I could find for the 12.2 ft m feature in Comet Wilson was a calcium 
tungstate! 


Walter Huebner 

How does an amorphous silicate become crystalline? 

Comment by Ray Russell: 

If grains of amorphous silicate are heated they can undergo a phase transformation to the 
crystalline form. Laboratory experiments with heating above 500 K, specifically to 700 K for 
hours, turns the amorphous material to crystalline form. There is a range of temperatures 
over different time periods that will accomplish this (i.e., 500 K for days would be expected 
to do a similar job as 700 K for hours). This has not been well documented, i.e., we have not 
made measurements of time required versus temperature. 

David Lynch 

We think this mechanism could happen for periodic comets whose volatiles would be cleaned 
from the surface, leaving a crust that can be heated repeatedly. 
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Abstract 


The 1667 MHz and 1665 MHz transitions of the OH radical were recently monitored in several 
comets with the Nangay radio telescope: P/Brorsen-Metcalf (1989o) (August 4 to October 31, 
1989), Okazaki-Levy-Rudenko (1989r) (October 3 to December 2, 1989), Aarseth-Brewington 
(1989al) (December 8 to 30, 1989) and Austin (1989cl) (February 15 to June 14, 1990). We 
present gas expansion velocities measurements obtained from the analysis of the line shapes and 
discuss the long term variations of the water production rate, as measured from the OH radio 
lines. On October 13, 1989, we fortuitously* observed the occultation of a. background source by 
comet Okazaki-Levy-Rudenko (1989r). A preliminary report of this observation is given. Further 
analysis of the data is continuing. 

Introduction 

Thanks to their high spectral resolution, the observations of the OH 18-cm radio lines are a 
powerful tool to probe the neutral gas expansion velocity in the outer coma. This parameter is a 
requisite to the conversion of cometaiy molecular signals into production rates and is fundamental 
for the understanding of the thermodynamic of the coma. In order to separate the contribution of 
the OH ejection velocity v d from that of the OH-parent velocity v p in the OH radio profiles, 
Bockelee-Morvan et al« (1990) developed a deconvolution method, which consists in fitting 
trapezia to the observed spectra. The H 2 0 velocity is then derived from the half large base B| of 
the best trapezium, which is assumed to be close to Vp-t-v^. 

Gas expansion velocities 

The trapezoidal method was applied to several comets, including P/Halley, observed at Nangay 
over a variety of heliocentric distances and production rates. The inferred v p were found to vary in 
the range 0.5 to 2 km s' 1 and to depend, as expected from thermodynamical models, both on rh 
and on the gas production rate (Bockelee-Morvan et al., 1990). 

As a complementary study, trapezoidal fitting has been performed to some selected integrated 
spectra of comets P/Brorsen-Metcalf (1989o), Okazaki-Levy-Rudenko (1989r), Aarseth- 
Brewington (1989al) and Austin (1989cl). Fig. 1 shows the spectrum of comet Austin observed 
from May 21 to 25, 1990 with the fitted trapezium. The results of the fits, i.e. the half larger base 
B| and the half smaller base B s of the best trapezia with their 1 -sigma formal errors, are given in 
Table 1, together with the observational parameters. Outflow velocities can be readily deduced by 
subtracting v d = 0.9 km s' 1 from Bj, according to our best observational estimate of the OH 
ejection velocity (one should note that the water photodissociation model of Crovisier (1989) 
predicts a somewhat larger value of 1 .05 km s' 1 ). 

The observations performed at small rh ( < 1 AU) confirm the global increase of the outflow 


velocity with decreasing heliocentric distances (Table 1). Outflow velocities as large as 1.6 km s' 1 
, . ijre found at 0.5 AU for production rates on the order of 10 29 s'*. On the other hand, when 
H observed at large heliocentric distances, comets with similar production rates exhibit outflow 
velocity always below 1 km s’ 1 . The observations performed near rh = 1 AU are in agreement 
with the general trend of the outflow velocity with increasing production rate for constant rh 
found by Bockelee-Morvan et al. (1990) : at rh * 1 AU and below Q(OH) = 5 10 28 s' 1 , v p is in 
the range of 0.7-0.8 km s* 1 ; a steep increase of the outflow velocity with increasing production 
rate then occurs at Q(OH) = 10 29 s' 1 , reaching values larger than 1 km s'* and on the order of 1.8 
km s* 1 for Q(OH) = 6 10 29 s'*. For large rh (> 1 AU) and small Q(OH) (< 5 10 28 s’ 1 ), most 
inferred HjO velocities are below 0.8 km s' 1 . 


Austin 1989c! 900521-900525 



OH production rates 


As detailed by Gerard (1990) and Crovisier and Schloerb (1990), the determination of the OH- 
parent production rate from the observed OH 18-cm flux density (Table 1) heavily relies on the 
density and excitation models assumed. We use here the model of Bockelee-Morvan et al. (1990), 
which brought radio OH production rates of comet Halley into reasonable agreement with UV- 
derived values. This model uses the Haser-equivalent formalism, in which the H 2 0 expansion 
velocities are those derived by trapezoidal fitting. In addition, the quenching of the OH maser 
lines by collisions is taken into account, using the quenching radius measured by Gerard (1990) 
in comet P/Halley and scaling the value for the production rate and heliocentric distance as 
discussed in Schloerb (1988). For comparison, OH production rates calculated with the set of 
Haser-equivalent parameters referred to as 1986a, assuming no quenching effect, are also given in 
Table 1. Both set of production rates are computed using the OH 18-cm maser inversion curve of 
Despois et al. (1981). 

From the radio measurements, the OH production rate in P/Brorsen-Metcalf varied little 
through August 1989, from about 1 to 2 10 29 s' 1 . Roettger et al. (this workshop) report water 
production rates between 4.2 10 28 and 6.7 10 28 s* 1 from July 27 to August 5, 1989, while 
DiSanti et al. (this workshop) measured about 9 10 28 s' 1 on August 8, from observations of the 
Oj 'D line. For comet Okazaki-Levy-Rudenko, we infer production rates between about 2 and 7 
10 28 s' 1 as the comet approaches the Sun from 1 to 0. 65 AU. For the pre-perihelion observations 
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Figure 2: Geometry of the 
occuitation. The origin of 
coordinates is the comet 
centre. The thick kine shows 
the relative position of the 
background source between 
the beginning (B) and end 
(E) of the observation. 
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Figure 3: The OH spectrum observed on October 13.54 (dotted line) superimposed on the 
comparison spectrum averaged over the two preceding and two following days (full line). 
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of comet Aarseth-Brewington at rh s* 0.5 AU, we derive an OH production rate of 2.2 10 29 s* 1 . 
From Table 1, it appears that the OH production rate in comet Austin was higher before perihelion 
than after, reaching 2 10 29 s'* at rh = 0.38 AU. A similar trend is present in the comet heliocentric 
magnitude. IUE OH observations lead to somewhat larger values than ours for the production rate 
of the end of May 1990 ( 6 10 2i * s' 1 on May 24, Budzien et al., this workshop) 

Comet Austin was hardly detected from April 6 to May 15 1990, mainly because the maser 
OH inversion rate was not favourable at that time. However, the observations performed during 
that period gave interesting results on the excitation mechanism of the cometary OH molecule: 
Schleicher and A’Heam (1988) predict a change from absorption to emission at the heliocentric 
velocity of 35.1 km s"’, while the model of Despois et al. (1981) predicts the turn-over at vh = 
36.5 km s' 1 . The observations from April 28 to May 9, which cover heliocentric velocities 
between 35.3 and 35.6 km s 1 , reveal a weak line in absorption and are in favour of Despois et 
al.'s model. On May 31 and June 1, comet Austin passed across the galactic plane, leading to an 
amplification by a factor of about 2 of its OH maser radio line. The exact amplification factor to be 
used for retrieving the OH production rate is still in study. Thus, in Table 1, we omit the 
production rate result. 

Oecultation 

On October 13, 1989, we fortuitously observed the oecultation of a background source by comet 
Okazaki-Levy-Rudenko (1989r) (Bockelee-Morvan et al. 1989). B2 1426+294, a point source of 
0.4 Jy at 18 cm identified with a blue stellar object, passed at about 1' from the comet centre 
(corresponding to 60 000 km). The geometry of the event is shown in Fig. 2. The observed lines 
were enhanced by a factor of two with respect to the average of the spectra obtained before and 
after the oecultation (Fig. 3). This is in agreement with the maser theory of the OH line excitation. 
There is a marginal evidence of a variation of the signal intensity during the oecultation, while the 
line-of-sight to the source swept a path of about 50 000 km in the comet frame, which suggests 
that the OH distribution (or the excitation of the OH radicals) is not homogeneous within the coma 
(Crovisier et al., in preparation). The observation of such occultations is potentially very 
promising to investigate the coma structure, but such events are rare (only one is expected every 
several ten years at Nan?ay). 
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33 MILLIMETRE OBSERVATIONS OF COMETS P/BRORSEN-METCALF (!989o) AND 
AUSTIN (1989cl) WITH THE IRAM 30-in RADIO TELESCOPE. 

P. Colom (1,2), D. Despois (3), D. Bockelee-Morvan (4), J. Crovisicr (4) and G. Paubcrt (1) 


(1) IRAM, Avenida Divina Pastora, 7, N.C., E-18012 Granada, Spain 

(2) Centro Astronomico de Yebes, Apartado 148, E- 19080 Guadalajara, Spain 

(3) Observatoirc de Bordeaux, BP 89, Avenue Pierre Semirot, F-33270, Floirac, France 

(4) Observatoire de Paris, Section de Meudon, F-92195 Meudon, France 


Abstract 

Millimetre observations with the IRAM 30-m telescope were conducted in comet P/Brorsen- 
Metcalf (1989o) on September 1989 and Austin (1989cl) on April and May 1990. The HCN J(l- 
0) and J(3-2) lines were detected in both comets. The HCN production rate relative to water in 
P/Brorsen-Metcalf is comparable to that previously measured in comet P/Halley, while that 
inferred in comet Austin might be smaller by a factor of two. The H 2 CO 3 12 -2| t transition, 
marginally observed in comet P/Brorsen-Metcalf, was firmly detected in May 1990 in comet 
Austin. Observations performed at offset positions suggest that the source of H 2 CO might be 
distributed. The H 2 CO abundance is on the order of 0.5% that of water for both comets, assuming 
a scalelength of 10 4 km at 1 AU from the Sun for the distributed source. During the May 
observing period of comet Austin, two new species were detected for the first time in a comet: 
hydrogen sulfide H 2 S through its 1 10 - lot ort bo line at 169 GHz, and methanol CH 3 OH through 
J(3-2) aK = 0 transitions at 145 GHz. Preliminary estimates of their abundances are 1.5 10 3 for 
H 2 S and 8 10 3 forCH 3 OH. 

Observations 

The observations were performed with the IRAM 30-m millimetre telescope on Pico Veleta in 
Spain. Comet P/Brorsen-Metcalf (1989o) was observed on September 2-7 1989 at rh = 0.51 AU 
from the Sun and A = 1 AU from the Earth. Comet Austin (1989cl) was first observed on April 2 
1990 (close to its perihelion) at rh = 0.41 AU and 'A = 1.19 AU, and secondly in the May 21-25, 
1990 period (around its perigee) at rh = 1.1 AU and A = 0.24 AU. We used four receivers in the 
double side band: 1. a 83-1 15 GHz SIS mixer receiver; 2. a 130-170 GHz SIS mixer receiver, 3. a 
209-270 GHz SIS mixer receiver; 4. a 330-345 GHz cooled Schottky receiver. It was possible to 
use simultaneously 3 receivers, which allowed us to search for different lines at the same time. 
The backends consisted in two banks of 128x100 kHz channels (corresponding to a velocity 
resolution of 0.34 km s' 1 and 0.11 km s* 1 for the HCN lines at 89 GHz and 266 GHz 
respectively), and two banks of 512x1 MHz channels. In order to study the extension of the coma, 
observations at 6 arcsec and 12 arcsec North, East, South and West from the nucleus were 
performed in comet Austin during the May observing period. The observed integrated line areas or 
3-sigma upper limits are shown in Table 1 . 

Hydrogen cyanide 

The HCN J(l-0) (89 GHz) and J(3-2) (266 GHz) transitions were marginally detected in comet 
P/Brorsen-Metcalf. The HCN J(l-Q) transition was only detected previously in comet Kohoutek 
1973 XII and P/Halley despite many searches in other comets (Crovisier and Schloerb, 1990). 
Comet P/Brorsen-Metcalf is the first comet in which the HCN J(3-2) transition has been detected. 
On April 2, the HCN J(l-0) line was not detected in comet Austin. However, during the May 


Figure I: The HCN J( 1-0) 
spectrum of comet Austin 
observed on May 21 to 23, 
1990 (centre and offset 
positions). Velocity scale is 
with respect to the comet rest 
velocity, for the F 2-1 main 
hypeifine component, 



Figure 2: The HCN J(3-2) 
spectrum of comet Austin 
observed on May 23-24, 
1990 (centre and offset 
positions), 
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Figure 3: H 2 CO 3 12 -2j 1 
spectrum of comet Austin 
observed on May 21, 22, 25 
1990 (centre and offset 
positions). 
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observing session, both HCN J (1-0) and J(3-2) transitions were firmly observed in the comet 
(Bockelee-Morvan et al., 1990b; Crovisier et al., 1990a): some selected integrated spectra are 
shown in Figs 1 and 2. 

In order to interpret the present observations in terms of HCN production rates, we used the 
observed J(l-O) intensities and the model of Crovisier (1987), The HCN abundance relative to 
water in P/Brorsen-Metcalf is in agreement with that of ~ 10 3 measured in P/Halley (Table 1), but 
the HCN abundance in comet Austin is smaller by about a factor of 2. A similar deficiency was 
found in comet Wilson (1987 VII) (Crovisier et al„ 1990b), suggesting a chemical difference in 
periodic and new comets. 

Formaldehyde 

The observations of the H 2 CO 3 12 -2, , transition at 226 GHz in comet P/Brorsen-Metcalf lead 
only to a marginal detection, with a signal-to-noise ratio of 4. However, it gave the first conclusive 
evidence for the possibility of observing cometary formaldehyde at millimetre wavelengths. 
Formadehyde was not detected in comet Austin in April. It was easily detected at 226 GHz in May 
(Bockelee-Morvan et al., 1990b), with a S/N ratio larger than 10 for the average spectrum 
obtained for May 21, 22, 25 shown in Fig. 3. 

It has been suggested that formaldehyde gas originates from a distributed source in the coma 
(Snyder et al. 1989) which could be polyoxymethylene (POM) associated with the dust (Bojce et 
al., 1990). Our observations of the 226 GHz H 2 CO line performed at offset positions also suggest 
a distributed source for H 2 CO. Indeed, formaldehyde is a short-lived molecule, with a 
photodissociation rate of 2.8 I O' 4 s' 1 at 1 All from the Sun, according to Huebner and Carpenter 
(1979). The line intensities expected for the 226 GHz transition at half-beam and one beam off the 
nucleus (6 arcsec and 12 arcsec corresponding to about 1 100 km and 2200 km) are about 0.7 and 
0.3 the intensity observed at the centre position, if one assumes that H 2 CO comes directly from the 
nucleus, whereas the observed ratios are 0.8 ± 0.2 and 0.64 ± 0.14. 

The density distribution of H 2 CO in the coma being ill-known, it is difficult to give accurate 
estimates of its abundance. We tentatively compute production rates assuming a daugther Haser 
density distribution, with a parent scalelength of 10 4 km at 0.9 AU from the Sun and scaled 
proportionally to rh 2 , We use a model where both collisional excitation and radiative excitation by 
the Sun are taken into account, and assume a kinetic temperature of 50 K in the coma. The derived 
production rates are given in Table 1: the H 2 CO abundances relative to water are on the order of 
0 5 % for both comets, i. e. much smaller than the 4.5 % found in comet P/Halley from the 
analysis of the 3.6 pm feature in the IKS-Vega spectrum attributed to f^CO (Combes et al., 
1988). For comparison, production rates computed using a parent distribution for H 2 CO are also 
given in Table 1 : as expected, much smaller abundances are derived for the May observing period 
of comet Austin. 

Hydrogen sulfide 

Hydrogen sulfide H 2 S was unambiguously detected in comet Austin on May 23, 24 and 25 1990, 
through its 1 io~ ^oi transition at 169 GHz (Crovisier et al., 1990a), with a S/N ratio of 7 for the 
average spectrum shown in Fig.4. This is the first time that H 2 S is observed in a comet, although 
it was highly suspected by the presence of H 3 S + in the coma of P/Halley (Marconi et al., 1990). 

A preliminary estimation of the H 2 S production rate is given in Table 1 . We assumed rotational 
LTE at 50 K and took a parent molecule density distribution. From the UV absorption spectrum of 
Lee et al. (1987) and the reference solar spectrum of Huebner and Carpenter (1979), we estimate 
the H 2 S lifetime to be about 4000 s at 1 AU. H 2 S is apparently a minor component of cometary 
atmospheres, with an abundance relative to water of 1.5 10 3 . 


Table I : Line intensities and production rates 


Molecule 

Date 

iv** 

|K km s -1 ] 

Q 

ls-'l 

Q(H 2 01 

[s-'l 

Q/OfHjO] 

PIBrorsen-Metcalf ( I989o) 
HCNJ(l-O) 89/09/04-07 

0. 1 8±0.03 

3.0E26 

2.5E29 b 

1.2E-3 

HCN J(3-2) 
H 2 CO 3, 2 -2,, 

89/09/02-03 

89/09/04-07 

0.96±0.2 1 
0.1610.04 

1.0E27 C 

2.5E29 b 

4.0E-3 C 



6.8E26 d 

2.5E29 b 

2.7E-3 d 

Austin (198901) 

HCN J(l-0) 

90/04/02 

< 0.10 

< 1.6E26 

2.0E29 b 

< 8.0E-4 

HCN J(l-0) 

90/05/23 

0.2010.07 

2.4E25 

6.0E28 e 

3.9E-4 

HCN J(3-2) 
H 2 CO 3 I2 -2 n 

90/05/23 

90/04/02 

1.2710.18 
< 0.07 

< 9. 1 E26 c 

2.0E29 b 

< 4.6E-3 C 

H2CO 3 1 2 v 2j j 

90/05/21-25 

0.2510.03 

< 6.7E26 d 
3.6E26 C 

2.0E29 b 

6.0E28 e 

< 3.3E-3 d 
6.0E-3 C 

^2$ 1 1 Q" 1 0 1 

90/05/24 

0.2010.07 

3.3E25 d 

9.0E25 

6.0E28 e 

6.0E28 e 

5.5E-4 d 

1.5E-3 

CH3OH (3,0)*(2,0)A 

90/05/25 

0,2510.05 

4.9E26 

6.0E28 e 

8.2E-3 


a. Observations aimed at the nucleus. 

b. From OH 18 cm observations (BockehSe-Morvan et al M 1990a). 

c. Assuming a daugther density distribution with a parent scalelength of 10 4 km at 0.9 AU from 
the Sun and scaled as rh 2 . 

d. Assuming a parent distribution. 

e. From OH IUE observations (Budzien et al„ 1990). 


Methanol 

Methanol was observed in comet Austin on May 25 1990 through its J(2-l) aK = 0 transitions at 
96 GHz, and its J(3-2) AK = 0 transitions at 145 GHz, using 512x1 MHz filter banks 
(corresponding to a velocity resolution of 2 km s' 1 at 145 GHz). Two J(2-l) transitions are 
marginally present in the 96 GHz spectrum while four J(3-2) lines are detected with signal to noise 
ratios between 2 and 5 in the 145 GHz spectrum (Fig. 5), insuring an unambiguous detection of 
methanol in comet Austin and giving the first evidence of the presence of this molecule in comets 
(Despois et al., 1990). 

Observing simultaneously several rotational lines of CH 3 OH is fascinating since it provides 
clues on its excitation conditions. As a preliminary study we assumed a unique temperature Trot 
for describing the population distribution among the rotational states. From the reference solar 


Austin { 1989c I) 
H2S 110-10! 


Velocity (tm/sJ 

Figure 4 ; The H 2 S l io” *oi spectrum of comet Austin observed on May 24-25 1990 (centre and 
offset positions). 


Austin ( 1 989c 1 
CH3QH 2) 


1 45080 


145100 145120 
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Figure 5 : J(3-2) aK = 0 transitions of CH 3 OH around 145 GHz observed towards comet 
Austin on May 25, 1990. Full circles and crosses indicate the line positions and the relative 
intensities expected for rotational LTE at Trot = 20 K and Trot -70 K, respectively. 
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spectrum of Huebner and Carpenter (1979) and the UV absorption spectrum of Nee et al. (1985), 
we estimate a CH 3 OH photodissociation rate of 1,3 10* 5 s‘ l at 1 AU. The line intensities computed 
with Trot s 70 K and Trot = 20 K are superimposed on the observed spectrum (Fig, 5). The 
observed transitions are much more easily explained by a rotational temperature of 20 K. 
Assuming this value and using the observed intensity of the (3,0)-(2,0) A line at 145,103 GHz 
(Table 1), we derive a methanol production rate of 5 10 26 s' 1 , corresponding to an abundance of 8 
10 3 with respect to water. 

Conclusion 

These new detections confirm that microwave spectroscopy is an efficient tool to study cometary 
parent molecules. The presence of H 2 S, H 2 CO and CH 3 OH, which are well-known constituents 
of the interstellar molecular clouds, supports the hypothesis that comets have retained, at least 
partly, some of the primitive matter of the Solar Nebula without alteration. 

Detailed studies of observations made at offset positions, of the excitation of the new 
molecular species, and of the coma kinematics from the observed line shapes, are in preparation. 
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ABSTRACT 

In 1989 September we conducted VLA searches for both formaldehyde (H 2 C0) 
and cyanoacetylene (HC 3 N) emission from Comet P/Brorsen-Metcalf (1989o). In 
early July, this comet was outgassing strongly at a brighter than average 
preperihelion rate, and our hope was that it would continue to do so 
through perihelion (1989 Sept. 11,9394). By observing the 1-, - 1, 0 
transition of formaldehyde (H 2 C0) at 4,829.659 MHz (6 cm, AJ during the 
phase of maximum outgassing, we expected to collect spectra with a high 
signal -to* Mtim. ratio which would show the temporal behavior of 
formaldehyde relative te other molecules. During this same observing 
period, we also' searched for the J-l-'d, F— 2-1 rotational transition of 
cyanoacetylene at 9098/3321 MHz (3.3 cm. A). 

I. INTRODUCTION 

In order to investigate the composition and outgassing behavior of 
comets, one should observe as many comets as possible in order to get a 
statistically representative sample. We started radio interferometric 
studies of comets several years ago, and after several discouraging 
attempts, we finally succeeded to image OH in Comet Halley in 1985, and in 
Comet Wilson (19861) in 1987 (de Pater, Palmer, and Snyder 1986; Palmer, de 
Pater, and Snyder 1989). We used the same techniques to detect H 2 C0 in 
comets Halley in 1985 and Machholz (1988J) in 1988 (Snyder, Palmer, and de 
Pater, 1989; 1990). We have other data on comets Liller (1988a), Machholz 
(1988j), and Okazaki - Levy - Rudenko (1989r) which are not yec fully reduced. 

In 1989 July we proposed observing Comet Brorsen-Metcalf with the VLA 
because at the time it was outgassing at a brighter than average 
preperihelion rate. It was our hope that it would continue to do so through 
perihelion on 1989 Sept. 11.9, producing cometary H 2 CO strong enough to 
image, rather than to simply detect. By observing H 2 C0 during the phase of 
maximum outgassing, we hoped to collect spectra with a high signal -to -noise 
ratio which would show the temporal variation of H 2 C0 relative to other 
molecules. The VLA referees' reports were generally favorable, except 
possibly for one comment: "H 2 C0 is a "known" comet parent molecule in the 
sense that it has been detected by a number of groups including this one. 

I think all of these detections are somewhat controversial, though, and it 
is important to verify them." Despite this, we were given time to search 
for the lj^ - I^q transition of formaldehyde (1I 2 C0) at 4,829.659 MHz and 
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for the J-1-0, F-2-1 rotational transition of cyanoacetylene at 9098.3321 
MHz (3.3 cm A). 


II. OBSERVATIONS 

We used the VLA to search Comet Brorsen-Metcalf for the 1-. - 1.. 
transition of H 2 C0 on 1989 September 2, 3, and 7, and for the 1-0 of HC 3 N 
on 1989 September 10 and 11. The initial data reduction for formaldehyde 
is farther along than for cyanoacetylene, so we will present our progress 
on the H 2 C0 observations in order to learn how they may correspond with 
other observations of P/Brorsen-Metcalf . The VLA was used with all 27 
antennas in the D (the most compact) configuration and our observational 
parameters were set up just as they had been for our earlier H 2 C0 
observations of comets Halley and Machholz (Snyder et al . 1989; 1990), As 
usual, D, Yeomans' cometary orbital elements were used in C. M. Wade's 
ephemeris program. Two identical spectrometers were used with a total 
bandwidth of 1.56 MHz each. Thirty one (Hanning smoothed) channels were 
recorded from each spectrometer; each channel had a width of 12.208 kHz 
(0.756 km s' 1 ). During the observations, the expected frequency of the 
cometary H 2 C0 line was positioned two channels above the center of the band 
in one spectrometer and two channels below the center in the other. This 
effectively served as a shift in the local -oscillator frequency when the 27 
overlapping channels from the two data sets were combined to form the 
spectra discussed here. The H 2 C0 absorption against NGC 2024 was observed 
at the beginning of the comet observing session to monitor the overall 
performance of the system. The untapered beam was 20" x 20", and the 
tapered beam was 40" x 40". A cometocentric velocity system is used in 
which v - 0 is the predicted velocity of the nucleus. 

III. DISCUSSION 

In our reduction procedure used for Halley and Machholz (Figure 1) , 

Figure 1. a) The 1 XI - 1 10 
transition of H 2 C0 (folded 
about V - -0.76 km s' 1 ) 
observed in Machholz. b) The 
averaged, folded, background 
spectrum made as in la, but 
at 2 random points away from 
the nucleus, c) The 1 1X - 1 10 
transition of H 2 C0 (folded 
about v - -0.76 km s' 1 ) 
observed in Halley, d) The 
averaged, folded background 
spectrum of 5 random points 
away from the nucleus . 

Ordinates : intensity in units 
of 10" 3 Jy (mJy) ; abscissae: 
predicted radial velocity wrt 
the rest frame of the nucleus 
(from Snyder et al. 1990), 


rj c~ 


§- a) Comat Machholz 


b) Average None Spectrum (Machholz) 
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spectra were formed from the data cube - the set of the Images 
of the comet for each velocity channel. By making a cut through the data 
cube at fixed position with respect to the cometary nucleus, a plot of 
intensity versus velocity for that position is obtained. Five spectra, 
taken in the direction of the cometary nucleus and at offsets 10" the 
half-power beamwidth) north, south, east, and west, were obtained and then 
averaged into one spectrum. Next, a version of Gdrard's <1987) SYMCOMET 
routine was applied to the data. In this routine, a search for a 
symmetrical spectral line is conducted by folding the spectrum about an 
axis which is sequentially stepped through the velocity channels. The noise 
in the folded channels is reduced by J2 from its original value; and when 
the axis passes through the center channel of the line, the standard 
deviation computed for all channels (excluding the center channel) will 
reach a minimum. Finally, the same averaging procedure was repeated with 5 
spectra around each of several randomly selected points away from the 
position of the comet, Figure la shows the H 2 C0 emission spectrum and lb 
shows the average of all of the off-comet comparison spectra for Comet 
Machholz. The value of a for all channels except the center channel of the 
line is the smaller of the two error bars located on the left-hand side. 

The larger error bar on the left-hand side represents J2a , the correct 
estimate of standard deviation for the unfolded data at v - -0.76 km s' 1 , 
the center channel and peak of the H 2 C0 emission line. Figures lc and Id 
show similar data for Comet Halley. Both comets had a line of H 2 C0 centered 
at the blueshift v - -0.76 km s' 1 . For both Machholz and Halley, the VLA 
was generally facing the sunlight hemisphere of the nucleus, and this small 
blueshift is consistent with anisotropic outgassing of the nucleus in the 
solar direction. This is illustrated in Figure 2. Also, in the case of 


Figure 2. The geometry of 
comets Halley and Machholz 
at the time of the 
observations in the sun- 
earth coordinate system, 
(from Snyder et al. 1990). 


1986 HjCO Observations of Comet Halley 



1988 H2C0 Observations of Comet Machholz 



Halley, a similar blueshift was found for other cometary species such as 
HCN and H 2 0. 

For Comet Brorsen-Metcalf , we used a new technique for reducing the 
data. Figure 3 shows a hypothetical representation of a data cube which is 
constructed from the interferometric observations of cometary H 2 C0. A 


pixel, represented by a shaded square, is 10 M xl0” in tapered data and 5 M x5" 

Figure 3* The data cube 
represents the set of 
images of the comet for 
each velocity channel as 
a function of a and 6 , 

By making a cut along a 
pixel (which is at fixed 
position wrt the cometary 
nucleus) through the data 
cube, a spectrum of 
intensity versus velocity 
for that position* is 
obtained. Averaging over 
the appropriate number of 
pixels simulates the 
synthesized interferometer 
beam. Here channel 14 is at 
0 km s” 1 cometoeentric 
radial velocity. 

in untapered data. For the data reductions, we used blocks which were 
either 3x3 pixels or 5x5 pixels. Note that 3x3 pixels untapered-lS^lS” , 

5x5 pixels untapered~ 25 , 'x25" , 3x3 pixels tapered-‘30 ,, x30 M , and 5x5 pixels 
tapered=50"x50 !! , so that the size of the untapered synthesized beam is 
bracketed by the untapered 3x3 and 5x5 pixels and the size of the tapered 
synthesized beam is bracketed by the tapered 3x3 and 5x5 pixels. This 
illustrates one of the strengths of using arrays for the radio observations 
of comets: within the field of view, different pixel clusters can be 
sampled to optimize the coupling of the synthesized beam to the gas 
distribution. This is particularly important in the case of comets, because 
small pointing errors caused by an inaccurate ephemeris can be corrected 
after the array observations have been made-a feat which is impossible with 
single dish radio observations. Figure 4 shows a composite spectrum 
displaying the spectra from the 5x5 tapered pixels (solid line), 5x5 
.untapered .pixels (lightly dashed line), 3x3 tapered pixels (dotted line), 
and the 3x3 untapered pixels (heavily dashed line) . There may be a weak 
H 2 CO emission feature at channel 13; since channel 14 is 0 km s' 1 in the 
cometoeentric rest frame, this would correspond to a redshift of 0.756 km 
s" 1 . However, we won't know the statistical significance (if any) of this 
feature until we finish the noise analysis of the data. 

IV. CONCLUSIONS 

We are finishing the noise analysis of our formaldehyde data in order to 
determine its statistical significance and later this year we hope to 
finish the cyanoacetylene data analysis, It will be important to compare 
our centimeter wavelength formaldehyde results with the millimeter 
wavelength formaldehyde results for Comet Brorsen-Metcalf presented by D. 
Bockelee-Morvan at this workshop (Colom et al. 1990) . In order to compute 
the formaldehyde production rate, it will be necessary to model an extended 
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formaldehyde scalelength with a nuclear source and a coma source. 

Both L. E. S . and P. P, wish to acknowledge NASA grant NAGW-1131 to the 
University of Illinois. I. dP. acknowledges support from NSF grant AST- 
8900156 to the University of California, Berkeley. 


Figure 4: Spectra 
obtained by averaging 
over different pixel 
combinations or cuts 
through the data cube. 
Solid line; 5x5 tapered 
pixels. Lightly dashed 
line: 5x5 untapered 
pixels. Dotted line: 

3x3 tapered pixels, ^ 
Heavily dashed line: ~ 

3x3 untapered pixels. 2 
Ordinate: intensity £ 
in mJy; abscissae: 
channel number, 
which runs in the 
direction opposite 
to radial velocity 
wrt the rest frame 
of the nucleus. In 
this rest frame, 
channel 14 * 

0 km s" 1 . 
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Dave Schleicher 

Since the H 2 CO is observed with the 5x5 grid, but not the 3x3 grid, can the 3x3 grid 
location be adjusted to search for the peak in the H 2 CO distribution? 

Lewis Snyder 

Yes, we are working on this, but it will take some time to complete. 
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SILICATES AND VARIABILITY IN 
COMETS OKAZAKI-LEVY-RUDENKO (l989r) AND AUSTIN (1989cl) 

Ray W. Russell and David K. Lynch 
Space Sciences Laboratory, The Aerospace Corp. 


Comet Okazaki-Levy-Rudenko (1989r) was studied over a wide range of wavelengths from 
the UV to the submillimctcr. Using visual magnitude estimates from the IAU Circulars and reported 
submillimctcr fluxes in conjunction with our thermal IR spectroscopy and filter photometry, the nature 
of the grains and the behavior of the comet can be assessed and compared to other comets studied 
under this program. In addition, preliminary results from a collaborative IR study of Comet Austin 
arc presented which address the nature of the dust in that comet. Progress on a model for cometary 
properties and brightness behavior is discussed. 

Introduction 

An extensive IR spectroscopic study of comets was begun with Comet G-Z, continued with 
Halley, and has been extended to include Wilson, Bradftcld, Brorsen-Metealf (BM), Okazaki-Levy- 
Rudenko (OLR), and Austin. Several others did not brighten up sufficiently or were unobservable 
due to weather, The goal of this on-going study is to develop a comprehensive cometary dust model 
that can account for the spectral signatures (or lack thereof) in the thermal IR, and predict, or at least 
explain in an objective fashion, the time variability (or lack thereof) in the cometary brightness curve. 
The current data base of spectroscopy in the IR is not large enough to fully accomplish this, but the 
situation is improving rapidly. This is indicated by the above list of comets, all of which have been 
studied in the last 6 years. This report represents the status of this work as applied to Comets OLR 
(1989r) and Austin (1989cl). 

Observations 
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IR circular variable filter wheel (CVF) spectroscopy and filter photometry of Comet OLR 
were obtained during 5-8 Nov. 1989 UT using the NASA Infrared Telescope Facility (IRTF) circular 
variable filter wheel spectrometer (ATI) and the bolometer filter system (BOLOl) with about 9.5" 
apertures. Magnitudes on 8.8 Nov. 1989 UT were L=6,7, M=4.4 and (10.3]=0.6, 

As part of a collaboration with M. Hanner and H. Campins, Comet Austin CVF and filter 
data were obtained at the 3 m IRTF on 6.8 May 1990 UT with the facility CVF spectrometer (ATI), 
with a resolving power of X/dX~ 50. A 5 mm (9.3 arc see) aperture and 45 arcscc N-S throw were 
used to allow for the extended nature of the comet at a distance from the earth of only .457 au. In 
addition to the usual 7-12.5 pm CVF wedge and moderate resolution filters used for the Halley watch 
observing program, two new filters were installed for this program and left in the ATI CVF dewar. 
They arc a 3.44 pm filler with a full width at half maximum transmission (FWHM) of 0.20 pm and 
an M filter which was measured cold in our lab to have a center wavelength of 4.73 pm and a 
FWHM of 0.61 pm. The fillers were then tested cold in the ATI dewar at the telescope with a 
monochromator and hot source to check for out-of-band spectral leaks. No leaks were found, lending 
credence to the technique of using data from these filters in conjunction with the longer wavelength 
data to derive color temperatures. 

The addition of these filters to ATI allows the characterization of the cometary dust 
continuum through the mid-IR without having to change instruments and calibrate two separate 
systems. In view of the difficulties inherent in (tying to observe close to the sun in the daytime, this 
can be a significant advantage. The ATI sensitivity at short wavelengths is reduced by the low 
quantum efficiency of the doped Si detector, but it was found to be adequate (about 10' 12 W/root Hz 
at 3.44 pm) for the study of a comet bright enough to be studied spectroscopically at 10 pm. The 
3.44 pm filter also permits a rough assessment of the strength of a possible "3.4 pm" emission 
feature above the nominal continuum. 
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The Cornel OLR spectra reported by Lynch and Russell ( B.A.A.S. 22 . 834, 1990) and shown 
in Figure 1 show a strong, structureless 10 pm silicate emission feature extending about 20% above 
the continuum. It was present on all three days, but was about a factor of 1.5 brighter on 7 Nov. At 
the same time, the visual magnitude estimates shown in Figure 2 (IAUC 4892, 4904, and 4908) 
indicated a .5 -.8 magnitude brightening as well. This comet seems to be more variable than, for 
example, Comet Wilson (19861), and of the order of or a little less variable than Comet p/Hallcy. 

The variability of OLR was reported to be about a factor of 3 during three days of observations in 
the submillimctcr (IAUC 4914) in the 18-23 Nov. time frame. 

The L-M-7.8 micron color temperature was 380 K, approximately 10% above the equilibrium 
temperature of 345 K. Figure 3 shows the spectroscopy data along with a normalized black body 
curve at the equilibrium and elevated temperatures for comparison. The grain temperature elevation 
compared to a black body at the same heliocentric distance indicates that the grains were inefficient 
radiators, and therefore small compared to the wavelength of the observation. This, in concert with 
the existence of silicate emission which is indicative of a wavelength dependent cmissivity, suggests 
that the dust grains were optically thin at these wavelengths. We can compare the estimates of dust 
mass derived from the IR and submillimctcr data, and find them to be discrepant by about three 
orders of magnitude. The source of the disagreement is being investigated. It may be model 
dependent or may simply reflect the variation in the dust production rate during this apparition of 
Comet OLR. 

Several spectra of Comet Austin were obtained in different modes in an effort to 
defeat the effects of thermal changes in the optics due to differential heating by (he sun shining on 
the telescope structure. The two best, complete spectra showed no 10 micron silicate emission 
exceeding 20% of the continuum and arc consistent with the absence of silicate emission. However, 
to accurately determine the presence or absence of a silicate feature with low contrast, the level of 
the continuum must be precisely determined. We chose to use the measurements of the continuum at 
3.44 and 4.73 microns (where blackbody flux distributions for typical cometary temperatures exhibit 
exponential dependence on temperature and thus even fairly large measurement errors translate into 
relatively small temperature uncertainties) and at the edges of the 8-13 micron window where the 
silicate has much lower cmissivity than in the 9-1 1 micron region. 

The IR magnitudes of Comet Austin obtained with the fixed fillers were 3.44 microns, 
6.4+/-0.35; M(4.73 microns), 4.7 +/-0.1; N(1Q.1 microns), 0.18 +/-0.05; 1 0.3 microns, 0.13+/-0.05. 
These data arc consistent with a slightly elevated color temperature and the presence of a 3.4 pm 
emission feature (Mumma ct al. IAUC 5012). 

Fitting a 300 K greybody to the ten micron spectrum only results in a curve which misses 
the 3.44 and 4.73 micron values completely. Figure 4 shows a 335 K greybody normalized to the 
4.73 micron point and just below the ends of the 10 pm spectrum, allowing for the fact that typical 
interstellar silicate cmissivity curves arc not zero at either 8 or 12.7 microns. The cometary spectral 
flux is normalized to the level of the 10.3 filter data, and is about 8 % lower than the other good 
quality, complete spectrum obtained the same day. If the comet spectrum in Figure 4 were raised to 
the level of the other narrowband spectroscopic data, the silicate feature would be correspondingly 
larger relative to the continuum. On the other hand, if one were to increase the comet flux at 4.73 
microns by one sigma (10%) and then fit the 335 K greybody to that value, no silicate feature would 
be present. 

We believe the intermediate normalization based on the 10.3 pm filter photometry shown here 
to be the most probable. Furthermore, considering the shape of the smoothest comet spectrum we 
obtained (Figure 4), the most likely choice is that there was a small silicate feature (about 10%) and 
that there may even be a break in the spectrum just past 1 1 microns which is seen in both complete 
spectra. This break was seen in Halley (Campins and Ryan 1989, Ap.J. . 341 . 1059) and Bradficld 
(Russell ct al. 1990, irt preparation) and identified by Campins and Ryan as due to a small admixture 
of crystalline olivine grains in the cometary dust. 

It should be pointed out that a small silicate feature like that possibly present here would not 
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be detected at all by broadband filter photometry. This emphasizes the value or the two new filters 
in ATI in conjunction with the CVF wedges for this type of program. The use of one dewar also 
means that all of the data can be taken through the same aperture with the same optics, thus 
removing another potential ambiguity from the data, 

These comets fit well in a picture where the properties of the dust and the behavior of the 
brightness during the course of the apparition determine the type of comet under study. This 
approach to studying comets is in contrast to one in which models for different types of comets are 
based on their orbits (short or tong period, or first time/one-visit comets). The comets which are not 
dusty seem (based on our observations of Brorsen-Mctcalf, for example; Lynch ct al„ this volume), to 
have a carbon rich grain population which emits greybody spectra with little or no silicate feature in 
the 10 micron region. They seem to lack rite dramatic variability seen in Comet O-L-R (1989r) or 
p/Hallcy, both of which had substantial silicate emission features. Austin may be an intermediate 
case between these two classes. The visual magnitude estimates from the IAU Circulars for Austin 
indicate, with the exception of one period of steady brightening, an essentially non-periodic variability 
of about a factor of 2 over the span of a few days. 

In a comet with a lot of volatile material on the surface, a more or less steady outflow is 
expected. In contrast, a comet with a crusty surface could behave cither in a minimalist fashion (i.c., 
very little sublimation at all, and thus fairly steady in brightness but at a low level) or in the very 
sporadic, outburst mode observed for p/Hallcy. This tough, crusty surface may well be mostly 
hydrocarbons which form a sticky layer, essentially black, which shields the inner parts of the 
comctary nucleus and emits the grey-body shape reported for B-M and closely approximated by the 
data shown here for Austin (silicate enhancement only about 10% above the greybody continuum). 

Additional spectra of a statistically significant sample will be required to test this model and 
flesh out the classification characteristics. We have such an ongoing program using our new 
broadband array spectrograph to span the 3-13.5 pm region simultaneously, thus avoiding problems of 
tracking and spatial structure which have plagued our CVF measurements. The first of these runs is 
scheduled for the next quarter on the 1RTF, 

Sum mary 

Thermal IR spectroscopy and mid-IR filter photometry of Comet OLR have been reported 
which show a strong silicate feature and significant temperature enhancement over a black body in 
equilibrium. Even though OLR was a "new" comet, it has many similarities in its spectral shape and 
temporal behavior to p/Hallcy. Similar observations of Comet Austin are suggestive of a small 
(about 10%) silicate feature superposed on a 335 K greybody continuum, and less temporal variation 
than Halley or OLR, but more than Wilson. The temperature derived for the dust in Austin is quite 
close to the blackbody equilibrium value of 314 K, and may be indicative of grains whose size is in 
the 1-10 micron range. The data arc also consistent with a carbon-rich grain population which would 
be expected to emit essentially greybody emission. As noted above, the comet exhibited variability 
in the visual and IR which was intermediate in nature bctwci»r> that of Halley and Brorsen-Metcalf, 
silicate- and carbon-rich comets respectively. 

This program clearly points out the need for IR spcctrographic data on a statistical sample of 
comets to take our understanding of comctary dust grains and nuclear sublimation forward. 
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DETECTION OF THE 3.4 MICRON EMISSION FEATURE IN COMETS fl< 
P/BRORSEN-METCALF AND OKAZAKI-LEVY-RUDENKO (1989r) AND 0 
AN OBSERVATIONAL SUMMARY 


T.Y. Brooke, A.T. Tokunaga 
Institute for Astronomy, University of Hawaii 
2680 Woodlawn Dr., Honolulu, HI 96822 
and 

R.F. Knacke 

Space Sciences Lab ES-63 

NASA Marshall Space Flight Center, Huntsville, AL 35812 

We are continuing a survey of the 3 fim spectral region of comets, including the broad 
3.4 /im emission feature due to C-H bonds in organic molecules. One goal is to determine 
how the strength of the emission band varies with heliocentric distance r in comets. This 
depends on both the production rate of the organic molecules and the emission mechanism, 
neither of which is well determined at present. 

In 1989, the 3.4 /tm feature was detected in two comets: the periodic comet P/Brorsen- 
Metcalf and new comet Okazaki-Levy-Rudenko (1989r), hereafter OLR. Spectra of the two 
comets obtained at the NASA Infrared Telescope Facility on Mauna Kea are shown in Fig. 
1. The feature was clearly evident in P/Brorsen-Metcalf down to r = 0.48 AU, the smallest 
heliocentric distance at which these cometary organics have been detected. 

The 3.4 /im features in P/Brorsen-Metcalf and OLR were spectrally similar to those 
observed in other comets. But the feature- to-continuum ratios in these two comets were 
higher than those expected from the trend seen in other comets to date. Fig. 2 shows 
the feature-to-continuum ratios at 3.36 /im as a function of heliocentric distance for all 
published 3 /im spectra of comets, including 3<r upper limits for spectra in which no feature 
was evident. Complete references for the data are given in Brooke, Tokunaga, and Knacke 
(1990). 

The spectral contrast of the 3.4 /im feature in comet P/Halley was generally lower 
at smaller heliocentric distances for r in the range 1—2 AU (Fig. 2). Comets Wilson 
(1987 VII) and Bradfield (1987 XXIX) appear to follow the Halley trend, though the data 
are limited. The feature-to-continuum ratios in P/Brorsen-Metcalf and OLR were higher 
than any reasonable extrapolation of this trend. This shows that the flux of the 3.4 /im 
emission is in general not coupled to the flux of the underlying continuum. 

Separate comparisons of the dust continuum and 3.4 /im band fluxes of comets suggest 
a plausible explanation for the higher feature-to-continuum ratios in P/Brorsen-Metcalf 
and OLR. Fig. 3a shows the dust continuum flux at 3.36 /im ratioed to the water produc- 
tion rate at the time of the observation for each of the observations represented in Fig. 2. 
Water production rates are used as a means of normalizing fluxes of comets with widely 
different levels of activity. Measured continuum fluxes were reduced to a common view- 
ing geometry assuming conservative radial outflow, i.e., flux proportional to the projected 
radius of the aperture. Details of the reduction and references for the water production 
rates are given in Brooke, Tokunaga, and Knacke (1990). 


Fig. 3a clearly shows that comets P/Brorsen- Metcalf and OLR, along with P/Encke, 
had lower 3.36 pm continuum fluxes relative to their water production rates than the 
values that would be expected for Halley, Wilson, or Bradfield at comparable heliocentric 
distances. 

Integrated 3.4 pm band fluxes were estimated by subtracting continuum fluxes from 
the spectra. Fig. 3b shows the band fluxes and upper limits ratioed to the water production 
rates for the same sample of comets shown in Fig. 3a. The same scaling to a common 
viewing geometry was used. In marked contrast to their continuua, the normalized 3.4 pm 
band fluxes of P/Brorsen- Met calf and OLR are consistent with the values of other comets. 

There is an approximate correlation between the 3.4 pm band flux and the water 
production rate at a given heliocentric distance for all of the comets in this sample (Fig. 
3b). The high feature-to-continuum ratios in comets P/Brorsen-Metcalf and OLR are 
thus most likely due to their weaker 3 pm dust continuua relative to water. These two 
comets also had weak dust continuua in the visible (Schleicher et al. 1989) and ultraviolet 
(Roettger et al. 1990) compared to their gas production rates; in fact, the two are among 
the most “dust-poor” comets. The present results suggest that comets with low dust-to-gas 
ratios will exhibit higher feature-to-continuum ratios at 3.4 pm. 

The heliocentric distance dependence of the normalized 3.4 pm band flux (Fig. 3b) 
can be used to constrain the emission mechanism if the production rate of organics relative 
to water is assumed to be the same in all comets. The present data are consistent with 
fluorescence by molecules excited by solar radiation, in which case the normalized band flux 
goes as r -2 . But simple models of thermal emission by a population of small, hot grains 
can also match the data (Brooke, Tokunaga, and Knacke 1990). Possibly both processes 
contribute. If the emission is from grains, it is interesting that the 3.4 pm flux from these 
grains is better coupled to the water production than the dust continuum, suggestive of 
two separate grain populations. 

The observations to date are consistent with cometary organics being present in all 
comets at comparable abundances relative to water. Loss of contrast in the 3.4 pm feature 
as comets get closer to the sun is then easily explained by dilution by thermal emission 
from the continuum grains, whose flux rises more steeply with heliocentric distance than 
that of the organics. The lack of a definite 3.4 pm feature in some comet spectra for r < 1 
AU could be due to insufficient signal-to-noise (Fig. 2). Upper limits to the band flux 
(Fig. 3b) suggest that 3.4 pm emission could have been present in these comets at levels 
consistent with those of other comets. 
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Fig. 1-Spectra of comets P/Brorsen-Metcalf and Okazaki- Levy- Rudenko taken with the 
Cooled Grating Array Spectrometer (A/AA = 150 — 200) at the NASA Infrared Telescope 
Facility on Mauna Kea. The spectra are composites of data taken on successive days 
with different grating positions. All spectra were obtained with a 2.7 arcsec diameter 
aperture centered on the infrared peak. The uncertainty in the absolute flux calibration is 
approximately 30 %. From Brooke, Tokunaga, and Knacke (1990). 
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j'ig. 2-Feature-to-continuum ratios as a function of heliocentric distance in comets. Down- 
ward arrows indicate 3 <r upper limits. 
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Fig. 3-(a) Dust continuum fluxes at 3.36 fj,m normalized by the water production rates. 
F* ont is the dust continuum flux reduced to a geocentric distance of 1 AU and a projected 
aperture radius of 980 km. (b) Integrated 3.4 um band fluxes normalized by the water 
production rates. If A is the integrated band flux reduced to a geocentric distance of 1 AU 
and a projected aperture radius of 980 km. Symbols are the same as Fig. 2. 
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COMET BRORSEN-METCALF IN THE 3.5 MICRON REGION 

A. D. Storrs (U. Texas), M. J. Mumma, and S. M. Hoban (GSFC) 

We observed comet Bror sen- Metcalf on UT August 21-24 1989, using 
the CRSP spectrometer and the 1.3m telescope at Kitt Peak National 
Observatory. The cometary continuum was detected on all nights. 

The data are very noisy, due to the short observation window and 
the untried nature of the instrument. Low resolution (0.0227 micron) 
spectra show the 3 .4 micron C-H stretch feature having a contrast of 
at most a factor of two to the neighboring continuum. High resolution 
(0.00331 micron) spectra between 3 .4 and 3.6 microns some 1 sigma 
features that might be attributed to the nu5 band of formaldehyde 
(H2C0). Similar spectra of the region between 3.2 and 3.4 microns 
show one 3-sigma line at 3 34 microns, which is as yet unidentified. 
Although the cometary spectra were more spatially extended than 
the spectra of standard stars, no extension of the line emission 
beyond the continuum was observed. 
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CARBON ISOTOPES IN COMETS 

Peter A. Wehinger 

Physics-Astronomy Department, Arizona State University 




This paper is a progress report based on high resolution spectra of selected bright 
comets with the goal of determining the carbon isotope abundance ratio, 12 C/ 13 C. The 
12 C/ ,3 C ratio has been determined for various Solar System objects (in the atmospheres 
of the giant planets, meteorites, the Earth, and the solar photosphere), where 12 C/ 13 C = 
89/1. In the interstellar medium, optical and radio observations give a range of 12 C/ 1S C 
= 43-67/1 depending on the observing techniques used and the specific interstellar cloud 
observed. For comet P/Halley, Wyckoff et al. (1989 Ap. /., 339, 488) found 12 C/ 1S C = 
65/l±8. Previous determinations of the 12 C/ 13 C ratio in comets by other investigators, 
involving lower resolution spectra, dealt with the C2(0,l) Swan system with band heads 
at 4737 A for 12 C 2 ( 0 , 1 ) and at 4745 A for 13 C 12 C(0,1). However, there is an NH 2 feature 
with four lines which coincide with the 13 C 12 C band head, thus making a determination 
of the 12 C/ 13 C ratio very uncertain based on the C 2 ( 0 ,l) band structure. 



Here we present echelle spectra of the CN(0,0) violet system in three additional comets: 
P/Brorsen- Metcalf (observed during August 1989), C/Okazaki-Levy- Rudenko (C/1989r) 
(December 1989), and C/Austin (C/1989cl) (May 1990). P/Brorsen-Metcalf has a period 
of 70 y (prograde) compared with P/Halley which is 76 y (retrograde). The similar periods 
made P/Brorsen-Metcalf of special interest for comparison with P/Halley. The other two 
comets were simply comets of opportunity. While not shown here, we have also acquired 
echelle spectra of C/Wilson, All four comets were observed when they were at or close to 
maximum brightness, ~5-6 mag for the integrated magnitude of the coma. Of the four 
comets discussed here, only P/Halley (~4.0 mag) has relatively high dust-to-gas production 
rates. Little or no scattered solar continuum was detectable in comets P/Brorsen-Metcalf, 
C/Okazaki-Levy- Rudenko, and C/Austin. 

All of the spectra were acquired with the Mount Stromlo Observatory's 1.9-m reflector 
near Canberra, Australia, except for P/Brorsen-Metcalf which was observed with the Kitt 
Peak National Observatory 1-m coude-feed telescope. Our observations consist of coude- 
echelle spectrograms of the CN(0,0) B 2 E + - X 2 E + violet system. Our immediate objective 
is to search for 13 C 14 N(0,0) lines in the R-branch (AJ=+1) and to compare these spectra 
with those of P/Halley at approximately the same spectral resolution. 

The instrumental details for the Mount Stromlo observations are: 1.9-m f/32 coude 
focus, spectral range: AA = 3862-3876 A, spectral resolution: A/AA = 80-100,000, AA 
= 40*50 mA(FWHM), typically ~2.5 pixels, with slit dimensions: 0.8 x 22 arc sec. The 
600 line/mm echelle grating was used with a cross disperser and an 81.3-cm focal length 
Schmidt camera. The two-dimensional detector is a photon counting array (PCA) which 
employs two Litton micro-channel-plate intensifies, an f/1 Nikon transfer lens, and a 
Fairchild 244 x 750 CCD. The PCA has inherently low noise characteristics and is read 
out continuously at tv-frame rates, thus making it an ideal system for detecting weak 
spectral lines at low contrast in extended objects such as comets. 
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The spectra of P/Brorsen-Metcalf were acquired with the 1-m f/30 coude feed tele- 
scope at KPNO. The slit dimensions were: 1.5 x 134 arc sec, while the spectral range was 
the same as the Mount Stromlo spectra. The detector was a Texas Instruments 800 x 800 
CCD with quantum efficiency of about 20% at 3870 A. The No. 5 Schmidt camera was 
used with an echelle grating and an order-sorting filter centered on the 3860-3880 A region, 
giving a spectral resolution of 100,000. The quality of the data was limited by the readout 
noise (12 e~ per pixel per read) which was comparable in intensity with the 13 C 14 N lines. 
Integration times were limited to 15 min because of frequent and numerous cosmic ray 
hits. 

In Figure 1 we present a moderate resolution spectrum of the CN(0,0) violet system 
in P/Halley showing both the P- and R-branches covering the spectral range 3850-3890 A 
obtained with the Mount Stromlo 1.9-m reflector, 81.5-cm camera, and a plane grating, 
with a spectral line width of ~200 mA (FWHM). Figure 2 shows the CN(0,0) R-branch. 
in P/Halley (1986 April 4). Figure 3 is a composite spectrum of the CN(0,0) R-branch 
in P/Brorsen-Metcalf coadded over five nights (1989 August 15-19). Figure 4 shows the 
CN(0,0) R-branch and part of the P-branch in C/Austin (1990 May 19 and 20), Detail of 
the weak 13 C 14 N R4, R5, R6 lines are shown for P/Halley (Figure 5), P/Brorsen-Metcalf 
(Figure 6), C/Okazaki-Levy-Rudenko (Figure 7), and C/Austin (Figure 8), 

In Figures 6, 7, and 8 possible detection of 13 C 14 N lines is indicated by the vertical 
arrows. However, we must stress that the identification of these weak features, as well as 
any resulting 12 C/ 13 C abundance ratios, depend on comparisons with computed synthetic 
spectra, which simulate the Doppler-shifted solar flux for each integration over the course 
of two nights to two weeks. In a separate paper presented at this workshop, by Marvin 
Kleine (Arizona State University), full fluorescence calculations for the CN violet system 
are described. Comparison of our observations and Kleine’s synthetic spectra are currently 
in progress. 

These observations, their reduction, and analysis represent the collaborative effort 
of Bruce Peterson (Mount Stromlo Observatory), Susan Wyckoff, Marvin Kleine, Maria 
Womack (all at Arizona State University), and Eric Lindholm (Boston University). We 
are grateful to the Directors of Mount Stromlo Observatory (Alex Rodgers) and Kitt Peak 
National Observatory (Dave De Young) for telescope time. This program is supported, 
in part, by NASA Grant NAGW-547 and the NSF U.S.- Australian Cooperative Science 
Program. I am most grateful to Bruce and Vickie Peterson for their kind hospitality on 
numerous my visits to Canberra. 
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Figure 1. P/Halley CN(0,0) violet system P- and Rebranch spectrum obtained with 
the 1.9-m Mount Stromlo reflector and a plane grating coude where A/AA as 25,000 (1986 
March 26). 



Figure 2. P/Halley CN(0,0) violet system Rebranch spectrum obtained with the 1.9-m 
Mount Stromlo reflector and the coude echelle spectrograph where A/AA = 100,000 (1986 
April 4). 
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Figure 3. P/Brorsen-Metcalf CN(0,0) violet system Rebranch spectrum obtained with 
the 1-m Coude Feed telescope at Kitt Peak National Observatory where A/AA = 100,000 
(1989 Aug 15-19). 
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Figure 4. C/ Austin CN(0,0) violet system R-branch spectrum obtained at Mount 
Stromlo (1990 May 19 and 20). 
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Dave Schleicher * - v 

1) What is the aperture size that was collapsed into the single spectrum? 

2) Is the apparent continuum actually a continuum or just the wings of the strong CN lines 
due to the instrument profile? 

Peter Wehinger 

The typical slit dimensions were ~ 20 x 1 arc sec. 

In our spectra of comet Austin and Comet Okazaki-Levy- Rudenko scattered solar continuum 
has not yet been removed. However, with the low dust production in these comets, corrections 
for solar continuum is small, compared to Comet Halley. 
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N 2 /NH$ Abundance Ratios in Star-Forming Regions 

A topic of importance in both cometary and interstellar astronomy is the N2/NH3 
abundance ratio in star-forming regions. A key test for solar system formation models is 
the N2/NH3 abundance ratio in comets, which are thought to be comprised of material 
from the protosolar nebula, and in molecular clouds, from which stars are thought to 
form. The N2/NH3 abundance ratio has been measured for comet Halley and found to 
be ~ 0.1 (Wyckoff and Theobald 1989 ). The N2/NH3 abundance ratio has never been 
measured in the ISM, but is estimated from theoretical models (i.e. Herbst and Leung 
1989 ) to be in the range of 100-200. If the theoretical models are truly representative of 
dense molecular clouds then the difference in the abundance ratios of N2/NH3 between 
comets and the ISM would imply that significant chemical processing must have occurred 
in the protosolar nebula that somehow altered the amount of N 2 relative to NH3. We have 
undertaken a systematic survey of star forming regions to establish N2 abundances and 
N2/NH3 ratios for comparison with cometary data. 

Abundances of NH3 in dense molecular clouds are well-established (e.g. Mauers- 
berger et at. 1986 ). Molecular nitrogen, however, is not observable with radio techniques 
due to lack of a permanent dipole moment. Protonated N 2 , (or N 2 H + ), is observable and 
therefore N 2 H + abundances can be measured and N 2 abundances deduced if we assume: 

1) All of the N 2 H+ is produced from the reaction 

n 2 + # 3 + -» n 2 h + + h 2 

2) All of the destruction of N 2 H + is due to electron recombination 

N 2 H+ + e" -4 N 2 + H 

3 ) Observed NH3 and N 2 H + line emission arise from the same volume. 

4 ) [N 2 ] ~ |N 2 H + ](k 2 /ki) where ki = rate constant for electron destruction of N 2 H + 
and k 2 = rate constant for formation of N 2 H + . 


Observations of the N 2 H+ J = 1 -4 0 (93 GHz), J = 3 2 (279 GHz) and N 15 NH+ 

J = 1 — ♦ 0 (91 GHz) transitions were performed with the NRAO 12-m telescope on Kitt 
Peak. The N 2 H + abundances were calculated using observed line temperatures of both 
the J = 1 -4 0 and J = 3 — » 2 transitions of the main isotope together with radiative 
transfer equations and a Large Velocity Gradient model. The spectra of the less abundant 
isotope transition were used to help constrain the optical depths of the clouds. 

The dense (n(H 2 ) ~ 10 4 to 10 6 cm** 3 ) molecular clouds observed are listed in Table 
1. Analysis of the observational data to determine the N2/NH3 abundance ratios in these 
sources is in progress. 


TABLE 1 

DENSE MOLECULAR CLOUDS WITH N 2 H+ 


Cold (£ 10K) 

Warm (£ 75K) 

L134N 

W3 

TMC-1 

NGC-7538 

ORION-3N1E 

W51 

p Oph 

W49 

IC443G 

DR21(OH) 

SGRB2 

ORION-S 

ORION-KL 


ill 


Two Clouds At Orion-KL 


A double velocity component was detected in N2H+ data obtained toward the Orion- 
KL (KL) region, which is located in the Orion Molecular Cloud, about 500 pc from the 
Sun. Orion-KL is considered to be a prototypical region of star formation, indeed it is the 
brightest infrared region of Orion (see Kleinmann and Low 1967). 

Evidence for a shift in radial velocity ranging from +8 to +10 km/s across the KL 
region has been observed (Ho and Barrett 1978, Bastien et al. 1981 and Plambeck et al. 
1985). Two possible explanations have been proposed: 

1) The velocity shift is due to rotation about- a central outflow source where the 
velocity differences are due to the redshift and blueshift of the rotating gas (Plambeck et 
al 1985). 

2) Two clouds with different radial velocities lie along the line of sight to Orion-KL. 
One to the north with V isr ~ 9.8 km/s and another to the south V LSR ~ 7.5 km/s (Ho 
and Barrett 1978, Bastien et al 1981). 

Results 

A spectrum of the N2H+ J = 1 — ♦ 0 transition is presented in Figure 1. The hyperfine 
structure of the AF\ = +1 (93.173809 GHz), -1 (93.176310 GHz) and 0 (93.171947 GHz) 
(Sastry et al 1981) transitions is seen. Two velocity components are apparent in the 
spectrum; one to the northeast (NE) of KL with an average velocity of 9.5+0.3 km/s and 
a weaker component to the southwest (SW) with an average velocity of 7. 2+0.3 km/s. 

In order to study the spatial extent of the clouds giving rise to the two velocity 
components, N2H + was mapped around Orion-KL with a 70 n beam. A schematic diagram 
of the spectra obtained are presented in Figure 2. Measured velocities, V lSRi from the 
individual spectra are plotted on a grid (a, 6) and two velocities are plotted for a given 
position if both velocity components were detected in a single spectrum. We have proposed 
boundaries of the two clouds with curved lines. Several compact infrared sources, including 
at least one protostar, have been found in the region and are plotted on top of our N2H + 
data. The infrared sources lie at the projected overlap of the two clouds. Therefore we 
suggest that this may be evidence for a cloud-cloud collision triggering star formation. 

Summary 

A survey of N2H + in 17 star-forming regions in the galaxy has been completed. N2 
abundances will be derived and N2/NH3 abundance ratios established in the molecular 
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clouds. The N2/NH3 abundance ratios will be compared with the same ratio in comets to 
assess the amount of chemical processing which may have occurred during the formation 
of the solar system. 

As part of this survey, a double velocity structure was discovered toward Orion-KL, 
a region of substantial star formation. We propose a two cloud model for the Orion-KL 
nebula, where two velocity components are unambiguously detected - one to the NE of 
Orion-KL at Vlsr ~ (9.5dt0.3)km/s and the other to the SW of Orion-KL at Vlsr ~ 
(7.2±0.3)km/s. Compact infrared sources (including at least one protostar) lie at the 
projected interface of the two clouds and suggest that a cloud-cloud collision may have 
triggered the observed star forming activities. 
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Vls* (km/s) 


Figure 1. The N 2 H+ J = 1 — ► 0 emission toward Orion-KL. Two 
velocities are seen for each AFj hyperfine component of N 2 H+ and are 
separated by ~ 2.5 km/s. 



Figure 2. Simple model for the two cloud components of Orion-KL 
with the measured LSR radial velocities of two clouds as shown. Locations 
of compact IR sources (Downes tt aL 1981) are plotted. The position of 
Orion-KL is noted with an *. Note that the IR sources are found within 
the projected overlap region of the two clouds, which we suggest i 5 evidence 
for a cloud -cloud collision triggering star formation. 
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Ray Russell 

Are there any other lines in this part of the spectrum that could act for the double set of lines? 
Maria Womack 

No. We were able to resolve the hyperfine lines of the N 2 H + (J = 1 — » 0) transition in one 
of our sources, L134N. We used this spectrum to identify the 1 -» 0 of N 2 H + toward Orion - 
KL. 


Walter Huebner 

Since you base your N 2 abundance on N 2 H + , it is important to understand the production 
of NaH + . How does one produce an ion in the interior of a dense interstellar cloud? If this 
is not understood, then the N 3 abundance in the center of the cloud could be quite different 
since the observed N 2 H + may be only on the surface of the cloud. 

Maria Womack 

N 3 H + is thought to be produced by ion-molecule reactions, which typically proceed at low 
temperatures. The strength of both the J = 1 -* 0 and J = 3 — ► 2 lines toward hot, dense 
regions was unexpected and we will work to identify the production mechanism of N 2 H + in 
these hot, dense clouds. It has been suggested that grain chemistry may be important. 


Timothy Brooke 

Your technique will only give you N3/NH3 in the gas phase in molecular clouds. For sompa*- 
ison to comets, don’t you also need a handle on the form nitrogen takes in grains? 

Maria Womack 

Yes, we will also consider the grains as a carrier of cometary nitrogen. 
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Introduction 

Measurement of the carbon isotope ratios in comets provides direct information about 
the chemical state of the interstellar medium (ISM) 4.5 billion years ago. All solar system 
objects (e.g. sun, planet atmospheres, meteorites and lunar rocks) for which the carbon 
isotope abundance ratio has been measured indicate 12 C/ 13 C = 89±2. This value is 
different from the ratio measured in the present-day local interstellar medium (ISM) where 
the carbon isotope ratio has been, measured to range from 12 C/ 13 C ~ 40 to 70. Models for 
the chemical evolution of the ISM (e.g. Audouze 1985) attribute the difference between 
the solar system and ISM carbon isotope ratios as due to a gradual 13 C enrichment of the 
ISM largely from red giant stars and novae over the past 4.5 billion years. However, a 
measurement of the carbon isotope ratio in comet Halley (Wyckoff et. al. 1989) indicated 
a value significantly different from the solar system value, namely, 12 C/ 13 C = 65± 9. If 
the comet formed coevally with the solar system, it is difficult to explain this abundance 
ratio. 

Because of the importance of the carbon isotope ratio in understanding the formation 
of the solar system and the long-term chemical evolution of the ISM, we have initiated 
a program to observe the CN B-X (0-0) and C 2 Swan (0-0) bands in a sample of bright 
comets. An important factor in determining the carbon isotope ratio is the fluorescence 
model used in the analysis. Here we present high resolution CN spectra of several comets, 
and describe the model presently in the process of development which will be used to de- 
termine the carbon isotope abundance ratios from observations of the R-branch rotational 
lines of the CN B-X (0-0) band. 

Swings (1941) was first to demonstrate that the CN band in comet spectra is excited 
by solar fluorescence. In fact he showed that the observed spectrum is extremely sensi- 
tive to the heliocentric velocity of the comet, because Doppler shifts of only a few km/s 
can cause large changes in the fluorescing solar flux due to the complex absorption line 
structure in the solar spectrum in the 3860 - 3880 A region. In the model described here 
we assume fluorescence equilibrium (established on timescales of minutes in comet comae) 
and compute the relative populations of the rotational energy levels in the three electronic 
states: X, B and A. In exciting the upper electronic states we account for the relative 
Doppler shifts of the comet and the sun, the so-called Swings effect. The model includes 
848 rotational levels of the molecule distributed between the three electronic states. The 
g-factors or fluorescence efficiencies are calculated from 

9line ~ 

for each observed line for each isotope, 12 C 14 N, 13 C 14 N, 12 C 15 N. The isotope abundance 
ratio is then computed from the observed line ratios, I/I,-, where I is the intensity of the 
normal isotope, and Ij is the line intensity of the less abundant isotope, from 

N/Ni = (rg.Wig)- 


11 s 


Fluorescence Calculations 


The evaluation of the fluorescence efficiency for each transition, requires the solution 
of a set of linear equations. These linear equations describe the equilibrium occupancy of 
the molecular energy levels. The set of equations are 


and 


where 


d N * m ’~ l 

= ^ “ BjiXjUy) — ^ ' (AjiXj — B, } X,U „) — 0.0 

imj+\ i=l 



1.0 


N = number of energy levels, 

A,j = probability of spontaneous emission, 
Bji = probability of spontaneous absorption, 
U„ = energy density of fluorescing radiation, 
z, = fractional occupancy of the energy level, 


Xj = fractional occupancy of energy level under consideration. 

This results in a N x N matrix that must be inverted to deterime the *<’*. The g- factors 
or fluorescence efficiencies are calculated from 


9linr — *i Aij 


for each line. 

To accurately model the Swings effect requires the use of a high resolution solar atlas. 
The atlas due to Kurucz et al. (1984) was utilized. The variability of the flux responsible 
for the pumping of the ground vibrational state of the electronic state B is shown in figure 
1. The values for the molecular constants of CN are taken from Ito et al. (1988) and 
Koltar et al. (1980). The oscillator strengths were taken from the work of Knowles et al. 
(1988). 

The input parameters to the CN spectral synthesis model are the heliocentric distance, 
heliocentric velocity, and parameters of the molecular isotope under evaluation. 

Summary 

The CN B-X (0-0) R branch fluorescence spectra for three comets (Austin, Bennett, 
Halley) were calculated and compared to the observe spectra (see figures 2,3,4). The 
synthesized spectrum in each case has been shifted to allow for the comparision of the 
individual line intensities. The Swings effect is clearly visible in each figure. In addition to 
the R(0-0) band in the CN spectrum, the P(l-l) band also contributes lines in the spectral 
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region of observation. This branch was calculated and compared to P/Halley (fig. 5). 
The overall agreement between observations and the CN spectral model is quite good. An 
exception to this statement would be the high J rotation lines intensities within the R(0-0) 
branch of P/Halley. Additional analysis is necessary to determine if this difference is due 
to mechanisms in addition to fluorescence (e.g. collisions). 

Further sensitivity analysis of the model is planed to place the comparison with 
observations on a quantitative level. This model will be applied to cometary spectra 
(Wehinger, these proceedings) to calculate 12 C/ 13 C abundance ratios from the analysis of 
the B-X (0-0) band of the CN radical. 

REFERENCES 

Audouze, J. 1985, Production and Distribution of C, N, 0 Elements , ed, I.J.Danziger, 
F.Matteucci, and K.Kjar (ESO Conf. Proc. 21) (Garching:ESO) p.373. 

Ito,H., et al.,J.Mol.Spectrosc, 127, 1988 

Knowles, P., et al.,J.Chem.Phys.,89,1988 

Koltar,A., et al.,J.Mol.Spectrosc,80,1980 

Kurucz,R., et al.,1984,ATattona/ Solar Observatory Atlas No. 1 

S wings, P., 1941, Lick Obs.Bull.,19,No.408,131 

Wyckoff,S., et al., 1989, Ap.J., 339, 488 


SOLAR FLUX (R. Kurucz et al.) 



118 



Relative Intensity Relative Intensity 


Cornel Austin 

Comparision Between Observation at 4 Model 

Observation I 

— Model' , , \ 

R * 1.06 AU, V « 33.S km/s ! , | $ 


386B 38’ 

Air Wavelength (A) 

l Spectrum shifted for comparision 

figure 2 


Comet Bennett 1970 II 

Comparision Between Observation and Model 

Observation 1 

Model* 

R = 0.7856 AU. V = 26.64 km/s 


Air Wavelength (A) 

1 Aikman et at., Icarus 21, 1974, convolved with gaussian 

2 Spectrum shifted for Comparision 


figure 3 
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ABSTRACT 


Comets P/Brorsen-Metcalf (1989o), Okazaki-Levy-Rudenko 
(1989r) and Austin (1989c1) were all observed to have well 
developed H2O+ structure close to the nucleus that showed changes 
on a time scale of only a few minutes. We present some preliminary 
data on the observed ion morphology and their evolution. 

Introduction 

Wide field groundbased imaging has provided a wealth of data 
on the spatial and temporal distribution of cometary ion tails 
resulting from the interaction of cometary and solar wind 
environments (e.g., Brandt, 1982; Ip and Axford, 1982). Although 
spacecraft experiments have measured the plasma distribution in 
Comets Giacobini-Zinner and Halley (e.g., Galeev, 1987; Brandt and 
Niedner, 1987), they refer only to the times of the encounters. On 
the other hand, groundbased time series imaging of cometary ion 
features may provide data on the large scale ion structure dynamics. 

Comets P/Brorsen-Metcalf (1989o), Okazaki-Levy-Rudenko 
(1989r) and Austin (1989c1) were all dust-poor comets presenting 
an unusual opportunity to record ion structures in the coma 
relatively free of contamination by solar continuum. We took 
advantage of this fact to explore the temporal changes that occur at 
high groundbased resolution (~300 km per pixel) and time scales of a 
few minutes. 

Qb.SflIttatifl.Q5 

These CCD imaging observations were obtained with the 1.5 
meter Catalina Observatory telescope using an RCA CCD at 0.7 
arcsec/pixel and either a broad red-band (650-950nm) filter, or a 
narrow band filter transmitting the H20+ (0,8,0) band.at 619nm. The 
field was slightly less than 4 by 6 arcmin. The length of the imaging 






sequences were limited by the observing window and need for 
complimentary spectroscopic observations. 

P/Brorsen-Metcalf (1989o) was observed on 1989 August 23 
and 24 UT when r - 0.66 AU, D • 0.75 AU. and the phase angle was 
91<>. A continuous series of 5 minute exposures of over nearly 1.3 
hours was obtained the first night, and less frequently over one hour 
on the second night. A one hour series of images of Comet Okazaki- 
Levy-Rudenko (1989r) was obtained on 1989 November 22 UT when r 
■ 0.68 AU, D.» 0.58 AU, and the phase angle was 102°.Short 
sequences of Comet Austin (1989c1) were obtained over the nights 
of 1990 May 6-8 UT when r - 0.73 AU, O - 0.44 AU and the phase 
angle was 107°. 

Image processing 

The relative lack of dust in these comets improved the H20+ 
visibility compared to most other comets. Large intensity gradients 
from the nucleus condensation to the outer tail, and some 
contamination by C2 still caused some difficulty in seeing low 
contrast features We found that the best way to see changes of the 
rapidly moving ion features was to difference images in the 
sequence. This required careful registration by nulling the central 
condensation after shifting both images about the same amount to 
0.01 pixel. The resulting derivative image shows changes in location 
of the features. Tracking errors sometimes show up as imperfect 
cancellation of the central condensation and causes a slight 
broadening of the feature edges, but do not significantly alter their 
location or shape. 

Ion Feature Morpholog y 

The August 23 sequence on P/Brorsen-Metcalf was the most 
extensive time series and will be discussed in detail; he sequences 
on the other comets showed a variety of H 20 + morphology consistent 
with that in P/Brorsen-Metcalf. The primary H20 + features 
consisted of a central tail "core" flanked by "turning" rays which, in 
time, folded into the central tail. Virtually all features extend off 
the field edge 150,000 km from the nucleus. The apex of the main 
central tail edges was very close to the nucleus and merged into the 
seeing limited central condensation less than 600 km from the 
nucleus. The folding tail rays were fairly symmetric about the 
central tail axis in P/Brorsen-Metcalf, but significant deviations 


were seen in the other comets. Modeling is underway to determine if 
the intensity profiles are consistent with either isolated rays, or 
conical sheets. If related, the viewing geometry of the projected 
neutral current sheet may complicate the appearance of the rays. 

Evolution of the Ion Features 

The figure shows an excerpt of the P/Brorsen-Metcalf 
sequence showing the formation of a set of folding rays. In this 
case, a "cloud” of ions accumulated in a uniform coma around the 
nucleus to a radius of at least 17,000 km before forming 
identifiable rays initially 65° apart. At this scale, the rays then fold 
back and form the tail core in 40 minutes yielding an average folding 
rate of about 45° per hour. The folding rate is slower as it 
approaches the tail core, and will yield a lower rate in wide angle 
images. The rays appear to be essentially parabolic forms with the 
central condensation at the focus. A second set of rays formed in the 
same way 40 minutes later. 

The central tail core may appear to continue to stretch out 
behind the nucleus, and in some cases "detach” from the central 
condensation in a disconnection event. The first three frames in the 
figure shows a receding central tail segment that separated 15-20 
minutes before the sequence started. The apparent velocity is about 
60 km/sec. 

Summary 

CCD observations of rapidly changing H20+ morphology in 
P/Brorsen-Metcalf with high spatial and temporal resolution show 
40 minute production cycles of symmetric folding ray structure, 
extremely narrow central tail core extending to within 600 km of 
the nucleus, and a disconnected tail segment receding 60 km/sec 
from the nucleus. Comets Okazaki-Levy-Rudenko and Austin 
exhibited more complex asymmetric ray structure, but the motion of 
folding rays (observed over shorter time intervals) appeared 
consistent with those in P/Brorsen-Metcalf. More analysis will be 
required to see if these observations can provide constraints to 
models of ion tail formation and evolution, but qualitatively, it 
seems that the observed structure is consistent with models of 
draping solar wind magnetic field lines. The close proximity of H2O+ 
to the nucleus will certainly provide constraints on the ion 
transport mechanism. 
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Sequence of 5-minute CCD difference images of P/Brorsen-Metcalf 
on 1989 August 23 UT. North is to the right, east is to the top, and 
the sun is to the upper left. The long edge of each frame is 195,000 
km at the comet. The midtime of the light component of each 
difference image corresponds to 10:53, 10:59, 11:05, 11:11, 11:17 
and 11:22 UT. Note the faint star trail near the central condensation 
in the first three frames. 
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Walter Huebner 

If the DE is caused by the passing of a sector boundary, should this not influence the tail rays 
also? 

Stephen Larson 

I would think so. It is not clear that I observed a “classical* DE, but I believe that I have 
seen published examples where the main tail broke off, but “new" rays were not 
affected. 


Susan Wyckoff 

What is the location of the tail ray vertex in your H 2 0 + images? Was it on the sunward or 
the tailward side of the center of the HjO + brightness? 

Stephen Larson 

The tail rays appear to approximate parabolas of different sizes sharing the same focus at the 
central condensation. 


Anita Cochran 

You observed that the H 2 0 + seems to be starting very near the nucleus. In the past, obser- 
vations of C0 + have indicated this ion also started very near the nucleus. Do you have a 
feel for whether H 2 0 + and C0 + start around the same distance and if not, which forms closer? 

Stephen Larson 

At the spatial resolutions that I normally work with, CO + and H 2 0 + are co-spatial. 


Paul Feldman 

Can you comment on the C0 + /H 2 0 + ratio in these comets and why you chose to image in 
H 2 0 + rather than CO + ? 

Stephen Larson 

The C0 + /H 2 0 + ratio was low in these comets. Images of Comet Austin in 1989o using the 
IHW (International Halley Watch) CO + filter were very weak. 


William Jackson 

Have you done any image processing? Do you expect to obtain more information with image 
processing? 

Stephen Larson 

Various image processing algorithms have been applied to improve visibility of features in 
individual images, but simple image differencing has been the most useful to show changes 
with time. 


David Lynch 

Your subtracted images are very nice and show lots of interesting structures. Can you tell us 
how you plan to analyze these and what you expect to learn from these images? Also, given 
the extreme sensitivity of the morphology in the image to minute registration shifts, how are 
you determining the best regestration? 

Stephen Larson 

The difference-images taken together will provide velocity, and perhaps acceleration maps. 
Registration on the central condensation is very important. I used the imcntr task in IRAF 
with a 15 pixel window to find the fractional pixel peak and shift the image the appropriate 
amount before differencing. Remaining asymmetries are due to imperfect telescope tracking 
over the 5 minute exposures. Registration was tested by shifting 0.1 pixel in each direction 
to establish the goodness of fit. 


I. Konno 

Comment: Takao Saito has found a case in which Comet Halley did not show a DE when it 
crossed the magnetic neutral plane. He also found DEs which are not related to sector bound- 
aries. He proposed other models to explain DEs. DEs are not always necessarily related to 
the sector boundaries or the solar magnetic neutral plane. 

Stephen Larson 

I agree that the DE mechanism is not fully understood, but they seem to be triggered by 
some disturbance in the solar wind. 
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A TAIL-WAGGING EVENT IN COMET AUSTIN 
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ABSTRACT: On the morning of May 21, 1990 a sequence of images of Comet 
Austin was obtained in the light of . The filter was centered at 602.5 nm 

with FWHM of 5.0 nm. We were able to follow at least two waves out through the 
main ion tail structure. During the course of the observations, 2.5 hours, we saw 
two new tail rays form and undergo similar wave motion. Several condensations, 
possibly representing disconnected material, were seen to develop and move out 
along the tail ray with a velocity of about 60 km/sec. 

INTRODUCTION 

For the past 17 years the Joint Observatory for Cometary Research, JOCR, 
has been studying the large-scale interaction of bright comets with the solar wind 
and the interplanetary magetic field using photographic images of cometary ion tails 
obtained with a 14-inch Schmidt telescope and unfiltered IIa-0 plates. (See Brandt 
et al. 1975; Niedner, Rothe, and Brandt 1978; Niedner and Brandt 1978; Niedner 
and Schwingenschuh 1987). Recently we have installed a small, wide-field imaging 
system that uses a Tektronics CCD detector and narrow-band molecular emission 
filters, Comet Austin is the first bright comet available since the equipment has 
been installed. We were able to mount an extended observing campaign from the 
16th of April through the end of June, during the time we were on the mountain we 
observed more than sixty percent of the nights when the comet was observable. We 
report in this paper the results from one night, May 21st, in which a series of ion 
tails formed and evolved. There are severed other nights during which the ion tail 
was active, this paper is just the first in a series of studies of the variable structure 
within the ion, gas, and dust features. 

INSTRUMENTATION AND OBSERVATIONS 

Our optical system consists of the 300-mm f/2.8 Nikon lens originally used by 
Parker, Gull, and Kirschner (1979). A set of full-aperture molecular emission-line 
and continuum filters is available, The filters, their central wavelengths, and full- 
widths at half maximium are listed in Table 1. The detector is a cooled Tektronics 


128 


512 by 512 CCD camera system built by Princeton Scientific Instruments. As with 
the Parker, Gull, and Kirschner observations our system is mounted piggy-back on 
another telescope. The field of view of the CCD is 1.9 degrees square. Each pixel 
has a width of 13.6 arc seconds. On the morning of May 21 Comet Austin was at 
a distance from the Earth of 0.25 AU which translated into 2470 km per pixel as 
projected onto the plane of the sky or 2645 km projected onto the prolonged radius 
vector. 

On the morning of May 21, 1990 a series of 13 images was taken through our 
H 2 O+ filter. All exposures were of 10 minutes duration. Table 2 below lists the 
frame number, image sequence number, and the UT at mid-exposure. 

DATA REDUCTION 

The basic data reduction was carried out at the NRAO-VLA Array Operations 
Center using IRAF running on a Convex C-l and Sun workstations. The data 
reduction sequence consisted of 4 steps. The first was to subtract the average 
of the overscan columns (columns 1 and 2) to remove the bias level. This was 
supplemented by subtracting an average bias frame to remove any residual low- 
level fixed-pattern noise. The second was to divide by an average twilight sky flat 
field obtained through the fyO* filter. Nineteen sky flats taken during a string of 
several clear nights were used to make the average fiat. The third was to subtract 
the sky background level. A 5-by-5 pixel mode filter was applied to each frame. A 
two-dimensional suiface was fitted to each smoothed image and subtracted from the 
initial image. The final images have a background level within 10 counts of zero, 
The final step was to correct for atmospheric extinction. The night’s extinction 
was determined from aperture photometry of 32 stars which appear unblended in 
all frames and whose images had more than 5000 counts above the sky. Fitting 
a straight line to the magnitude as a function of airmass we obtained a slope of 
0.058 ± 0.023. Using this slope we corrected all of the pixels within each frame 
to relative intensity per pixel outside of the earth’s atmosphere. We have not yet 
transferred the data to an absolute system. The extinction correction ranged from 
0.08 magnitudes when we first obsered the comet to 0,05 magnitudes at dawn. 

DESCRIPTION OF THE INDIVIDUAL FRAMES 

The 13 frames taken with the filter have been made into a time-lapse 

movie showing the formation and evolution of the ion tail. Since, within the printed 
page, it is impossible to present a movie we provide a series of still photographs and 
a short description of the individual frames shown in Figures 1-1 through 1-13. The 
sequence number appears in the upper left corner of each frame. The images have 
been rotated so that the direction of the prolonged radius vector (PRV) is from 
left to right. The amount of rotation from a north up - east left orientation was 
23.3 degrees counter clockwise. The width of each frame is about 860,000 km when 
projected onto the PRV. We have placed the brightest spot within the head at the 
same relative position in each frame. The displayed intensity is the logarithm of the 
number of counts, permitting display of details near the head of the comet as well 


as in the tail. During the course of the observations three tails are seen, they will 
be referred to as T1 through T3. Some dust is being seen through the filter. 

The dust cloud developing to the northwest of the comet head, above tail Tl is the 
beginning of the dust tail that will dominate Comet Austin by the end of May. 

We will describe motions of the tails as if they are happening in the plane of the 
sky. It should be noted that we have limited knowledge of the true three-dimensional 
nature of the motions, 

• Frame It In this frame Tl dominates and is already curved, concave up, near 
the head. T2 is very faint, below and well separated from Tl. There is a hint of a 
condensation, Cl, in Tl at a distance of 146,300 km from the head. 

• Frame 2: Tl is bending down toward the position angle of the PRV. The con- 
densation Cl has moved away from the head an additional 8000 km. T2 appears 
slightly longer and is moving up toward the PRV. 

• Frame 3: Tl continues to bend down toward the PRV. It appears to have sharp- 
ened and split into two parallel tails at a distance of 186,200 km. We believe that 
this is the point of the condensation Cl seen in the previous two frames. T2 is 
lengthening and continues to move toward the PRV. 

• Frame 4: Tl continues to bend down toward the PRV, parts of it are now parallel 
to the PRV near the head. It has a diffuse appearance with a very faint hint of the 
condensation Cl seen at 252,700 km from the head. Again the point Cl appears 
to be the beginning of two very narrow and separate rays. T2 is longer still and 
continues to move toward the PRV. 

• Frame 5: This frame has some electronic noise but shows Tl bending down toward 
the PRV. T2 continues to lengthen and move toward the PRV while remaining 
straight. Due to the noise we are not certain of the location of Cl. 

• Frame 6: Tl is now parallel to the PRV for about one-third of its length, with 
the remaining length continuing to bend down. Another condensation, C2, in Tl 
has been identified 286,600 km from the head. The outer end of Tl is diffuse in this 
frame. T2 has begun to kink upwards above the PRV position angle near the head. 

• Frame 7: Tl and T2 appear to merge near the head. One gets the impression 
that Tl is lying down on top of T2. No condensations are obvious in this frame. 

• Frame 8: The declination tracking failed during this exposure but we can still 
see merging of Tl and T2. As a result of the tracking failure no condensations are 
visible. 

• Frame 9: Tl and T2 are merged for most of their lengths. A condensation, C3, 
has appeared in Tl at 401,000 km from the head, just ahead of the two bright stars. 
The merged portion of Tl and T2 are close to the PRV. 

• Frame 10: Tl and T2 continue to merge as T2 lengthens. The condensation C3 
has moved down the tail with a velocity on the order of 60 km/sec and will continue 
to do so for the next 3 frames. Another tail, T3, has begun to form on top of Tl. 

• Frame 11: Tl is beginning to bend upwards toward T3 and away from the PRV 
near the head of the comet. T3 is lengthening. C3 continues to move out along Tl, 
T2 has merged with Tl. 
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• frame 12: Dawn is beginning to affect our observations, but Ti and T3 are clearly 
bending up and away from the PRV. T3 is lengthing. 

• frame 13* The T1/T3 tail combination bends downward then upward near the 
head. It wo i appear that the process is about to begin all over again. 

COMET HEAD PHOTOMETRY 

We have measured the brightness, in instrumental magnitudes, of the comet 
head by doing aperture photometry on the individual frames. We present results 
at two extreme radii in Figure 2< The plus symbols show the brightness of the 
comet within a 54 arc second (4 pixel) diameter aperture as a function of tirhe, 
the crosses show the brightness within a 272 arc second aperture. The magnitudes 
have been scaled to two different, arbitrary zero points. The changes with time are 
approximately the same. We wonder if the apparent fall in brightness is associated 
with the formation of tail T2 and if the later increase in brightness is associated 
with the formation of tail T3. 

SPECULATION - TAIL WAGGING VERSUS DISCONNECTION? 

While watching the movie of these frames one is reminded of a flag waving in 
a stiff breeze. A better analogy might be a flag which has been released from its 
flag pole but continues to wave as it is accelerated in the direction of the breeze. 
Like the flag, the tail of the comet continues to wave, raising the question of the 
cause of the waving. Magnetism almost certainly provides the restoring force, but 
what causes the displacement? Is it due to irregularities in the solar wind, and can 
those irregularities cause nearly periodic waving? Or is the tail waving because of 
eddy shedding, just like the flag? 

High-resolution, large-scale displays of frames 9, 10, and 11 show discontinu- 
ities and condensations in the tail. Any of these frames viewed independently might 
be taken as evidence for a disconnection event. This leads us to speculate that the 
origin of tail wagging and disconnections could be the same. Could the trigger of 
tail wagging also trigger instabilites in the plasma tail, possibly associated with 
eddy shedding in the solar wind, waves in the plasma, or reconnection of the field? 
We hope that modeling of the tail motions in the movies of Comet Austin will allow 
some progress toward an answer. 
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Figures 1-1 to 1-3 (left column) and 1-4 to 1-6 (right column) - H 2 0+ Images of 

Comet Austin. Each exposure is 10 minutes. Frames are d-^cribed in the text. 
Frame 1-5 is dominated by noise. 
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Figures 1-7 to 1-9 (left column), 1-10 to 1-12 (right column), and 1-13 (bottom) 
Images of Comet Austin. 
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Table 1. Filter* for Wide- Field Imaging 


filter 

central A (nm) 

FWHM (nm) 

blue continuum 

445.0 

10 


514.0 

8 

h 7 o+ 

620.5 

5 

red continuum 

647.0 

8 


Table 2. 

HiO+ Frames from May 21 

sequence number 

JOCR archive 

UT at mid exp. 

1 

802 

08:28.72 

2 

803 

08:43.35 

3 

804 

09:01.17 

4 

805 

09:16.63 

5 

806 

09:31.22 

6 

807 

09:44.90 

7 

808 

09:56.37 

8 

809 

10:07.80 

9 

810 

10:20.37 

10 

811 

10:32.35 

11 

812 

10:43.82 

12 

813 

10:57.05 

13 

814 

11:08.03 
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Time, UT Hay 21, 1990 


Figure 2 - Comet Head Photometry. The plus symbols show the magnitude of the 
comet head inside a 54 arcsecond aperture, the crosses show the magnitude within 
a 272 arcsecond aperture, both at arbitrary zero points. 






Walter Huebner 

Are there any data for coincidence of sector boundaries with the DEs? If sector boundaries 
cause DEs, why do they not disrupt the tail rays? 

David Westpfahl 

The data are being tabulated now. 


Ray Russell 

What happens to the center of the coma (< 5 arc minutes diameter) during a DE? Does 
it brighten or not? If so, before, during, or after the DE? This information would be most 
helpful in the interpretation of our spectroscopy in the thermal IR. 

David Westpfahl 

The CCD Frames have been calibrated but brightnesses have not yet been determined. We 
plan to do this in the next few months. 
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Summary 

Outstanding disturbances of the plasma tails were observed in 1989 in three comets, 
Brorsen-Metcalf, Okazaki-Levy-Rudenko, and Aarseth-Brewington. Time variations of the 
tails were obtained from photographs provided by many astronomers. A two-dimensional 
MHD simulation was performed varying the speed and the direction of the solar wind flow. 
The simulation agreed quite well with the observations. Solar flares were identified as the 
sources of these disturbances. It was found that the sudden change in direction of the 
plasma tail axis occurs when the comet crosses a discontinuity surface of the solar wind 
structure accompanied by solar flares. 

1. Introduction 

A study of disturbances in plasma tails of comets is very useful in understanding the 
structure of the solar wind (Saito and Oki, 1987). We calculated how the plasma tail of a 
comet deforms with changes in the solar wind speed and direction, then we traced the time 
variation of the plasma tail disturbances using photographs sent to us from throughout 
Japan and compared them with the simulation. 

2. Tail Bending Phenomena and the Variation of the Solar Wind Flow 

During the last half of 1989 three active comets, Brorsen-Metcalf (B-M), Okazaki- 
Levy-Rudenko (O-L-R), and Aarseth-Brewington (A-B) appeared. Since the sun wets very 
' active during this time many disturbance phenomena of the plasma tails were observed. 

In many observed disturbance phenomena there are cases for which the tail is bent or 
the axis of the tail changes greatly in a short time as seen in series of photographs. Such 
phenomena are thought to be caused by the change of the speed and direction of the solar 
wind flow crossing the comet (Brandt et al., 1980). 

It is known that, in general, when there is discontinuity in the flow direction of the 
solar wind and the discontinuity surface is inclined relative to the flow, the flow direction 
is shifted as in Figure 1. The discontinuity surfaces inclined relative to the flow in the solar 
wind may be caused by a corotating stream or a solar flare. In" the case of the corotating 
stream, the change of the flow appears as in Figure 2 (Dessler, 1987). A corotating stream 
and a solar flare are similar to each other in that their discontinuity surfaces are both 
inclined relative to the flow so that they may be thought to produce a similar change in 
the flow direction relative to the discontinuity surface. 

3. Two Kinds of Plasma Flows 

In the case of our earth’s magnetosphere, the plasma flow is different on the inside 
and outside the magnetosphere. On the outside it is a solar wind flow which is nearly anti- 
sunward. On the inside it is plasma convection, which flows sunward, due to the Je x Bz 


force produced by the tail sheet current Je and the tail field Bz- 

In the case of a cometary magnetosphere, i.e., a cometary plasma tail, the outer flow 
is also the solar wind. The inner flow is^the plasma convection due to the Jc x Be force 
where Jc is the tail sheet current and Be is the magnetic field component perpendicular 
to Jc- However, if Bz and Be are in the same direction, then Je and Jc are in opposite 
directions (See Fig. 3 of Saito et al., 1987). Therefore the inner flow of the cometary 
magnetosphere is in the anti-sunward direction. A cometary plasma is accelerated due to 
the inner flow from the speed of a few km/s near the nucleus, eventually reaching the solar 
wind velocity (Minami and White, 1986; Tomita et al., 1987). When the direction and 
the speed of the solar wind changes, the shape of the cometary or earth’s magnetosphere 
will change and the deformation' travels leeward nearly at the solar wind speed. 

In this paper we solve the problem of deformation of a cometary magnetosphere, i.e., 
a plasma tail, and its transmission by the solar wind using simulation. Deformation of the 
magnetosphere is represented by deformation of the tail axis. Therefore, for the sake of 
computatinal convenience, we assume in this paper that the solar wind plasma “starts at 
the cometary nucleus” and flows with the given solar wind velocity. 

4. Tail Bending Simulation 

Pizzo (1989) has published a model of the variation of the solar wind parameters in 
the interaction region for a corotating stream (Figure 3) using a two-dimensional MHD 
simulation. We performed calculations to see how the plasma tail deforms, using the 
windsock theory, when a comet crosses the region where the solar wind has the speed and 
the direction as in Figure 3. In this model we assume that (1) a comet crosses a space with 
a structure in which the solar wind has the velocity as in Figure 3, and that (2) the plasma 
from the comet has the speed and direction of the solar wind at that moment, and that it 
flows away from the head keeping that velocity (Figure 4). We connect the the positions 
of the particles at a certain time that left the nucleus at a different time to determine the 
form of the plasma tail at that time. Then we trace how the plasma tail changes its form 
in time for a given velocity structure. In this way, we obtained the manner in which the 
deformation occurs for the cases with and without shocks (Brandt, 1970; Pizzo, 1989). 

Figure 5 shows the simulation for the case when shocks formed. We see that the tail 
axis shifts greatly at the shock pair and at the discontinuity surface between them. Also 
even a slight change of the direction of a few degrees with the distance from the comet 
head causes a large structure with an acute angle. Similarly, Figure 6 shows the simulation 
for the case when no shocks formed. Although not as rapid of a change as in the case of a 
shock, it shows a fairly large deviation of the tail axis. We also see that a large structure 
with an acute angle appears similarly to the shock case. 

5. Comparison with the Observational Results 

We chose examples (Figure 7) from the three comets and compared them with the 
simulations. The first example is a phenomenon in Comet B-M on August 13, 1989, in 
which the direction of the tail axis changes suddenly from north to south (Figure 8A). 
There is a solar flare which is thought to correspond to the disturbance on the above date 
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as in the case of Comet Bradfield, 19791 (Niedner et al., 1983). From the relative positions 
of the flare, the comet, and the earth, the shock and the discontinuity surface produced by 
the solar flare are expected to cross through the comet with an angle shown in Figure 9 
when the comet is seen from the earth. Thus we can interpret that the change of the tail 
axis from north to south in this phenomenon was caused by the comet crossing the solar 
wind discontinuity surface produced by the flare. The way the tail deforms resembles very 
much the results at the time the tail crosses the discontinuity surface. 

The next example is a phenomenon in Comet O-L-R on November 16, 1989, in which 
the tail is seen bent upward at 19:42 (Figure 10). Although there is a flare which is 
thought to correspond to this disturbance, from the relative positions in this case there is 
the possibility that we are looking at the discontinuity surface from a rather large angle 
from the earth, thus it did not look as pronounced as in the case of Comet B-M. 

The last example is a phenomenon in Comet A-B on December 10, 1989, in which we 
see how the tail axis moves from upward to downward (Figure 11 A). This disturbance is 
similar to the first example (See Figure 11B). 


6. Conclusion 

The results obtained by tracing the deformation of the plasma tails of comets using 
the solar wind structure with a two-dimensional MHD model agreed very well with the 
observed tail disturbances. There exist solar flares which caused the disturbances in the 
plasma tails of comets. From the relative positions, the sudden change in the direction of 
the tail axis is considered to happen when the comet crosses the discontinuity surface of 
the solar wind structure produced by the solar flare. 
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FORMATION IN THE COMA FROM PHOTODISSOCIATION OF ICE 

GRAINS 
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ABSTBAgT 

Long ago visual observations of comets suggested that 
there are jets in comets but it has only been recently that 
A 'Hearn and his coworkers have proven that some of these 
jets are due to emission from the CN radical [1,2,3]. 

Recent studies in our laboratory have shown that CN radicals 
can be ejected directly into the gas phase from the 
photolysis of frozen vapors if the parent molecule has been 
excited to a repulsive excited state. This later 
observation suggest that the jets that have been observed 
may be due to photodissociation of icy grains in the coma. 

In this presentation a theory of radical formation from icy 
grains will be presented. It will be shown that direct 
formation of free radicals in the coma is an effective way 
to produce radicals from icy grains in the coma. The model 
predicts that icy grains could produce from 6 to 800,000 OH 
radicals/s per grain depending upon whether the radius of 
the grain is 0.3 to 100 urn. 
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Lewis Snyder 

What is the structure of NCCN? 

William Jackson 
It is linear. 


Ray Russell 

Would small silicate grains (strings of 100 - 200 A spheres, stuck together in long twisted 
strings like IDP’s) coated with organics provide a viable source for the radicals? 

William Jackson 

If we had 100 A spheres each individual sphere would not give you enough of a yield to explain 
the jets. If thousands of these spheres are together then one might be able to obtain a high 
enough yield to explain the observations. 


Carey Lisse 

Were the production rates quoted in the final table for 1 grain or a jet of an estimated mass, 
made entirely of one given particle size? Also, what ambient grain temperature and pressure 
have been used for these calculations? 

William Jackson 

The production rates were total production rates for all grains of a given size. We used the 
flux measured in Halley for a given grain size. We used a grain composition of 2 % HCN, 2 % 
CjNj, and 90% HjO. 

No ambient grain temperature or pressure was needed for these calculations. 
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' AN ALTERNATIVE TO THE DIRTY SNOWBALL MODEL OF COMETS 



Tom Van Flandern 

VF Associates / 6327 Western Ave, NW / Washington DC 20015 


A. Comparison of Models 

To be viable, a scientific theory or model must not be 
contradicted by observations or experiments. To be useful, a model 
must explain some phenomenon and aid our understanding of it. To 
be valuable, a model must make successful predictions. 

By these criteria, the "Dirty Snowball" model for the nature 
of comets is now on shaky ground. By contrast, the "Satellite 
Model" easily passes all these measures of the worth of a theory. 
Although the Dirty Snowball model is familiar to all comet 
investigators, the newer Satellite Model may not be. 

Briefly, in the Satellite Model, 113 comets consist of a cloud 
of asteroid-like debris of all sizes, which is gravitationally 
bound together. Typically, the comet's sphere of influence, within 
which it can keep material in stable satellite orbits against the 
attraction of the Sun, is dominated by a primary "nucleus", but 
also includes bodies of kilometer, meter, millimeter, and dust 
size, all of which contain volatiles. Comet tails are formed when 
solar radiation pressure drives gas and dust away. In the Dirty 
Snowball model, that material must be driven off the nucleus. In 
the Satellite Model, most of it is spread throughout the sphere of 
influence, and needs only to escape. 

Let me briefly review each important physical characteristic 
of comets, and how it fits into each model. Surely a good model 
will automatically explain important physical properties. A 
failing model often needs to be patched up to accommodate new 
findings. 


B. Albedos and Spectra 

In the Dirty Snowball model, the chief ingredient of comets is 
water-ice. In the Satellite Model it is carbonaceous rock. In 
either model the surface is coated with an extremely black 
carbonaceous residue. The Dirty Snowball model does not expect 
this black coating; it merely accommodates it. The Satellite Model 
anticipates that the same carbonaceous material which coats 
asteroids and so many solar system moons also coats comets. Comets 
have a mean albedo of about 4%, primarily from this coating. Other 
solar system moons which are icy underneath their carbonaceous 
coatings have mean albedos of about 20%; whereas asteroids, which 
^re rocky underneath, have mean albedos of about 4%, the same as 
comets . 

The spectra of comets are very similar to those of asteroids. 
In the Dirty Snowball model that is coincidental, due to the 
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surface coating. In the Satellite Model, it is because comets and 
asteroids are made of the 6ame basic materials throughout. The 
absence of ices in the reflection spectra of comets 123 is a 
noteworthy feature, arguing against the Dirty Snowball model. In 
Comet Halley, heavy molecules to 120 amu were detected 153 : fairly 
heavy stuff if comets are primitive objects, but quite expectable 
if they are like asteroids. 

C. Densities and Strengths 

We have no direct measurements of the densities or tensile 
strengths of comets. The Dirty Snowball model requires that both 
be low, partly to permit the solar non-gravitational forces to be 
as strong as is observed; and partly to account for the ease with 
which comets split, and the low accelerations and high surface 
brightnesses of the fragments. It is concluded that comets must be 
either flat and pancake-like (ruled out by spacecraft missions), or 
fluffy dustballs M3 . The Satellite Model expects asteroidal 
bodies with densities in the 2-3 g/cc range. In it, the high non- 
gravitational forces are brought about by the large surface area of 
the sphere of influence, typically 10,000 times greater than the 
area of the nucleus alone. All of the satellite bodies and dust 
within the sphere of influence (not just the nucleus) feel the push 
of solar radiation and affect the motion of the center of mass. 
Since a split event is just the gravitational escape of a 
satellite, and nothing actually breaks up, comet materials of 
normal strength are expected. The indirect observational evidence 
which is available tells us this: Micro-meteoroids collected in 
the upper atmosphere during meteor showers associated with comets 
are primarily of carbonaceous chondritic (normal density) 
material- 53 . And fireballs associated with comets (because their 
orbits have aphelia far beyond Jupiter) survive to low altitudes in 
the Earth's atmosphere , C6] again suggesting they have the strength 
of rock, not fluffy dustballs. 

D. Radar Reflectivity, Lightcurves, Comas 

The radar reflectivity of comets is around 10%. This implies 
a non-porous surface, probably rocky, C73 as expected by the 
Satellite Model. The Dirty Snowball model specifically expects a 
porous surface. 

Comet lightcurves are often complex. The well-observed 
lightcurve for Comet D 'Arrest is of a distinctive triple-maximum, 
triple-minimum type associated with binary stars and asteroids. 1 -® 3 
In a few comets, multiple or "granulated" nuclei have been directly 
observed even ift cases where no "split" event followed. 

The Dirty Snowball expects the comas of comets to be formed as 
they approach the Sun, and grow in size as they get closer. The 
Satellite Model predicts that the comet's sphere of influence, and 
therefore the size of its ever-present coma, will generally shrink 
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as the comet approaches the Sun. Observations show that comet 
comas are already formed even when comets are first detected in 
excess of 10 au from the Sun, and generally diminish in linear size 
as they approach the Sun. 

E. New versus Old Comets 

Comets "new" in the Oort sense lose 2-3 magnitudes in absolute 
brightness on their first passage near the Sun. Old comets show no 
further absolute brightness losses, in a few cases even after 
dozens of revolutions. The Dirty Snowball model predicted that the 
brightness losses would continue, but at a reduced rate, due to the 
formation of a surface “frosting"; and that Comet Encke would be 
considerably diminished if not totally evaporated by this time. 
The Satellite Model predicts exactly the observed behavior: As 
comets approach the Sun for the first time, their spheres of 
influence shrink to a smaller diameter than ever before, causing 
the gravitational escape of much material which was previously 
gravitationally bound inside the coma. But on subsequent 
revolutions, as long as the comet's perihelion distance is not 
reduced further, a comet can lose only gas and small dust particles 
from solar radiation pressure. The millimeter and larger size 
bodies are virtually unaffected, insuring that the comet will 
maintain its general appearance at least until solar radiation 
pressure successfully removes all gas and dust from the coma. 

Because there are numerous smaller spheres of influence within 
the primary one in the Satellite Model, this process of losing dust 
and gas can take many thousands of years. Moreover the volatiles 
are not lying in layers of ice near the surfaces of these satellite 
bodies. Instead, they are interstitial, just as is found in 
carbonaceous meteorites on Earth, which are about 15% interstitial 
water by volume. The Satellite Model obviously predicts that "new" 
and "old" comets will not differ in albedo. The Dirty Snowball 
model expects an albedo change when the surface "frosting" forms 
just after the first perihelion passage. Observations should soon 
tell us which prediction is correct. 

F . Jets 

In the Dirty Snowball model, volatiles must be ejected from 
the surface of the nucleus in high-speed jets to account for their 
abundances and relative independence of solar distance. Although 
the Satellite Model can accommodate jets, it does not predict them. 
Spacecraft encounters with Comet Halley photographed what appeared 
to be jets streaming sunward from the nucleus. However, unlike the 
expectations of the Dirty Snowball model, dust in the jets appears 
to be unaccelerated, and there is no evidence for pressure in 
cometary jets above the ambient pressure, and no supersonic flow in 
them. Moreover, bursts of dust particles were observed to come 
from sources in the coma other than the nucleus, as the Satellite 
Model would certainly predict. 



Are the spacecraft-observed features really "jets"? The 
Satellite Model suggests another explanation of the pictures. The 
"jet" features eminate from surface spots of high albedo* One 
might have expected true jets to eminate from the hotter , darker 
regions, rather than the cooler, brighter ones. In the composite 
picture of Comet Halley, the "jets" look instead very much like 
sunlight back- scattered into the dust. But in the single frames, 
the "jets" seem concentrated in narrow, sunward beams. This seems 
likely to be the result of the well known "opposition effect" for 
asteroidal material, in which most of the reflected light is 
concentrated in a cone with a divergence of just plus or minus 3 
degrees from the Sun direction. So the "jets" may be just dust 
backlit by sunlight over surface bright spots. A single "jet" not 
pointing exactly sunward, or not eminating from only the brightest 
surface spots, or visible in direct (not backlit) sunlight, would 
provide a counter-example to this proposal. 

G. Split Comets 

The Dirty Snowball model did not expect comets to split, but 
accommodates that by assuming low strength materials and explosive 
jets of volatiles. The Satellite Model expects major satellites to 
drift off on independent solar orbits if the comet's sphere o,f 
influence is reduced to smaller than the satellite's distance from 
the primary nucleus. These events would have the appearance of 
"splits", but could occur only on the inbound leg, or while the 
comet was close to perihelion. While the comet's distance from the 
Sun was increasing rapidly, its sphere of influence would also, 
preventing gravitational escapes or "splits". Observations show us 
that all 14 split events where the date of split is uniquely 
determinable occurred in the range from three years before 
perihelion to 10 days afterwards . C9] Of 18 other cases with too 
few observations to make a split date determination, 8 appear to 
have happened after perihelion. But considerable uncertainly 
exists in each of these cases; and it is impossible to distinguish 
a true split event from a chance alignment of well -separated comet 
components in the line of sight due to viewing geometry. 

Perhaps the single strongest observation which distinguishes 
these two models is that of the relative fragment velocities 
resulting from split events. The square of relative fragment 
velocities must be proportional to the energy available no 
accomplish the split (from 2E = m v 2 ) . In the Dirty Snowba .l 
model, if the energy comes from solar forces which are inverse 
square (light, heat, gravity, radiation), then relative fragment 
velocities will be proportional to solar distance to the minus one 
power (b = -1). If the energy comes from entirely internal sources 
unaffected by solar energy, then relative fragment velocities will 
be independent of solar distance (b = 0). The Satellite Model 
makes a very specific prediction about relative fragment 
velocities, since they are gravitational escapes from the sphere of 
influence: relative velocity is proportional to solar distance to 
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the minus one-half power (b = -0.5). Observations (see the 
previous Sekanina reference) over two orders of magnitude in solar 
distance for comet splits show that b - -0.57 ± 0.10. Whereas this 
result is in precise agreement with the Satellite Model, it would 
seem to rule out the Dirty Snowball model with more than four 
standard deviations of statistical significance. 

H. Conclusion 

In conclusion, the Dirty Snowball may no longer be a viable 
model, even after its patches have been applied. Although It 
appeared for a while to aid understanding, its predictions have not 
been noted for their success; and it is now flatly contradicted by 
one line of evidence. The Satellite Model provides a whole new 
understanding of comets which happens to describe all of their 
physical properties rather well. 

We cannot make progress in understanding comets until we get 
the model right. Research interpreted with the wrong model may be 
doomed to early obsolescence. So the question is one of 
importance. Let us therefore give the matter our most critical 
thinking, and our severest judgment. Can we move forward and gain 
understanding of comets without a change of models? 
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William Jackson 

How do you continue to have H, OH, CH, Cj, C 3 , etc., in short period comets, e.g., Comet 
Encke? 

Tom Van Flandern 

The volatiles are not from surface ices, as determined by the reflection spectra. Since one must 
go below the surface to "boil off” volatiles, this is easily accomplished for rocky, chondritic 
material. Even carbonaceous chondrite meteorites still have about 15% interstitial water by 
volume. 


Ray Russell 

How does the satellite model explain the amount and size distribution of dust being released 
after so many apparitions as P/Halley and P/Encke have had? Wouldn’t all the dust be 
swept out or wouldn’t your model require such a high density for collisions to occur so that 
the whole gravitational sphere would be optically thick (in which case the pieces would be seen 
in mm radio and in radar observations) and they would grind themselves to dust very quickly? 

Tom Van Flandern 

The preservation of the dust is no more difficult than in the Dirty Snowball model, because 
one has thousands of source bodies gravitationally binding dust, not just the nucleus. Just as 
for zodiacal dust, the time scale to remove dust depends on how close the comet gets to the 
Sun, how long it stays there, and the size distribution and lifetimes for the dust particles. In 
the Satellite Model, a continuum of sizes of dust particles are present, including sizes large 
enough to be very little influenced by solar radiation pressure. Russell’s question seems to 
suppose that the gravitational binding forces in comets are too small to hold the dust for long 
against radiation pressure. But some of these gravitational binding forces are strong enough 
to hold larger dust particles indefinitely. As for the grinding of the smaller dust, Russell’s own 
paper noted that the dust is not in thermal equilibrium, and Konno and Huebner discussed 
evidence of dust fragmentation in the coma, such as an increase in the number of small grains 
away from the nucleus. 


General Discussion: 

Many participants objected to the Satellite Model that had just been presented. It was felt 
that this model was untenable for several reasons. Some of the objections to the model in- 
cluded: (1) The lack of evidence of satellites of Comet P/Halley in the images from the Halley 
Multicolour Camera or the Vega cameras, (2) the difficulties that the model has in explaining 
the observed gas and dust production as a function of heliocentric distance which are easily 
explained with Whipple’s Icy Conglomerate Model, (3) the statistics used for comet splitting, 
(4) the inconsistency with the persistence of dust over many revolutions of the comet around 
the Sun, and (5) inconsistencies with Sun-grazing comets. Tom Van Flandern replied with 
qualitative counter arguments to these objections, but the real test has to be a quantitative 
model based on physics and chemistry and its comparison with observational comet data. 
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Summary 

We have developed a hydrodynamic model of a spherically symmetric dusty gas flow in 
a cometary atmosphere assuming a single-fluid, inviscid, perfect gas. The hydrodynamics for 
gas and dust, which involves the gas drag force (momentum transfer), heat exchange between 
gas and dust, photodissociation energy for H 2 O gas, and radiative heating and cooling terms 
for dust particles, are solved using the Gear method for stiff, coupled differential equations. 
Calculations have been done with a dust size distribution for radii a = 0.01/zm to 10 cm with 
densities variable with the size. A nucleus size of 4.0 km radius with a density of 0.5 g cm -3 
and a toteil dust-to-gas mass ratio y = 1 were adopted. 

There are indications from in situ observations that dust particles fragment into smaller 
ones. Fragmentation of dust partilcles has been incorporated into the model. This is done 
by adding source and sink terms in the continuity equations for the dust. Lifetimes for the 
decay of dust particles were assumed as a function of particle size. It is also assumed that dust 
particles always fragment only into the next smaller size. 

1. Introduction 

Recent observations of dust in comets suggest fragmentation of dust particles at some 
distances from the nucleus. Fragmentation of dust is implied mostly from the relative increase 
in the number of smaller grains or the decrease of the intensity of scattered light less rapid than 
l/R, which is expected for dust with constant speed and no fragmentation. From Giotto’s PIA 
and DIDSY data, McDonnell et a!. (1987) found some evidence for possible fragmentation 
or evaporation of grains. Vaisberg et al. (1987) obtained several pieces of evidence for an 
evolution of the mass spectrum and the supply of new particles in the low-mass range while 
dust propagates from the source. There are also some ground-based observations suggesting 
dust fragmentation. It is thought that fragmentation may arise from the sublimation of a gluing 
material such as POM (polyoxymethylene) in the grain agglomerates of the dust (Boehnhardt 
et al., 1990). 

We study the radial distribution of dust by solving hydrodynamic equations. In order to 
model the dust fragmentation in a comet, we assume the dust size distribution at the source 
(surface of the nucleus) and the lifetime of grains which is a function of the particle size. 

2. Modeling 

We solve sets of hydrodynamic equations for a dusty gas flow in a comet coma which 
assumes single-fluid, inviscid, perfect gas with the Gear method, a stiff differential equation 
solver. The equations contain the gas drag force on dust, heat exchange between gas and dust, 
and photodissociation energy for H 2 O gas. 


Calculations have been done with dust particles of sizes from O.Ol fiiti to 10 cm in radius 
with a continuous distribution approximated by 22 discrete particles sizes seperated by a 
logarithmic scale so that a particle always fragments into ten particles of next smaller size. It 
is assumed that the smelliest particles have infinite lifetime against fragmentation. The size 
distribution at the nucleus surface is given by 

Pa ~ Pmaxifll Q-max)^ r 

where p a is the spatial mass density of the particles of radius a and p m ax is the density of the 
particles of the maximum radius a max that can be entrained by the gas. The bulk density of 
dust particles is a function of the dust size: 

pbuik = 2.2 - lAa/(a 0 -fa) g cm -3 , 

where a 0 = 10 -4 cm (Lamy et al., 1987). We assumed for the nucleus of Comet Halley that 
the radius R n — 4.0 km, the density p n = 0.5 g cm -3 , and the dust-to-gas mass ratio \ = 1. 

3. Results and Discussion 

Figures 1 to 3 are obtained from calculations with a = —2.0,0 = 1.0, and r max = 300 
s. The gravity between the nucleus and the dust has been taken into account but the results 
show that the gravity of a Halley-class nucleus has very little effect on particles smaller than 
about 1 cm. 

Figure 1 shows the average velocity of dust particles of different sizes as a function of 
cometocentric distance. The farther out in the coma particles fragment, the less acceleration 
the fragmented particles experience. Therefore the mass-averaged velocity of the fragmented 
particles starts decreasing at the distance of the scale length of their parent particles. 

Figure 2 shows the dust size distribution at various cometocentric distances. The number 
density is plotted as a funcion of dust particle radius. The relative density of larger particles 
decreases with the distance because of fragmentation. The number density of small particles 
could increase with the distance because of rapid fragmentation despite the R~ 2 dependence 
of the number density when there is no fragmentation. 

In Figure 3, the cumulative number density of dust particles at the surface of the nucleus 
( R = 4.0 km), and at R = 1,000 km and R = 3,000 km is shown. The result (Lamy et al., 
1987) from observations by the Vega 2 spacecraft is also shown. The spacecraft results are 
deduced from observations made at much farther distances over- several thousand kilometers 
assuming an R~ 2 dependence of the density. Therefore the spacecraft data may not give a true 
distribution at R = 1,000 km if there is fragmentation within the these distances. 

In case of dust fragmentation, the dust size distribution at the nucleus of a comet may not 
be uniquely determined without the knowledge of the lifetimes of all kinds of dust particles. 
Spatial profiles of narrow band continua should give information on the scale lengths of the 
parent particles. Lifetimes of dust particles may be obtained from the scale lengths. 



COMETOCENTRIC DISTANCE (km) 

Fig. 1.- Mass-averaged dust velocity. Average velocity decreases because the particles 
fragmented at larger distances have lower velocities (typically that of their parents). 
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Fig. 2 - Dust number density at various cometocentric distances as a function of the 
particle size. Because of rapid fragmentation, the number density of small particles can increase 
at larger distances despite the R~ 2 dependence of the density in the case of no fragmentation. 







Fig. 3.- Cumulative number density of dust particles as a function of mass. The number 
density shown indicates the total number of dust particles of the equal or larger mass at the 
point. The profile ded.uced from the Vega 2 observations is shown by a dotted-line. 
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D. Boice 

What is the Gear method? 

I. Konno 

It is a method to solve stiff differential equations. 


P. Weissman 

The Vega 2 encounter was at 8,000 km. Why do you compare it with your model curve at 
1,000 km? 

I. Konno 

Lamy et al. [Astron. Astrophys. 187, 767 (1987)} deduced the dust size distribution at 1,000 
km from the Vega 2 data. 


Carey Lisse 

Do you plan on extending your dust size distribution model to cometocentric radii layer then 
10 3 km? 


I. Konno 

Yes, we plan to calculate much farther with the one-dimensional model, but we need to use 
a two-dimensional model when solar radiation pressure on dust becomes significant. 


Paul Feldman 

Does the Giotto data show a variation in the dust size distribution with distance from the 
nucleus? 

Ichishiro Konno 

It may, but the published Giotto data gives only the dust size distribution averaged over a 
few thousand kilometers. 


Ray Russell 

Do you calculate the IR emission expected from your model? Thermal images may provide 
the constraint you need to get a unique determination. 

Ichishiro Konno 
No, not yet. 
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Abstract. The spacial distribution of neutral species in the outer coma of a comet is de- 
scribed using a modified version of the “fountain” model. This model accounts for neutrals 
that are produced in multiple generations after release from the nucleus (or production in 
the collision zone of the inner coma) by an iterative application of the “fountain” model. 
A preliminary description of the model is given; however, results will be presented in a fu- 
ture paper. Suitable applications of the model include the following: (1) distributed coma 
sources of neutral gas-phase species from anisotropic dust emission, (2) the description 
of cometary comae at large heliocentric distances, and (3) day /night asymmetry of gas 
production from the nucleus. 

Model Description 

Two papers that develop the theory of the “fountain” model by Haser (1957, 1965), 
form the basis of the present multigeneration model. The neutrals flow freely outside 
of the collision zone, the extent of which defines the inner boundary of our model, R c , 
approximately 10 4 km for Comet Halley at the Giotto encounter. The outer boundary of 
the model, R g , occurs at approximately 10 7 km in this case, where significant curvature 
of the particle’s trajectory is introduced by the effects of gravity. The inner boundary 
conditions for the densities and velocities are given by our spherically symmetric, fluid 
dynamic model with detailed gas-phase chemistry (Schmidt et al. 1988). 

Each neutral species is characterized by a lifetime against photo destruction (photodis- 
sociation, photodissociative ionization, find photoionization) and charge exchange with the 
mass-loaded solar wind. The excess energies of the dissociation products are included to 
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describe the formation of the following generation of neutrals. We assume that the coma 
is optically thin to visible and UV sunlight, that the neutrals do not undergo collisions in 
the outer coma (except with the mass-loaded solar wind), and that the conditions can be 
described in the steady-state approximation. The effect of radiation pressure from sunlight 
on the different neutral species is incorporated into the model. Each species is character- 
ized by a radiative acceleration parameter, a, that depends on its unique ability to interact 
with light, folded with the solar spectrum. This effect results in an axisymmetric distri- 
bution for the neutrals. Typical reciprocal lifetimes against photodissociation, Tpp, and 
photoionization, 7>;, from Schmidt et al. (1988) and values of a from Konno (1990) for 
important cometary neutrals are presented in Table 1. Also given are stand-off distances, 
z max (the maximum distance to which the species can travel along the Sun-comet axis), 
assuming an initial velocity of 1 km s -1 . 

In addition to the nucleus and inner coma sources of neutrals, each point within the 
region bounded by the model acts as a source for a subsequent generation when a neutral 
of the previous generation is destroyed. Just prior to destruction, particles at any point 
within the outer coma have two peculiar velocities with respect to the nucleus, as given by 
the original “fountain” model. This results in a skewing of the parabaloid which bounds 
the set of all possible trajectories of the isotropically emitted, next generation species. 
When all of the important reaction pathways that lead to the creation of a given neutral 
are considered, each spacial distribution that results can be superimposed to give the total 
density distribution. This feature of the model is particularly useful for isolating the most 
important reactions that contribute to a given species as a function of cometocentric and 
heliocentric distances. 

Computational Grid 

The computational grid is axisymmetric, With the Sun-comet line as the axis of symme- 
try. It has been described by Schmidt and Wegmann (1980) with parameters appropriate 
to the bounded region discussed above. This grid was chosen for computational efficiency 
and compatability with our coma chemistry code with solar wind interaction (Schmidt et 
al. 1988). A complete solution can be obtained along the Sun-comet line with a double 
integration in a half-plane bounded by this axis of symmetry. The solution for a general 
point off-axis can be found with a triple integration in a half- volume bounded by this plane. 


Discussion 


Whenver possible, we will incorporate data for determining radiative acceleration and 
photo lifetimes for interesting cometary neutrals. This includes the water group species 
(such as, OH and O) and other species (such as, C2, C3, NH 2 , and CH) that have short 
stand-off distances. Preliminary results along the Sun-comet axis have compared favor- 
ably with our previous radially symmetric calculations. Possible applications of the model 
will be the following: (1) distributed sources of neutral gas-phase species in the coma 
from anisotropic dust emission, (2) cometary comae at large heliocentric distances, and 
(3) day/night asymmetry of gas production from the nucleus with possible compositional 
anisotropies. The results can be applied to observed isophotes for cometary neutrals as 
well as used for computing two-dimensional, ionic source terms for cometary plasma mod- 
els. 
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the computational grid used in the model and gave valuable comments. This study was 
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Table 1 


Radiation Acceleration and Photolifetimes of 
Important Cometary Neutrals 
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2.1 - 10” 3 


1.7 • 10®km 

5.5 - 10 6 

1.0 • 10 * 

8.9 - 10 7 

6.6 - 10* 

2.8 • 10® 

9.1 • 10 4 

1.9- 10® 

6.2 • 10 4 

3.6 • 10® 

2.9- 10® 
1 . 6 - 10 ® 

1.7 ■ 10® 
1.8-10® 

1.1 • 10 ® 

2.4 • IQ 7 


1.17 • lO ^ s" 1 

5.00 - 10 “® 

1.22-10“* 

3.17 - 10“ 7 
1,76 • 10 “® 

1.00 - 10 “® 
2 . 10 - 10 -4 
1.40 • 10“ 7 
4,16 • 10“ 3 
2.00 • 10 "® 
7 . 11 - 10 “® 
2.78 - 10 “* 
4.00 - 10" 4 


1 . 10 - 10 “® 
6.55 • 10~ 7 
3.80 • 10“ 7 
6.10 • 10“ 7 


9.10 • 10" 7 
7.56 • 10“ 7 

3.60 • 10" 7 
4.04 • 1G“ 7 

7.31 • 10“® 
3.58 • 10“® 
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CO is hardly a third generation species, particularly for R > 10 4 km. 

Dan Boice 

I considered CO a third generation species because of the distributed source found by Eber- 
hardt et al. [{Astron. Astrophys. 187, 481 (1987)] using Giotto NMS data that existed out 
to about 20,000 km. I assumed that CO comes from a heavier molecule released by the dust. 
Additionally, on a larger distance scale (R ~ 10* km) COj photodisociates into CO. 
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^ MONTE CARLO MODEL OF A COMET COMA 

(j y Anthony J. Ferro 

Department of Physics and Astronomy 
Arizona State University 


Introdugtign 

A steady-state multi-generational Monte Carlo model has been developed to describe 
the distribution of a neutral coma species which has been produced after several photodis- 
sociation steps. This distribution can be compared with either line profiles derived from 
long-slit spectroscopic observations or narrow band imaging, and be used to determine the 
elemental abundances present in the comet coma. A Monte Carlo model was chosen due 
to limitations in standard forms of more analytic models, such as those of Haser (1957) 
and Festou (1981). In contrast to the analytical models, the Monte Carlo method can be 
modified easily to include new physical parameters. 

The Multi-Generational Model 

The model can be described using the example of the production of NH in comets by 
the photodissociation sequence, 


NH 3 -♦ NH 2 -♦ NH -♦ N. (1) 

The grandparent molecules, NH 3 , are ejected radially from the nucleus, either symmet- 
rically, or asymmetrically. After a characteristic time, Tnh 3 ~ 7700 s, the solar EUV 
radiation produces NH 2 , with the excess photon energy partitioned into the kinetic and 
internal energies of the fragments, NH 2 and H. In the current model, the ratios of the 
energies are set at 9:1. Momentum is conserved in each dissociation, while the directions 
of the NH 2 fragments are randomized, in the center of mass (NH 3 ) frame. This process is 
repeated for NH 2 and NH, using calculated photodissociation rates. 

Figure 1 shows a flowchart outlining the model program. Note that the basic structure 
consists of two nested loops, one for each test particle and one for each dissociation. 
Typically, each model is ru,n for 10 6 to 10 8 iterations, with three or four dissociation steps 
for each particle. All calculations are done in three dimensions, centered on the comet 
nucleus. The final output is a two dimensional map of the observed molecular species. 
Spatial profiles can be taken from this map to compare with observed spatial profiles, or 
the map may be directly compared to narrow-band imaging data. 

The current mode) allows simulation of various types of asymmetric ejection of the 
initial parent from the comet nucleus. The asymmetry can range from simple hemispheric 
ejection to ejection from very narrow regions, simulating jets. The model currently also 
can simulate accelerations due to solar radiation pressure. For certain molecules this 
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radiation pressure, derived from the g-factors, can play a major role in the morphology of 
the distribution of that molecule and of their daughter products. 

In cases where there are several pathu for producing a particular molecule, a suite of 
models may be run, and the results scaled then added. For instance, to model the NH 
distribution, the sequence in equation (1) is the primary process, with a branching ratio 
of 97%, 100% and 100%. However, the photodissociation sequence. 


NH 3 NH -* N (2) 

with a branching ratio of 3% and 100%, can also be considered. A model can be run for 
each sequence, then the resulting maps multiplied by the appropriate ratio and added. 

Result? 

All results presented were run on the Arizona State University Supercomputing Fa- 
cility CRAY-XMP/18se. Typical run times for a single model ranged from 20 to 120 cpu 
minutes. 

Figure 2 shows a map of NH 2 generated with the Monte Carlo model in which so- 
lar radiation pressure has been taken into consideration. The solar acceleration is 0.17 
cm/sec 2 . Each side of the map represents 5 x 10 5 km. 

Figure 3 shows a map of NH generated with the model. Again solar radiation pressure 
has been considered. Each side of the map represents 2 x 10 6 km. The wide scale of NH 
is due to its long lifetime. The lifetimes used in this model are given in Table 1. 

In figure 4 we show the NH distribution observed in P/Halley, compared with the 
Monte Carlo model. Data are from S. Wyckoff (see article in these proceedings). In this 
case the model was symmetric and azimuthally averaged to increase signal to noise ratio. 
The fit, especially at distances greater than 10 6 kms. is quite good. 

Future Work 

Planned future enhancements to the model include: timed release to model a rotating 
source, inclusion of collisions by suitable modifications to the expected lifetimes close to the 
nucleus, addition of distributed sources, and the inclusion of non-gravitational forces. It is 
also hoped that the code can be modified to run in reasonable times on smaller computers, 
such as the current generation of high speed workstations. 



Festou, M.C. 1981, Astr. Ap ., 95 , 69. 

Haser, L. 1957, Bull. Acad. Roy. de Belgique, Classe des set., 5th ser., 43, 740. 
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Initialize Variables 


For Each Particle: 

N«1 to # particles, do: 


Collapse from 3-D to 2D. 
Output Image. 


Eject particle: 

Assign random ejection time 
Assign random ejection directic n 
Assign ejection speed 


Convert from spherical to rectangular 
coordinate system. All coordinates 
are cometocentric 


Assign Directions: 
Random, but opposite. 


Calculate Disassociation Energy: 

Lambda based on cross section. 


For each disassociation: 

D=1 to # of disassociations, do: 


Calculate Change in Velocity: 

Momentum is conserved. 


Determine Particle State 
at Time of Observation: 

May be in any parent or 
daughter state. 


Add Velocities. 


If Particle is in Target State: 


Determine Location (Voxel): 
Increment number of particles 
in that volume element. 


Figure 1 : Flowchart of Monte Carlo Model 




Figure 2: NH2 vith solar acceleration 


Table 1: Timescales used in Monte Carlo model 
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Dan Boice 

This is very interesting and important work. You mentioned that a future extension will be 
the inclusion of distributed sources. Gould you elaborate on how* you will do this and what 
your goals will be? 

Tony Ferro 

With this model, I believe, that a distributed source could be modeled by simply assuming 
a very large parent writh a mass similar to that of dust grains, and modifying the ejection 
velocity. This would act as another step in the model. 
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A QUASI -LINEAR THEORY TO EXPLAIH 
ION ACCELERATION H I THE DISTART 

COMXTAEY EHVHUMW. 


Raj 1 Sinha (University of New Mexico ), S. Peter Gary (Las-Alamoe National 
Lab) , Norman Roderick (University of New Mexico ) . 


Spacecraft observations at comet halley as well as computer simulations 
have demonstrated thar pitch -angle scattering of newborn cometary ions 
proceeds at a relatively fast rate, leading to relatively isotropic shell- 
like velocity distribution functions, Energization processes whereby shell 
distributions become more Maxwellian- like and a few ions are accelerated to 
high energies appear to proceed more slowly. This paper will describe our 
research on the latter, slower process, in which the scattering is assumed 
to be due to the resonant, growing magnetic fluctuations driven by the non- 
Maxwe Ilian nature of the ion distribution. 

We begin with an ion shell distribution which is isotropic in the wave 
frame and use a Fokker-Plank equation derived from quasilinear plasma theory 
to describe the broadening and energization of the cometary ion 
distribution. 
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MODELS OF THE SPIN STATE CHP THE COMET HALLEY NUCLEUS 

William H. Julian ~ f 

Department of Mathematical Sciences, New Mexico State University (f^ * 
Las Cruces, New Mexico 88003 USA 


Abstract 

Eight rotation-precession models of the comet Halley nucleus have been proposed by 
various authors: (1) Julian <1987), Sagdeev et <z/.(1989), and Abergel et al.( 1989); (2) Sekanina 
(1987); (3-5) Belton (1990); (6) Mdhlmann (1989); (7) Festou et a/.(1987); and (8) Wilhelm 
(1987) and Watanabe (1989). Julian (1990) has evaluated the eight models in relation to the 
constraints imposed by: (1) the observed long-axis directions at the Vega 1, Vega 2, and Giotto, 
encounters, (2) the ground-based emission periods harmonically related to 7.4 days (Belton 
1990), (3) the need for a two-day spin period in the analysis of the jet morphology, (4) the 
Smith et al.( 1987) constraint on the net iong-axis roll between Vega 2 and Giotto, (5) the 
resistance of the spin state of the nucleus to change due to the torque from the jets, and (6) the 
7.4-day repetition of the spatial orientation of the nucleus. 

Constraint 1. Sagdeev et al.( 1989) analyzed spacecraft-encounter images to determine 
the directions of the long axis of the comet Halley nucleus. They found a solution involving a 
2.2-day long-axis precession (of the long axis around the spin axis). Belton (1990) 
re-examined the motion of the nucleus with a 3.7-day long-axis precession and with the ends 
of the long axis reversed at Vega 1 and reflected in the plane of the sky as seen from Giotto. 
The eight models (Julian 1990) were constructed to satisfy either Sagdeev's solution (Models 
1, 2, 6, 7, and 8) or Belton's solution (Models 3, 4, and 5); all are marked Y (Yes) under 
constraint 1 in the Table. 

The inertial motion of the nucleus consists roughly of the long-axis precession 
combined with a roll of the nucleus around its long axis. The roll period is close to 7.4 days 
in Models 1, 2, 3, and 5; to 2.5 days in Model 4; to 5 days in Model 6; and to 15 days in 
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Models 7 and 8. Models 1, 3, and 8 are SAMs, while the rest are LAMs. In SAMs the roll 
motion is oscillatory, while for LAMs it is a complete rotation. 

Constraint 2 . Belton (1990) analyzed observations of the emission from Halley with the 
windowCLEAN algorithm to determine the dominant observed emission periods. He found 
that the dominant periods were constrained to be harmonically related to roughly 7.4 days (i.e, 
7.4, 3.7, 2.5, and 1.9). He ruled out a 2.2-day periodicity in the observed emission. This calls 
into question Models 1, 2, 6, 7, and 8 which were constructed with a 2.2-day long-axis 
precession. Models 3, 4, and 5 satisfy Belton's constraint and are accordingly marked Y under 
constraint 2 in the Table. 

Constraint 3 A two-day spin period derives from the analysis of jet morphology 
(Hoban et at. 1988; Keller and Thomas 1988; Rabinowitz 1988; Schulz and Schlosser 1989). 
Models 3 and 5 require a re-analysis of the jet morphology with a three-day spin period; the 
other models are marked Y under constraint 3 in the Table. 

Constraint 4. Smith et a/.(1987) compared the similar profiles of the nucleus that were 
seen from the Vega 2 and Giotto spacecraft. They deduced that the nucleus rolled an integral 
number of times between the two images which were taken 4.7 days apart. The models for 
which the residual net roll computed by Julian (1990) is less than 45° are marked Y under 
constraint 4 in the Table. Models 3, 4, 6, and 8 are favored; the others are not. 

Constraint 5. Constancy of the spin state is required by various analyses. Yeomans and 
Kiang (1981) have suggested that the spin state of the comet Halley nucleus has been constant 
for over 27 apparitions. Schleicher and Bus (1990) have determined that the emission periods 
seen in the 1910 apparition are consistent with a negligible change in spin state between that 
apparition and 1985/1986. All SAM models (1, 3, and 8) have a spin state that is readily 
changed by the torques from nine randomly-placed jets (Julian 1990); the other models are 
marked Y under constraint 5 in the Table. 

Constraint 6. Schleicher et al. (1990) have concluded that "the data demand that the 
nucleus returns to essentially the same orientation with respect to the sun each 7.4 ± 0.2 days. 


except for longer-time-scale seasonal change." Models 3, 4, and 5 are based on a 3.7-day 
long-axis precession and on roll periods of 7.4, 2.3, and 7.34 days, respectively. These periods 
are nearly harmonics of 7.4 days. Thus, the models return to nearly the same inertial 
orientation each 7.4 ± 0.2 days and are marked Y under constraint 6 in the Table. The other 
models do not, since 2.2 days is not a harmonic of 7.4 days. 


TABLE. CONSTRAINTS. 


Constraint 

Model 

1 

2 

3 

4 

5 

6 

1 

Y 


Y 




2 

Y 


Y 


Y 


3 

Y 

Y 


Y 


Y 

4 

Y 

Y 

Y 

Y 

Y 

Y 

5 

Y 

Y 



Y 

Y 

6 

Y 


Y 

Y 

Y 


7 

Y 


Y 


Y 


8 

Y 


Y 

Y 




Conclusions . All eight models satisfy constraint 1. If we require that the emission periods are 
harmonics of 7.4 days (constraint 2) then we are left with Models 3, 4, and 5. These are 
exactly the models that satisfy constraint 6. Model 3 is sensitive to inter-apparition spin state 
change, as are all the SAM models. Model 4 with a 3.7-day long axis precession and 2.5-day 
roll satisfies all six constraints. Its emission spectrum is weighted towards the fifth harmonic 
of 7.4 days, while that of Model 5 is weighted towards the third. Model 3 with 3.7-day long 
axis precession and 7.4-day roll satisfies constraints 1, 2, 5, and 6. If we accept Model 5, then 
we must re-examine the Smith et al.( 1987) long-axis roll constraint, and re-analyze the jet 
morphology with a three -day spin period. 


The roll and long-axis precession periods of Model 4 stand in a 2:3 resonance so the 
nucleus rolls three times around its long axis in the 7.4 days required for the long axis to 
process twice around the spin axis. For Model 5 the resonance is 2:1 so the nucleus rolls once 
around its long axis in the 7.4 days required for the long axis to process twice around the spin 
axis. The next step is to fit the observed ground-based emission waveform (Millis and 
Schleicher 1986; Schleicher et al. 1990) and the jet morphology (Hoban et al. 1988; Keller and 
Thomas 1988; Rabinowitz 1988; Schulz and Schlosser 1989). 
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THE CRAF MISSION AND EARTH-BASED COMET OBSERVATIONS 




Paul R. Weissman and Marcia Neugebauer 
Earth and Space Sciences Division 
Jet Propulsion Laboratory 
Pasadena, CA 91109 


The Comet Rendezvous Asteroid Flyby (CRAF) mission was approved for a New Start in 
the 1990 NASA appropriation, along with the Cassini (Saturn Orbiter/Titan Probe) mission, 
to be built simultaneously. Detailed design and planning of the CRAF spacecraft, mission, 
and scientific payload, are now underway. The CRAF mission will perform detailed studies 
of a typical short-period comet as it moves around its orbit, building on the information 
gained from the Giotto, Vega, and Suisei fast flybys of Comet Halley in 1986. 

CRAF will be launched in August, 1995 on a Titan I V/Centaur rocket, and will match orbits 
with short-period Comet Kopff in July, 2000, 850 days before perihelion at a heliocentric 
distance of 5.05 AU. After a series of close, slow flybys, the spacecraft will be placed in 
orbit around the cometary nucleus at an altitude of - 63 km, yielding better than 1 
meter/line-pair imaging resolution on the nucleus surface. One year after arrival, CRAF 
will fire an instrumented penetrator into the surface of the cometary nucleus to make 
detailed compositional and thermo-physical measurements over a period of about eight days. 
Instruments onboard the orbiter will collect and analyze cometary dust and gases, study the 
interaction of the cometaiy coma with the solar wind, and image and spectrally analyze both 
the nucleus surface and the coma, as the comet moves closer to the Sun and activity grows. 
Current plans call for the spacecraft to fly 50,000 km down the comet’s tail, shortly after 
perihelion passage in late 2002. On the way to Comet Kopff, CRAF will fly by the 88-km 
diameter, C-type, main belt asteroid 449 Hamburga in January, 1998, at a closest approach 
distance of about 3,500 km, adding to our knowledge of these primitive bodies which are 
believed to be the source of many meteorites on Earth. In addition, CRAF will perform a 
gravity assist flyby of the Earth in July, 1997. CRAF is a joint project between NASA and 
the Federal Republic of Germany. 

Earth-based observations of the target comet, either by ground-based facilities or by Earth- 
orbiting satellites, can provide essential information for CRAF mission planning, and are 
needed in support of the rendezvous mission. Estimates of the nucleus size, albedo, shape, 
rotation period, and rotation pole orientation are needed for preliminary mission planning, 
resource allocation, and sequence design. Measurement of the target asteroid’s rotation 
period and shape, as well as a detailed spectral reflectance curve, are also needed. 
Astrometric observations of both the target comet and asteroid can result in improved 
ephemerides, helping to optimize trajectory design. Ground-based observations of the 
comet during the rendezvous will allow comparison with in situ measurements by the 
spacecraft, thus improving the interpretation of ground-based observations of other short- 
period comets. 
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Opportunities also exist for active participation in CRAF data analysis by scientists not 
currently associated with the mission. A program of Participating Scientists is planned which 
will supplement currently selected science teams during flight. Announcements of 
Opportunity will be issued by NASA approximately one year prior to the Hamburga flyby, 
and later approximately one year prior to rendezvous with Comet Kopff 

The attached table provides orbital data on the future perihelion passages of Comet Kopff, 
as well as that for two other short-period comets which are prime backup targets for the 
CRAF mission, Comets Tempel 2 and Wild 2. The CRAF science office is interested in 
receiving information on any observations of these three comets. Observations of interest 
include measurements of gas production rates, dust production rates, and visual magnitudes 
at any point in the comets' orbits, and observations of the bare cometary nuclei at large 
solar distance. Evidence of cometaiy activity at large heliocentric distance, beyond 3 AU, 
would be particularly valuable. Please send any information, preprints, reprints, abstracts, 
etc. to the authors above. 

Comet Kopff is particularly well placed for Earth-based observations during its 1996 
perihelion passage, as is Comet Wild 2 in 1997, and Comet Tempel 2 in 1999. Comet Wild 
2 is well placed for post-perihelion observations in 1990, 

A more complete description of the CRAF mission can be found in: Neugebauer, M., and 
Weissman, P. R., CRAF mission, EOS 70, 633, 1989. 


CRAF Target Comets* 


Comet 

T 

Peri. 

Node 

Incl. 

q 

e 

Kopff 

1990 1 20.3725 

162.8291 

120.2891 

4.7205 

1.585125 

0.542991 

1996 7 2.1840 

162.7631 

120.2823 

4.7245 

1.579523 

0.544079 


2002 12 12.0646 

162.6860 

120.2977 

4.7222 

1.583575 

0.543302 

Wild 2 

1990 12 16.9253 

41.5729 

135.5763 

3.2494 

1.578058 

0.540983 


1997 5 6.6533 

41.6992 

135.5280 

3.2481 

1.582620 

0.540225 


2003 9 25.9857 

41.6832 

135.5130 

3.2457 

1.590364 

0.538794 

Tempel 2 

1994 3 16.8162 

194.8570 

117.5760 

11.9783 

1.483536 

0.522447 


1999 9 8.4247 

194.9972 

117.5388 

11.9798 

1.481679 

0.522818 


2005 2 15.0423 

195.5328 

117.1761 

12.0208 

1.426862 

0.535387 


’Orbital data from: Yeomans, D. K., and Wimberly, R. N., Cometaiy apparitions: 1990 - 
2010, in Comets in the Post-Halley Era, eds. R. L Newbum, J. Rahe, and M. 
Neugebauer, Kluwer, Amsterdam, in press, 1990. 
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CONFERENCE SUMMARY 



W. F. Huebner 
Southwest Research Institute 
San Antonio, TX 78228*0510 


We had hoped that Jurgen Rahe would prepare and present the summary of this work- 
shop. Unfortunately he is unable to attend. Therefore, I will try to summarize thie meeting 
in his stead. 

We have heard reports on the observations of four comets: Brorsen-Metcalf (1989o), 
Okazaki-Levy-Rudenko (1989r), Aarsetb-Brewington (1989al), and Austin (1989cl). Al- 
though only preliminary results are available at this time, the path of research that has 
been initiated is clear: In many regards it follows closely the directions established by the 
new Comet P/Ilalley findings. 

The relationship of minor species with each other and with possible parents as well as 
with dust are being pursued in a number of investigations. Of particular interest are the 
abundance ratios of CH 4 to CO and NH 3 to Nj, because they give information about the 
place of origin of comets. These determinations are very difficult because there are no direct 
observations of CH< or N 3 and NH 3 may be detectable only in the radio range and apparently 
no such observations were made. Thus the determination of these ratios is indirect* Closely 
related to this line of research and to the place of origin of comets are the isotope ratios of 
12 C/ 13 G. Whereas the 13 C isotope showed a relatively strong spectrum in Comet Halley, this 
is apparently not the case with the comets under investigation here. 

The relative constancy from comet to comet of some, presumable chemically unrelated 
minor species may also yield information about the origin of comets. But first it will be 
necessary to establish whether these species originate in the icy or the dusty component of 
the nucleus. The fact that some have an almost constant ratio with respect to OH seems to 
indicate that their abundances are closely correlated to that of water, for example, they may 
be trapped in water ice to their point of saturation. Establishing a correlation of the minor 
species to the dust is a much more difficult problem that still needs further investigation. 

Other species, such as S 3 , appear to vary much more widely. In fact, only upper limits 
have been established for S 2 in the comets under investigation here. Unidentified lines and the 
new species CH 3 OH and H 2 S, reported for the first time at this workshop, give new impetus 
to comet research. 

Observations of dust in the IR and the correlation of the dust with visual brightness 
of a comet and its chemical composition require close cooperation between various observers. 
As reported at this workshop, a trend is being established for simultaneous, multi-instrument 
observations by closely collaborating observing teams. In this regard simultaneous radio and 
optical observations are also emerging as a new vehicle. 

Similar to multi-instrument observations, the time may be ripe for closer collaboration 
between observing teams and modelers. The Haser model, although very useful for short-lived 
species and for quick analyses, cannot be used for species involving several, nearly equally 
effective parents, as may occur in parallel photolytic and chemical reactions. Similarly, the 
dust size distribution as a function of cometocentric distance needs to be analyzed in closer 
collaboration between observers and modelers. 

This has been a very instructive meeting -with many very open discussions. We had 
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international participation with two poster papers (Kozuka et al. and Watanabe et al.) from 
our Japanese colleagues and two papers (Bockelee-Morvan et al. and Colom et al.) from our 
French colleagues, presented by Bockelee-Morvan who made a special trip to present these 
papers at our workshop. I am very grateful to Dr. Bockelee-Morvan for making this effort and 
want to thank all of you for attending and participating. I also want to thank Ms. Staglik from 
the SwRI for taking care of the many details in arranging this workshop and Mr. Sommerfeld 
from the University of Albuquerque for operating the projection equipment. Finally, I would 
like to thank my colleagues, Jurgen Rahe, Peter Wehinger, and Ichishiro Konno for their help 
in organizing this workshop and the NASA Planetary Atmospheres Program for supporting 


